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'^Second-echelon  interdiction  is  closely  tied  to  the 
close-in  battle.  In  reality,  it  is  part  of  the  integrated 
battle  which  says  that  all  means  are  used  to  fight  the  battle 
at  all  distances.  The  objective  of  this  research  effort  was 
to  examine  what  effect  the  combination  of  air  sorties  and 
Army  missiles  would  have  against  a  motorized  rif  division 
(MRD)  moving  on  a  hypothesized  three-route  marcl  A  portion 
of  a  selected  East  German  road  network  was  chos  as  the 
hypothesized  area  through  which  an  MRD  would  be  ^.veling 
on  its  way  to  a  release  point. 

A  model  of  a  MRD  broken  down  into  44  units  was  con¬ 
structed  using  the  SLAM  computer  simulation  language.  Units 
were  broken  down  to  battalion  size,  with  a  few  exceptions. 

A  flight  of  two  aircraft  is  used  to  interdict  the  road  on 
which  convoys  are  traveling.  Once  a  convoy  is  blocked  by 
the  air  strike,  a  retargeting  is  conducted  against  the  convoy 
using  a  generic  Army  missile  loaded  with  wide-area  anti¬ 
armor  munitions  (WAAM) .  Both  the  sortie  interdiction  and 
missile  attack  were  modeled  explicitly  using  Monte  Carlo 


simulation . 


STOPEM:  A  SIMULATION  INTERDICTION  MODEL 


OF  A  MOTORIZED  RIFLE  DIVISION 


I  Introducticn 


A  topic  of  general  interest  that  has  sparked  much 
attention  recently  is  the  outcome  of  a  possible  general  war 
in  Europe.  Since  the  United  States  maintains  its  largest 
military  contingent  outside  its  continental  borders  as  part 
of  its  commitment  to  NATO,  the  possibility  exists  for  a  mili¬ 
tary  confrontation  between  NATO  and  the  Warsaw  Pact  that 
would  involve  the  two  superpowers.  While  a  general  war 
brings  visions  of  worldwide  destruction  involving  the  nuclear 
arsenals  of  both  superpowers,  one  scenario  woulo  keep  the 
conflict  at  a  conventional  level  within  NATO. 

The  system  of  war,  with  its  many  interactions,  has 
three  escalatory  levels  above  conventional  war:  (1)  chemical- 
biological  (C-B),  (2)  tactical  nuclear,  and  (3)  strategic 
nuclear.  While  no  one  can  predict  when  and  how  C-B  or 
nuclear  weapons  would  be  used,  the  rational  person  would 
argue  that  use  of  these  weapons  would  probably  be  preempted 
by  use  of  conventional  weapons  to  decide  a  conflict  or  con¬ 
ventional  superiority  could  force  negotiations.  As  part  of 
this  same  rationale,  the  fact  that  both  sides  possess  one 
or  both  capabilities  would  be  a  deterrent  for  either  side 
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to  employ  them  first.  Understanding  the  Warsaw  Pact's  con¬ 
ventional  forces'  vulnerabi 1  it ies  and  weaknesses  within 


this  system  context  is  essentia]  to  stopping  them  convention 
ally  so  that  the  next  threshold  of  war  is  not  breached. 
Although  this  is  one  possible  scenario,  another  one  in  the 
literature,  "A  Selective  Nuclear  Policy  Strategy  in  Europe: 

A  Selective  Targeting  Doctrine?'',  mentions  the  use  of  Warsaw 
Pact  selective  targeting  through  use  of  chemical  and/or 
nuclear  munitions  to  facilitate  capturing  Europe  relatively 
intact,  since  their  conventional  forces  are  equipped  to 
operate  in  a  chemical  or  nuclear  environment  (Ref  19). 

While  little  is  openly  published  about  Soviet  doc¬ 
trine,  The  Offense  (A  Soviet  View)  discusses  employment  of 
nuclear  weapons  as  an  extension  of  firepower  combined  w-ith 
maneuver  (Ref  53 : v , vii-viii ) .  Because  of  this  lack  of 
information,  little  is  known  about  Warsaw  Pact  doctrine. 

On  the  other  hand,  more  is  being  written  in  U.S.  military 
journals  about  targeting  and  employment  of  nuclear  weapons 
to  interdict  enemy  forces  while  they  are  deep  in  their 
territory  (Ref  30:2-6,  33:2,  57:35).  While  the  Soviet's 
military  force  is  capable  of  overwhelming  conventionally 
any  adversary,  with  the  added  capability  for  employment  and 
operation  in  a  C-B  and  nuclear  environment,  the  U.S.  has 
responded  to  this  threat  with  increased  technological 
enhancements  in  equipment  and  weapons  and  with  more  of  a 
willingness  to  possibly  use  nuclear  weapons  to  offset  the 
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numerical  oviet  advantage.  The  clear  danger  exists  that 
if  a  conflict  starts  conventionally,  the  side  that  is  frus¬ 
trated  first  or  feels  that  it  is  losing  the  initiative 
might  escalate  into  the  next  level  of  war.  Once  the  war 
escalates  into  the  nuclear  environment,  the  costs  to  both 
sides  might  far  outweigh  the  gains. 

Within  the  last  10  years,  the  U.S.  Army  has  examined 
how  to  fight  the  next  war.  This  examination  has  stemmed 
from  the  misunderstanding  of  Field  Manual  (FM)  100-5,  pub¬ 
lished  in  1976,  to  its  rewriting  in  1982.  The  many  articles 
published  in  military  journals  have  had  a  great  influence 
on  the  revision  of  this  key  FM.  The  study  of  adversaries' 
weaknesses  and  strengths,  and  the  questioning  of  U.S.  doc¬ 
trine,  must  continue  if  the  Army  and  Air  Force  expect  to  be 
ready  to  fight  and  win  any  future  conflict.  Because  of  the 
introduction  of  many  new  weapon  systems  now  and  in  the  near 
future  (see  Figure  1.1),  Army  commanders  will  have  at  their 
disposal  the  means  to  see  and  attack  deep  on  the  battlefield. 
Systems  such  as  the  corps  support  weapon  system  ( CSWS ) , 
ground-launched  cruise  missile  (GLCM),  multiple  launch 
rocket  system  (MLRS),  and  family  of  scatterable  mines  (FASCAM) 
give  the  commander  a  capability  to  attack  near  and  far  tar¬ 
gets.  Other  systems  such  as  all-source  analysis  system 
(ASAS),  tactical  satellite  (TACSAT),  stand-off  target 
acquisition  system  (SOTAS),  and  tactical  fire  direction 
(TACFIRE)  provide  the  detection  means  and  the  integration 


capability  for  employment  of  all  weapon  systems.  Having 
the  weapons  without  the  means  to  effectively  employ  them 
reduces  their  worth.  Combined,  all  systems  in  Figure  1.1 
present  a  formidable  capability.  However,  owing  to  war's 
complexities,  a  need  exists  to  examine  how  these  numerous 
systems  will  interact  and  how  they  will  damage  or  delay 
attacking  forces. 

Figure  1.2  is  a  graphical  representation  of  the 
Warsaw  Pact  second-echelon  threat.  As  the  figure  shows,  the 
enemy  units  are  deployed  to  achieve  depth  and  dispersion. 

Not  only  is  the  bunching  of  units  prevented,  but  the  majority 
of  the  forces  are  uncommitted  to  exploit  success  or  reinforce 
other  enemy  formations.  The  echelonment  of  forces  is  an 
exploitable  weakness  upon  which  NATO  must  capitalize  if  the 
likelihood  of  victory  on  the  battlefield  is  to  be  increased. 
While  "echelonment  is  neither  axiomatic  nor  uniform,"  this 
thesis  assumes  an  attack  on  NATO  w’ould  use  the  idea  of 
echelonment  (Ref  4:40).  General  Starry,  former  Training  and 
Doctrine  Commander  (TRADOC),  keyed  on  the  importance  of 
echelonment  of  forces.  He  stated  that  the  Warsaw  Pact 
keeping  a  significant  portion  of  its  forces  uncommitted 
allowed  it  to  retain  the  advantage  of  the  initiative  to 
commit  at  its  choosing  either  to  reinforce  success  or  to 
bypass  friendly  forces.  Furthermore,  to  upset  the  Warsaw 
Pact  plan,  this  initiative  must  be  taken  away  from  the  Warsaw 
Pact  and  then  retained  by  NATO  to  gain  victory  (Ref  57:34). 
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Fig  1.1.  A  Substantial  Step  Toward  Future 
Capabilities  (Ref  57:33) 


Fig  1.2.  The  Second-Echelon  Threat 
(Ref  57:35) 
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A  simulation  model  of  second-echelon  interdiction 
would  give  some  insights  about  interdiction  capability  as 
well  as  providing  possible  answers  to  how  long  enemy  forces 
could  be  kept  out  of  a  sector  and  what  their  combat  power 
would  be  upon  their  arrival  at  the  point  of  commitment.  A 
simulation  model  of  interdiction  of  second-echelon  forces 
would  give  policy  makers  an  evaluation  tool  to  determine 
where  possible  conventional  interdiction  shortfalls  might 
exist.  An  analytic  approach  to  addressing  this  problem  was 
not  considered  because  the  system  being  examined  is  large, 
complex,  and  stochastic.  Thus,  simulation  can  best  be 
employed  to  provide  insights  to  this  particular  problem. 
Having  a  credible,  conventional  interdiction  capability 
should  lower  the  U.S.  nuclear  response.  However,  the  uncer¬ 
tainty  of  the  Warsaw  Pact  initial  response,  or  its  response 
if  thwarted  in  any  attack,  still  remains. 

Problem  Statement 

What  is  the  structure  of  the  interdiction  subsystem 
contained  within  the  overall  system  of  conventional  war, 
how  can  this  particular  subsystem  be  captured  in  a  dynamic 
model,  and  how  can  this  model  be  used  to  evaluate  specific 
policies?  More  specifically: 

1.  What  are  the  significant  relationships  in  Corps- 
directed  aircraft  sortie  interdiction  of  terrain-restricted 
roads  and  the  effects  on  second-echelon  forces  moving 
through  the  road  network  interdicted? 


2.  How  would  a  relook  capability  that  enables 


retargeting  of  the  interdicted  point  with  missiles  loaded 
with  wide  area  anti-armor  munitions  OVAAM)  reduce  the  combat 
power  of  delayed  units  or  increase  their  delay  time? 

3.  How  can  these  relationships  be  incorporated  into 
a  model  that  will  enable  military  planners  to  evaluate 
capabilities  and  effectiveness  of  selected  conventional 
munitions  against  second-echelon  threats? 

Objectives 

The  primary  objective  of  this  research  is  to  provide 
a  validated  model  of  interdicting  a  selected  road  network 
in  East  Germany  to  gain  insights  into  how  long  a  Warsaw  Pact 
motorized  rifle  division  (MRD)  could  be  delayed  from  being 
committed  into  a  corps  sector  and  what  its  combat  power 
would  be  upon  reaching  its  destination.  Intermediate 
objectives  are: 

1.  Develop  a  basic  model  containing  a  typical  MRD 
broken  into  battalion-size  units  making  a  100  kilometer  (Km) 
road  march  on  three  separate  routes. 

2.  Add  an  interdiction  scheme  that  shows  how  these 
effects  alter  the  MRD’s  arrival  and  combat  power,  where 
combat  power  refers  to  its  main  fighting  systems  (armor, 
infantry,  artillery). 

3.  Verify  and  validate  the  model. 

4.  Use  the  model  to  evaluate  effects  of  aircraft  and 


missile  sortie  allocation  on  the  MRD. 


System  Structure 

This  research  is  directed  at  understanding  and 
modeling  the  flow  of  convoys  through  a  road  network  and  how 
to  impede  their  flow  so  that  their  arrival  into  a  corps 
sector  would  either  be  delayed  for  a  short  time  with  high 
losses,  be  delayed  for  a  long  time  with  low  losses,  be 
delayed  for  a  long  time  with  high  losses,  or  not  be  affected 
at  all.  The  ideas  above  represent  four  possible  enemy  force 
statuses  of  combat  capability.  For  example,  if  a  force 
could  not  be  delayed  out  of  sector  for  as  long  as  a  pre¬ 
specified  time,  say  12  hours,  then  a  possible  acceptable 
alternative  might  be  to  keep  it  out  less  time,  but  with  a 
reduced  percentage  of  combat  vehicles,  such  as  tanks,  infan¬ 
try,  and  artillery,  by  diverting  more  assets  to  interdiction 
To  achieve  such  goals  requires  that  all  resources  be  effec¬ 
tively  integrated  on  any  future  battlefield  so  that  the 
likelihood  of  success  is  increased. 


Background 

To  set  the  problem  in  perspective,  a  description  of 
the  integrated  battlefield  is  given  first.  A  corps  sector 
is  described  and  an  explanation  of  how  the  Corps  commander 
views  the  battlefield  is  given.  Also,  the  conduct  of  the 
integrated  battle  is  described  along  with  how  deep  attack 
affects  enemy  force  arrivals  at  the  forward  line  of  own 
troops  (FLOT).  Next,  an  examination  is  made  of  the  last 
decade's  questions  aoout  Army  doctrine.  Specifically,  the 


debate  is  traced  that  questions  how  the-  Army  will  fight  the 
next  war  and  what  the  implications  are  of  the  active- 
defense. 

System  Description.  Figure  1.3  is  a  causal  loop 
diagram  that  depicts  the  complex  interactions  of  deep  attack. 
Such  diagrams  are  used  in  system  analysis  as  tools  to 
graphically  represent  and  to  aid  in  visualization  of  the 
system  structure  and  key  relationships.  Feedback  between 
or  among  relationships  denotes  that  one  relationship  can 
influence  another.  However,  even  when  relationships  do 
exist,  feedback  may  not  exist.  Connections  between  variables 
are  depicted  as  solid  lines  with  an  arrow-tipped  end  that 
either  has  a  positive  or  a  negative  sign.  A  positive  sign 
indicates  that  an  increase  in  one  variable  will  result  in  an 
increase  of  another  variable  in  the  same  loop.  A  minus  sign 
denotes  that  increasing  one  variable  will  decrease  the  other 
one  at  the  arrow  tip.  By  multiplying  all  the  signs  within  a 
loop,  a  net  sign  is  obtained.  A  positive  net  sign  indicates 
a  positive  or  reinforcing  loop.  An  outside  type  of  control 
prevents  a  continued  increase  of  relationships  when  acted 
upon  by  some  external  influence.  The  relationships  within 
the  loop  continue  increasing  unless  restrained  by  this 
external  factor.  On  the  other  hand,  a  negative  or  goal¬ 
seeking  loop  tends  toward  equilibrium  or  a  balance  when 
acted  upon  by  an  outside  element.  The  net  effect  of  this 
process  is  that  causal  loop  diagrams  aid  the  models  during 
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the  development  phase  by  assisting  in  conceptualizing  and 
formulating  key  relationships.  Thus,  a  better  system  under¬ 
standing  is  the  result. 

Figure  1.3  shows  part  of  the  complexities  associated 
with  deep  attack.  As  the  figure  shows,  interdiction  of  red 
forces  decreases  blue's  interdiction  resources,  but 
increases  red  equipment  destruction,  red  delay,  red  target 
acquisition,  diversion  of  red  assets,  red  vulnerability,  red 
attrition,  and  blue's  reconstitution  time.  As  the  figure 
also  shows,  these  same  variables  affect  other  variables 
within  the  loop.  For  example,  attrition  can  have  an  affect 
on  destruction,  delay,  or  disruption.  But,  destruction, 
delay,  or  disruption  also  affect  red  attrition  because  as 
convoy  elements  are  destroyed,  delayed,  or  disrupted  on 
their  routes  of  march,  they  present  better  targets.  This 
bunching  of  vehicles  on  the  roadway  increases  red's  vulner¬ 
ability  as  well  as  facilitating  acquisition.  One  fallacy 
with  stating  that  target  acquisition  is  always  improved  is 
that  terrain  masking  may  make  it  more  difficult  to  detect 
targets.  So  in  actuality,  target  acquisition  and  other 
variables  could  be  decreased  due  to  an  inability  to  see 
targets.  This  is  a  function  of  such  variables  as  terrain, 
degree  of  cover,  time  of  the  year,  weather,  and  detection 
equipment . 

A  portion  of  a  corps  sector  was  selected  for  this 
problem  as  being  the  area  where  an  MRD  would  penetrate  the 


sector  and  become  the  concern  of  the  commander.  Also,  this 
same  sector  is  part  of  an  overall  area  of  interest  being 
examined  to  provide  insights  to  the  CSWS  office,  Fort  Sill, 
Oklahoma.  Figure  1.4  shows  the  whole  corps  sector  organiza¬ 
tion.  A  brief  explanation  follows  of  the  key  elements  on 
the  figure  as  well  as  tying  them  to  the  concept  of  the 
extended  battlefield. 

Briefly  stated,  the  extended  battlefield  idea  is  a 
more  descriptive  term  that  captures  how  to  view  the  battle¬ 
field  in  terms  of  time,  space,  enemy  forces,  friendly  forces 
and  weapons  employment.  Enemy  forces  are  to  be  engaged  whil 
not  in  contact  to  frustrate  their  command  and  control,  and 
to  strip  away  their  initiative.  Secondly,  all  current 
actions  are  interrelated  in  time  in  that  plans  for  the  close 
in  battle  as  well  as  attack  of  follow-on  forces  are  tied  to 
winning  at  the  FLOT .  Finally,  the  integrated  assets  of 
"higher  level  Army  and  sister  services"  are  carefully 
employed  in  interdicting  the  battlefield  (Ref  57:32).  The 
key  point  is  to  carefully  employ  all  assets  and  resources  so 
that  interdicting  deep  targets  is  tied  to  the  close-in 
battle . 

As  enemy  forces  penetrate  the  corps  sector  and  are 
located  about  96  hours  out  from  the  FLOT,  they  are  within 
the  corps  area  of  interest.  Within  this  area  the  commander 
must  monitor  enemy  movement  so  that  he  can  determine  what 
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Fig  1.4.  Organization  of  the  Defense  -  Corps 


forces  could  affect  his  future  operations.  The  decision  to 
interdict  is  made  here. 

Closer  in  exists  the  area  of  influence  where  the 
commander  is  also  capable  of  acquiring  the  fighting  enemy 
units  with  assets  organic  to  or  in  support  of  his  command. 

For  instance,  a  division  may  be  maneuvering  units  against 
the  second-echelon  regiments  of  the  first-echelon  division. 

In  order  to  increase  the  chances  of  success  against  the  enemy 
forces  and  prevent  their  being  reinforced  by  second-echelon 
divisions,  the  second-echelon  divisions  would  have  to  be 
destroyed,  delayed,  or  disrupted  for  a  prespecified  period 
of  time  and/cr  decreased  level  of  combat  power.  This  is  an 
area  where  intelligence  estimates  and  advance  planning  would 
pay  off.  Tables  1.1  and  1.2  show  areas  of  influence  and 
interest  as  functions  of  level  of  command  and  time  distance 
from  the  FLOT  in  hours. 

TABLE  1.1 
Areas  of  Influence 

Distance  from  FLOT 

Level  of  Command  (time  in  hours) 


Brigade 

Division 

Corps 


0  to  12 
0  to  24 
0  to  72 


TABLE  1.2 


Areas  of  Interest 


Level  of  Command 


Distance  from  FLOT 
(time  in  hours) 


Brigade 

Division 

Corps 


0  to  24 
0  to  72 
0  to  96 


(Ref  16:53) 


The  remainder  of  Figure  1.4  depicts  the  location  of 
friendly  forces.  The  point  where  the  forwardmost  friendly 
elements,  the  covering  force,  are  located  is  defined  as  the 
FLOT.  From  the  FLOT  back  to  the  forward  edge  of  the  battle 
area  (FEBA)  is  designated  the  covering  force  area  (CFA)  where 
the  corps  covering  force  (CF)  will  provide  early  warning  to 
the  forces  in  the  main  battle  area  (MBA)  about  the  enemy's 
intent  and  direction  of  main  attack,  "to  develop  the  situa¬ 
tion,  and  to  delay  or  defeat  the  enemy's  leading  fighting 
force"  (Ref  30:10-8).  As  the  battle  develops  in  the  CFA, 
commanders  make  final  preparations  in  the  MBA.  Also,  as  the 
CF  fights  back  to  the  MBA,  the  CF  hands  the  enemy  over  to 
the  forces  in  the  MBA,  where  the  majority  of  the  fighting  is 
anticipated  as  taking  place.  Consequently,  it  is  in  the 
MBA  where  the  majority  of  friendly  units  are  located. 

Finally,  the  rear  area  is  to  the  immediate  rear  of  the  rear 
boundary  of  the  MBA,  where  command  and  control  elements, 


reserve,  and  support  units  are1  located.  This  quick,  simpli¬ 
fied  description  of  how  a  corps  sector  is  organized,  is  more 
complicated  than  explained.  However,  the  points  covered  are 
the  highlights. 

As  Figure  1.5  shows,  the  advancing  enemy  force  moving 
in  is  detected  and  engaged  in  the  area  of  interest.  Within 
this  area,  the  enemy  forces  are  subjected  to  attack  by  all 
available  resources.  The  result  of  continued  interdiction 
all  along  its  movement  path  toward  the  FLOT  combined  with 
maneuver  create  the  result  in  Figure  1.6.  Finally, 

Figure  1.7  shows  the  reconstituted  FEBA,  after  destruction 
of  the  first-echelon  forces.  The  reconstituted  FEBA  is 
better  prepared  to  engage  the  new  units  entering  or  that 
have  entered  corps  sector. 


Fig  1.5.  The  Deep  Battle  (Ref  57:39) 


A  plot  of  this  effect  of  interdiction  is  shown  in 
Figure  1.8.  The  upper  portion  of  the  plot  of  enemy  front¬ 
line  strength  against  time  shows  that  no  interdiction  allows 
the  enemy  forces  to  maintain  a  fairly  constant  density  of 
forces  at  the  FLOT.  This  translates  to  more  combat  power 
to  maintain  the  initiative  and  to  overwhelm  the  friendly 
force.  The  bottom  curve  shows  the  effect  of  interdiction 
and  attack.  As  interdiction  destroys,  delays,  or  disrupts 
enemy  forces,  window's  for  action  are  created  where  a  time 
period  exists  that  f riendly-to-enemy  force  rat ios  are  favor¬ 
able  to  attack.  It  is  within  these  windows  created  by  inter 
dieting  that  maneuver  and  firepower  will  produce  destruction 
of  enemy  forces.  This  curve  shews  why  it  is  absolutely 
necessary  to  destroy,  delay,  or  disrupt  any  enemy  forces 
that  could  interrupt  the  maneuvering  friendly  units  before 
they  get  to  finish  their  destruction  mission.  The  resul+s 
from  the  figure  were  generated  from  simulation  comparisons 
conducted  by  the  Army's  Field  Artillery  School,  Fort  Sill, 
Oklahoma,  of  1980  European  corps  battles  (Ref  57:42). 


Fig  1.8.  Why  Deep  Attack?  (Ref  57:45) 


Doctrinal  Refinement.  Two  discernible  periods  of 
Army  doctrine  refinement  can  be  traced  from  the  writing  of 
FM  100-5  in  1976,  to  the  time  between  its  initial  publica¬ 
tion  and  its  republication  in  1982.  The  original  FM  version 
was  "notably  flexible  and  deliberately  nonrestrictive" 

(Ref  60:3).  The  active  defense  called  for  lateral  movement 
along  the  FEBA  to  achieve  concentration  of  forces  at  the 
threatened  point.  Depth  on  the  battlefield  was  achieved  by 
a  series  of  delaying  actions  conducted  all  the  way  back  to 
the  rear  boundary.  These  actions  were  the  maneuvers  asso¬ 
ciated  with  the  active  defense.  Offensive  action  was  inter¬ 
preted  to  be  practically  nonexistent.  At  the  rear  boundary, 
the  defense  ended  up  being  linear.  "It  was  this  rigidly 
limited  form  of  defense  which  attracted  most  of  the  attention 
of  critics  and  prompted  the  long  series  of  debates  which 
still  goes  on."  (Ref  60:3,4)  Colonel  Tate  and  Lieutenant 
Colonel  Holder  in  their  article,  "New  Doctrine  for  the 
Defense,"  went  on  to  say  that  the  doctrine  in  the  FM  100-5 
was  interpreted  to  be  dogmatic,  when  in  fact  it  was  supposed 
to  be  flexible  and  nonrestrictive.  Hence,  the  ensuing  debate 
and  the  eventual  rewriting  of  this  manual . 

The  author's  perception  is  that  a  critical  transition 
occurred  with  an  enlightening  article  written  by  Colonel 
Wagner,  "Active  Defense  and  All  That."  Colonel  Wagner  gave 
his  astute  conceptualization  of  how  the  active  defense  was 
to  be  fought.  Using  his  command,  the  11th  Armored  Cavalry 


Regiment,  as  an  example,  he  reintroduced  oliense  into  the 
active  defense.  Through  his  example,  he  described  his 
command  performing  a  covering  force  mission  using  all  avail¬ 
able  assets,  terrain,  and  maneuver  to  defeat  an  attacking 
force  (Ref  63).  The  scenario  and  tactics  he  described 
represent  the  essence  of  the  active  defense.  This  particular 
article  did  more  to  crystalize  the  flexibility  and  intent  of 
active  defense  than  any  other  article  up  to  that  point. 

The  article  by  General  Starry,  former  TRADOC  com¬ 
mander,  entitled,  ''Extending  the  Battlefield,”  was  a  bench¬ 
mark  in  further  clarifying  the  integration  of  all  assets  and 
thought  to  the  complete  depth  of  the  battlefield.  Commanders 
and  staff  had  to  view  the  battlefield  as  closely  interrelated 
parts  of  the  close-in  battle  and  deep  battle.  This  careful 
integration  of  all  assets  in  fighting  the  extended  battle¬ 
field  reminded  all  commanders  that  success  in  the  close-in 
battle  was  tied  to  the  deep  one  (Ref  57). 

Other  journal  articles  before  and  since  the  Wagner 
and  Starry  articles  also  deserve  some  brief  comment  to 
reflect  the  breadth  and  intensity  of  attempting  to  clarify 
doctrine  and,  in  some  cases,  introducing  innovative  ideas. 
Articles  such  as  Brit t ingham ' s ,  "Use  the  Lightning,"  describe 
an  expanded  role  for  attack  helicopters  to  capitalize  on 
their  potency  as  very  mobile  anti-armor  forces.  The  mobility 
and  flexibility  in  this  weapon  system  gives  the  commander  an 
excellent  counter  penetration  force.  Other  ideas  include 


doctrinal  questions  previously  discussed  as  well  as  articles 
by  some  Air  Force  authors  from  defining  the  mission  of  battle¬ 
field  air  interdiction  (BAI)  to  the  effects  of  European 
weather  on  round-the-clock  air  operations.  In  True's 
article,  "The  Tourniquet  and  The  Hammer,"  the  author  suggests 
in  his  title  the  analogy  of  deep  interdiction  using  two  his¬ 
torical  raids  on  the  Romanian  oilfields  at  Ploesti  (Ref  62). 
Frizzo  writes  about  some  unique  ideas  in  the  rapid  deploy¬ 
ment  of  U.S. -based  forces  to  Europe  being  used  as  infantry 
reinforced  with  an  anti-tank  capability.  These  forces 
represent  mechanized  units  or  light  infantry  units  deployed 
without  heavy  equipment  to  be  employed  against  the  second 
echelon  forces  (Ref  22).  Porreca  advocated  a  rethinking  of 
tactics  to  use  the  indirect  approach  as  espoused  by 
B.  H.  Liddell  Hart.  One  of  his  major  points  is  that  military 
school  students  have  become  so  enamored  of  the  most  current 
terminology  to  describe  current  doctrine  that  they  are  just 
parroting  the  party  line  without  offering  any  new  insights 
(Ref  48).  An  examination  of  the  bibliography  gives  an  indi¬ 
cation  of  some  of  the  professional  thought  that  has  served 
to  enlighten,  to  question,  to  analyze,  and  to  present  ideas 
that  are  relevant  to  this  crucial  issue.  This  brief  excur¬ 
sion  plus  the  synopsis  of  the  following  studies  will  serve 
to  illustrate  the  intensity  and  concern  being  shown  by  mili¬ 
tary  professionals  and  others  in  both  branches  of  the  service 
and  outside  agencies. 
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Air-Land  Forces  Application  (ALFA)  Agency,  head¬ 
quartered  at  Langley  Air  Force  Base,  Virginia,  is  a  joint 
service  organization  examining  the  problems  associated  with 
conduct  of  the  deep  battle.  Since  the  Air  Force  possesses 
more  assets  to  conduct  deep  interdiction,  to  include  second 
echelon  forces,  both  the  Army  and  Air  Force  have  had  to 
reexamine  the  doctrine  of  interdiction  beyond  the  FLOT  in  the 
Army's  defined  areas  of  interest  and  influence.  Tradition¬ 
ally,  the  Air  Force  had  had  almost  exclusive  say  in  what 
was  interdicted  beyond  a  certain  point  in  the  battlefield. 
But,  now  Army  commanders  want  a  bigger  input  in  influencing 
action  in  what  has  been  traditionally  an  Air  Force  domain 
because  affecting  the  arrival  of  echeloned  forces  has  a 
direct  impact  on  maneuver  schemes  in  the  close-in  battle. 
Since  both  services  are  working  toward  the  same  end,  it  is 
desirable  that  better  cooperation  and  compatible  doctrine 
exist  between  the  services  in  order  to  purposefully  employ 
limited  assets  to  get  the  greatest  effect. 

Two  Army  agencies  that  are  also  examining  the  problem 
of  second  echelon  interdiction  are  located  at  Fort  Sill, 
Oklahoma,  and  Fort  Leavenworth,  Kansas.  Fort  Sill  is  the 
home  of  the  Corps  Support  Weapon  Systems  (CSWS),  a  Department 
of  the  Army  task  force.  The  group  there  is  working  with 
Sandia  National  Laboratory  and  their  network  interdiction 
model  to  evaluate  nuclear  strikes  targeted  against  follow-on 
forces.  The  model  is  a  computerized  simulation  that  depicts 
real  time  movement  of  two  Soviet  divisions  moving  on  two 
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routes  from  assembly  areas  into  their  objective  areas.  The 
personnel  working  on  the  model  are  still  adding  refinements 
to  the  model  that  will  capture  all  the  effects  of  strikes 
of  this  nature.  This  work  is  ongoing  at  Sandia  National 
Laboratory,  Livermore,  California  (Ref  14). 

The  Fort  Leavenworth  group,  The  Combined  Arms  Combat 
Development  Activity  (TCADA)  (Ref  42),  is  examining  how 
revived  offensive  capability  in  the  form  of  maneuver  is 
affecting  outcomes  in  the  close-in  battle.  The  ideas  being 
examined  are  similar  to  the  ones  put  forth  by  Wagner  in  his 
articles  on  the  active  defense.  Results  of  these  studies 
stress  the  creation  of  lighter  and  more  mobile  units  with 
the  capability  to  fight  in  a  very  dynamic  battlefield 
environment . 

The  level  of  interest  in  this  second  echelon  inter¬ 
diction  problem  extends  to  other  agencies,  student  thesis 
efforts,  and  numerous  interesting  articles  in  professional 
journals.  The  DCUBE  Model  (Destroy,  Disrupt,  and  Delay), 
an  Air  Force  analytical  model  developed  for  it  by  A.  T. 
Kearney,  Inc.,  evaluates  mobility  disruption  by  air  strikes 
by  both  Blue  and  Red  air  forces  against  the  opposing  side’s 
ground  forces.  The  model  uses  an  Arrival  Rate  (AR)  sub¬ 
model  to  determine  force  arrival  rates  at  the  FLOT.  These 
results  are  then  fed  into  the  Ground  Battle  (GB)  sub-model, 
which  uses  a  modified  Lanchester  Square  Law  in  its  simulation 


The  model's  measure  of  effectiveness  (MOE)  is 
battle  outcome.  Parametric  studies  were  conducted  to  deter¬ 
mine  the  model’s  sensitivity  to  weather  effects  and  force 
variations.  Plots  against  time  were  then  made  that  show 
the  arrival  of  reinforcements  and  combat  forces  remaining 
on  the  battlefield  for  these  parametric  variations  (Ref  59). 
These  model  outputs  verify  the  same  results  the  Fort  Sill 
study  obtained  in  Figure  1.8.  "A  Simulation  of  Second 
Echelon  Air  Interdiction"  is  an  AFIT  student's  thesis  attempt 
at  correcting  weakness  in  the  DCUBE  model  to  develc  •  his 
own  air  interdiction  model.  Bennett  identified  weaknesses, 
such  as  constant  number  of  attacking  aircraft,  no  incorpora¬ 
tion  of  variance,  and  constant  kill  rates  for  airplanes  and 
trucks  (Ref  59),  that  he  thought  should  be  corrected  by  his 
formulation  in  an  attempt  to  get  a  better  representation  of 
interdiction.  His  thesis  compares  his  model's  arrival  rate 
output  against  results  of  the  DCUBE  model  to  arrive  at 
similar  conclusions.  Identified  weaknesses  in  Bennett's 
effort  include  the  following:  the  second  echelon  forces 
were  modeled  as  trucks  only  and  many  replications  are  needed 
to  get  results  (Ref  5). 

Dees'  thesis  (Ref  17)  gives  a  good  background  dis¬ 
cussion  on  second  echelon  interdiction,  how  similar  problems 
are  being  examined  in  the  Army  community,  and  on  his  recom¬ 
mendation  to  apply  Queueing-Graphical  Evaluation  and  Review- 
Technique  (Q-GERT)  in  the  development  of  his  proposed 


Model  Q-STAR  (Queueing  supplement  to  STAR  (Simulation  of 


Tactical  Alternative  Responses)  combat  model).  He  stresses 
the  utility  of  how  O-GERT  can  be  readily  used  to  model  sec¬ 
ond  echelon  forces  in  Q-STAR.  This  thesis  sets  a  foundation 
for  later  development  of  Q-STAR  (Ref  17). 

Interdiction  Planning  is  an  unpublished  paper  ob¬ 
tained  from  Major  Starner,  formerly  of  CSWS  at  Fort  Sill. 
Ihis  source  discusses  a  general  approach  to  interdiction 
planning  and  then  applies  this  methodology  to  a  particular 
sector  in  Germany.  This  article  reinforces  General  Starry's 
message  of  practicing  interdiction  planning  now  with  current 
resources  so  that  expertise  is  developed  at  all  levels. 
Targeting  is  examined  for  both  division  and  corps  levels 
(Ref  54) . 

An  examination  of  the  remainder  of  the  bibliograph¬ 
ical  references  shows  other  articles,  books,  and  studies  on 
work  done  by  both  Army  and  Air  Force  personnel.  Included 
are  also  works  by  other  governmental  agencies  or  contractors 
The  level  of  interest  displayed  is  indicative  of  the  impor¬ 
tance  attached  to  the  particular  problem  of  second  echelon 
interdiction . 

Methodology 

A  systems  approach  was  the  methodology  applied  in 
this  research  effort.  Since  the  subsystem  described  pre¬ 
viously  is  part  of  a  large  and  more  complex  system,  a 


systems  approach  can  best  be  used  for  this  policy  analysis 
because  it  aids  in  understanding  the  system  more  fully. 

The  Systems  Approach.  By  using  the  process  described 
by  Schoderbek,  Schoderbek,  and  Kef alas  in  Management  Systems 
Conceptual  Considerations,  a  more  complete  understanding  can 
be  gained  about  that  system  so  that  the  perceived  problem  or 
required  policy  to  be  studied  can  be  isolated  into  its  com¬ 
ponent  parts  for  analysis:  "the  input(s),  the  process(es), 
the  output (s),  and  the  feedback  control"  (Ref  50:14).  A 
mathematical  model  can  then  be  formulated  to  represent  this 
isolated  part  of  the  system  so  that  an  analysis  can  be  made 
between  feedback  structures  and  system  components  as  they 
interact  over  time.  Once  the  model  is  built,  studies  can 
be  conducted  on  the  model  to  obtain  insights  on  how  changes 
affect  the  system  or  to  analyze  new  problems.  The  value  of 
this  effort  does  not  lie  solely  in  information  obtained,  but 
also  in  gaining  insights  about  the  system  all  during  the 
model  development. 

Thus,  systems  analysis  can  provide  a  valuable  tool 
in  correctly  specifying,  delimiting,  and  understanding  the 
problem  to  be  studied  prior  to  time  being  wasted  on  answering 
the  wrong  question. 

Applied  Methodology.  How  system  analysis  was  used 
in  the  research  phase  follows.  Since  second  echelon  inter¬ 
diction  is  a  subsystem  of  the  main  system  of  war,  a  narrowing 
of  the  overall,  general  system  was  necessary  to  isolate  key 


relationships  essential  to  understanding  how  the  subsystem 
operates  and  interacts.  A  cyclic  process  was  used  through¬ 
out  the  model  development  until  the  model  was  finally 
accepted  as  an  accurate  system  representation  of  the  inter¬ 
diction  process.  This  iterative  process  was  a  constant,  on¬ 
going  thing  to  insure  system  accuracy.  Figure  1.9  depicts 
this  process  as  it  relates  to  the  simulation  process. 

Initial  system  conceptualization  and  ideas  were 
obtained  on  visits  to  the  CSWS  office  at  Fort  Sill,  Okla¬ 
homa,  and  Sandia  National  Laboratory  at  Livermore,  California. 
Further  readings  and  research  aided  in  further  understanding 
the  system  process.  From  those  efforts,  system  causal 
relationships  were  developed. 

A  thorough  map  analysis  of  a  road  network  in  East 
Germany  identified  choke  points  on  the  routes  hypothesized 
for  movement  of  a  motorized  rifle  division  (MRD) .  From  the 
identified  choke  points  emerged  an  idea  of  where  to  cut  the 
roads  with  general  purpose  bombs  so  that  convoy  passage 
would  be  blocked.  The  MRD  was  selected  as  moving  on  three 
separate  routes  to  minimize  the  convoy  length  and  to  repre¬ 
sent  WP  doctrine  of  maintaining  unit  dispersal.  Unit 
characteristics,  such  as  number  and  type  of  units,  types  of 
vehicles,  convoy  lengths,  were  obtained  from  Opposing  Forces 
Europe  (Ref  31).  The  terrain  restricted  points  on  the  route 
imply  that  the  terrain  restricts  the  movement  of  the  units 
to  the  roadway  only  because  high  ground  to  either  side  of 
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the  road  Goes  not  allow  bypass.  Also,  selecting  these 
restricted  points  included  insuring  that  no  other  alternate 
roads  in  the  immediate  area  would  allow  easy  routine  around 
the  obstacle.  Targeting  an  area  other  than  these  points 
would  not  impose  appreciable  delay  or  create  unit  bunching 
because  units  could  bypass  the  obstacle  by  going  around  it 
either  by  use  of  an  alternate  route  that  easily  puts  them 
back  on  their  original  route  or  by  simply  going  off  the 
road  and  back  on  the  road  again.  Once  the  convoys  are 
blocked,  a  relook  and  retargeting  capability  would  allow  for 
retargeting  of  the  area  with  one  or  more  missiles  loaded 
with  wide  area  antiarmor  munitions  (WAAMf.  ) .  This  targeting 
process  should  tell  a  commander  where  to  make  the  necessary- 
cuts  so  that  he  can  destroy,  delay,  or  disrupt  the  flow  of 
incoming  enemy  forces  and  thereby  create  a  tactical  situa¬ 
tion  that  allows  for  destruction  ol  the  first  echelon  while 
delaying  the  arrival  of  the  second  echelon.  The  above 
implementation  is  described  in  subsequent  chapters. 

Order  of  Presentation 

The  next  six  chapters  present  the  implementation  of 
the  research  methodology,  policy  evaluation,  and  recommenda¬ 
tions  for  further  research. 

Chapter  II  is  the  model  Description  and  Development. 
Within  it,  key  points  discussed  include  how  parameters  were 
developed,  what  assumptions  were  used,  how  the  interdiction 
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was  formulated,  and  how  all  parts  of  the  model  were  inte¬ 
grated. 

The  Simulation  Model  is  next  in  Chapter  III.  Within 
this  chapter,  how  the  model  was  computerzied  is  described. 

Validation  and  Verification,  Chapter  IV,  disucsses 
the  process  of  verifying  the  functioning  of  the  model. 

Data  Collection,  Chapter  V,  discusses  experimental 
design  and  sample  size  determination. 

Chapter  VI  presents  the  data  analysis  for  the  experi¬ 
mental  design.  Finally,  Chapter  VII  presents  the  Conclusions 
and  Recommendations,  and  Recommendations  for  Follow-on  Study. 

Summary 

Presented  in  Chapter  I  has  been  the  problem,  the 
research  question,  a  background  on  research,  and  the  applied 
methodology  for  this  research  project.  Chapter  II  follows 
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System  Structure 


The  model  STOPEM  was  developed  to  analyze  the  problem 
stated  in  Chapter  I.  The  purpose  of  this  chapter  is  to  pre¬ 
sent  the  methodology  involved  in  developing  the  structure 
for  STOPEM.  The  main  components  of  this  structure  are: 
convoy  movement,  sortie  attack,  missile  attack,  damage  assess 
ment,  and  delay  assessment.  STOPEM' s  initial  development 
consisted  of  modeling  the  physical  convoy  movement  of  the 
44  motorized  rifle  division  (MRD)  units  through  three  sepa¬ 
rate  routes.  Once  the  physical  movement  development  was 
completed,  the  actual  interdiction  and  damage  assessment  of 
the  convoys  were  added.  This  particular  method  aided  in 
model  verification.  The  following  will  trace  the  develop¬ 
ment  of  these  ideas. 

Convoy  Movement 

Both  Opposing  Forces  Europe  (Ref  31)  and  Soviet  Army 
Operations  (Ref  32)  were  consulted  to  obtain  the  necessary 
information  to  build  this  initial  portion  of  the  model. 

Soviet  Army  contains  the  basic  Soviet  Doctrine  on  convoy 
movement,  convoy  intervals,  rates,  and  use  of  multiple 
routes  to  maintain  dispersion.  Appendix  A  of  Opposing  Forces 
lists  descriptions  of  units  from  company  to  Warsaw  Pact 
FRONT,  their  equipment,  and  personnel  assigned. 


From  Opposing  Forces  Europe  (Ref  31),  a  MRD  was 
broken  out  into  44  different  units  of  approximately  battal¬ 
ion  size.  The  exceptions  to  this  breakdown  are  division 
headquarters  (modeled  as  two  ent ities--main  and  alternate, 
both  of  equal  size),  the  regimental  headquarters  (five), 
three  anti-aircraft  batteries,  a  target  acquisition  battery, 
and  the  servi ces-rear  security.  Table  2.1  shows  the  break¬ 
down  into  the  various  elements,  their  route  designation,  and 
the  average  length  of  the  convoy.  The  three  anti-aircraft 
batteries  are  detached  from  their  regiment  unit  to  provide 
air  defense  within  the  three  columns.  The  only  MRD  unit  not 
considered  in  this  formulation  is  the  division  reconnaissance 
battalion.  This  unit  is  omitted  because  it  would  be  operat¬ 
ing  about  50  kilometers  in  front  of  the  MRD  (Ref  32:9-4). 
Thus,  it  would  not  be  cost-effective  to  delay  this  unit 
since  it  could  radio  back  to  the  main  force  about  blockages 
ahead  and  that  would  allow  the  MRD  an  early  option  to  select 
an  alternate  route. 

Sampling  from  a  triangular  distribution  is  conducted 
for  convoy  rate,  convoy  length,  and  convoy  interval  because 
variables  are  considered  to  be  stochastic  in  nature.  Since 
data  for  these  parameters  in  Soviet  Anny  Operations  (Ref  32) 
is  given  in  minimum  and  maximum  values,  a  triangular  distri¬ 
bution  is  selected  as  the  sampling  distribution  because  the 
description  of  the  process  and  the  absence  of  real  data  fit 
the  circumstances  when  this  distribution  could  be  used.  In 
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TABLE  2.1 


Motorized  Rifle  Division  Units 


Unit 

Number 

Description 

Route 

Average  Le 
(KM) 

1 

1st  MR  Regiment  Advanced 

Guard 

Blue 

3.6 

2 

1st  Regiment  HQs . 

Blue 

1 . 13 

3 

Tank  Bn  (-) 

Blue 

1 . 35 

4 

MR  Bn 

Blue 

1.65 

5 

MR  Bn  (-) 

Blue 

1.5 

6 

Services  -  Rear  Security 

Blue 

3.75 

7 

Tank  Bn 

Blue 

2.1 

8 

Division  HQs  (alternate) 

Blue 

1.31 

9 

Target  Acquisition  Battery 

Blue 

0.6 

10 

Artillery  Bn  (152  mm) 

Blue 

2.4 

11 

Artillery  Bn  (122  mm) 

Blue 

2.4 

12 

AAA  Battery 

Blue 

0 . 75 

13 

Anti  -  Tank  Bn 

Blue 

2 . 85 

14 

FROG  Bn 

Blue 

1.95 

15 

Transportation  Bn 

Blue 

7.65 

16 

2nd  MR  Regiment  Advanced 

Guard 

Green 

3.6 

17 

2nd  Regiment  HQs 

Green 

1.13 

18 

Tank  Bn  (-) 

Green 

1.35 

19 

MR  Bn 

Green 

1.65 

20 

MR  Bn  (-) 

Green 

1.2 

21 

Services  -  Rear  Security 

Green 

3.75 

22 

Artillery  Bn  (122  mm) 

Green 

2.0 

23 

Artillery  Regiment  HQs 

Green 

2.25 

24 

AAA  Regiment  (-) 

Green 

2.25 

25 

Engineer  Bn 

Green 

3.75 

26 

Signal  Bn 

Green 

2.25 

27 

MLRS  Bn 

Green 

2.7 

28 

Chemical  Bn 

Green 

1.65 

29 

Maintenance  Bn 

Green 

2.55 

30 

AAA  Battery 

Green 

0. 75 

31 

Medical  Bn 

Green 

1.35 

32 

Tank  Regiment  Advanced  Guard 

Red 

2.25 

33 

Tank  Regiment  HQs 

Red 

0.75 

34 

Tank  Bn 

Red 

1.5 

35 

Tank  Bn  (-) 

Red 

1.2 

Note:  (-)  =  Unit  minus 

some 

forces 
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TABLE  2.1  (Continued) 


Un  i  t 
Number 


Descript i on 


Average  Length 
Route  (KM; 


other  words,  a  triangular  distribution  "is  used  when  a  most 
liekly  value  can  be  ascertained  along  with  minimum  and  maxi¬ 
mum  values,  and  a  piece-wise  linear  density  function  seems 
appropriate"  (Ref  49:30).  The  mode  for  the  distribution  is 
calculated  as  the  average  of  the  minimum  and  maximum  values 
given.  Once  these  three  parameters  are  available,  sampling 
can  be  conducted.  A  flexibility  of  this  distribution  is  its 
ability  to  vary  the  mode  to  meet  changing  conditions.  For 
example,  the  mode  for  the  rate  would  be  expected  to  decrease 
over  time  because  of  the  effects  of  degradation  of  road 
surfaces  and  interdiction. 

Next,  the  parameters  for  convoy  interval,  rate,  and 
length  are  obtained  as  follows.  The  minimum  and  maximum 
values  for  the  above  parameters  are  obtained  from  Soviet 
Army  Operations  (Ref  32:3-20,3-21).  The  interval  between 
all  units  is  taken  to  be  between  three  and  five  kilometers, 
with  a  mode  of  four.  The  exception  is  the  advanced  guard 
element  on  a  route.  The  interval  between  it  and  the  main 
body  of  the  convoy  is  between  20  and  30  kilometers  ,  with  a 
mode  of  25.  A  mixed  convoy  rate  for  daytime  conditions  is 
selected  because  a  day  road  march  is  assumed.  A  day  road 
march  is  assumed  due  to  constant  convoy  movement  forward. 

The  values  for  the  minimum,  mode,  and  maximum  are  20,  25,  and 
30  kilometers  per  hour,  respectively.  Minimum  and  maximum 
convoy  lengths  are  computed  by  multiplying  the  number  of 
vehicles  in  the  convoy  (discussed  in  subsequent  paragraph) 
by  the  minimum  and  maximum  interval  between  vehicles 


(25  and  50  meters,  respectively)  plus  the  minimum  and  maxi¬ 
mum  intervals  between  companies  within  a  battalion  (25  and 
50  meters,  respectively).  This  method  of  computation 
assumes  a  measurement  interval  from  mid-vehicle  to  mid¬ 
vehicle.  Once  these  minimum  and  maximum  unit  lengths  are 
computed,  an  average  value  (the  mode)  is  computed.  This 
value  appears  in  Table  2.1  as  the  average  length. 

To  compute  unit  minimum  and  maximum  lengths,  the 
number  of  vehicles  in  a  convoy  has  to  be  determined. 

Opposing  Forces  Europe  has  what  appears  to  be  an  incomplete 
listing  of  the  number  of  vehicles  under  the  MRD  table 
(Ref  31:A-12).  For  example,  the  engineer  battalion  has 
385  personnel,  but  only  10  wheeled  vehicles  to  carry  these 
personnel  and  equipment.  This  is  resolved  by  going  to 
another  page  of  the  same  appendix  to  find  the  same  unit  with 
a  complete  listing  of  vehicles.  In  other  cases,  units  listed 
do  not  include  any  vehicles  at  all.  In  these  cases,  a 
similar  organization  found  under  a  separate  listing  is  used 
as  a  comparable  substitute.  For  example,  the  data  for  the 
medical  battalion  was  found  under  divisional  services  for 
the  MRD.  In  Figure  2  of  Soviet  Army  Operations,  a  regimental 
headquarters  is  shown  being  augmented  with  another  unit  such 
as  a  regimental  artillery  group.  Since  this  reoresents  a 
particular  attachment  for  combat  and  not  a  general  case  as 
this  model  portrays,  regimental  headquarters  in  the  model 
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are  portrayed  without  special  attachments.  Of  course, 
having  additional  attachments  implies  that  the  unit  length 
parameters  would  increase.  What  these  examples  attempt  to 
portray  is  that  exact  figures  for  all  units  are  not  available 
so  cross  checking  made  above  is  used  to  calculate  the  number 
of  vehicles.  The  sources  examined  represent  the  most  com¬ 
plete  unclassified  listing  available. 

All  track  vehicles  (tanks  (all  types),  infantry 
fighting  vehicles  (BMP),  self-propelled  artillery  pieces, 
and  engineer  vehicles)  and  all  wheeled  vehicles  (towed 
artillery  pieces,  trucks,  and  reconnaissance  vehicles)  are 
grouped  into  either  track  or  wheeled  vehicle  categories. 

Once  total  number  of  vehicles  per  unit  is  determined,  unit 
length  is  computed.  For  example,  unit  one  has  57  track 
vehicles  and  32  wheeled  vehicles.  Units  1,  16,  and  39  are 
all  advanced  guard  elements  for  three  different  motorized 
rifle  regiments.  This  element  is  a  motorized  rifle  battalion 
reinforced  with  an  armor  company,  a  reconnaissance  company, 
mortars,  an  engineer  detachment,  an  artillery  battery  and 
various  combat  support  vehicles.  To  compute  unit  minimum, 
maximum,  and  average  lengths  for  this  unit,  the  following 
is  done : 

min  length  =  (number  of  vehicles) (minimum  interval)  + 
(number  of  companies) (minimum  interval) 

=  ( 89 ) (25 )  +  (6 ) (25 ) 

=  2375  meters 

=  2.375  a  2.4  kilometers 
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max  length  =  (number  of  vehicles) (maximum  interval) 
(number  of  companies ) (maximum  interval) 


=  ( 89 ) ( 50  )  +  ( 6 ) ( 50 ) 

=  4750  meters 
=  4.75  *  4.8  kilometers 


un ^ t fo  _  minimum  length  +  maximum  length 

length  2 

2.4  +  4.8 
2 

=  3.6  kilometers 

A  similar  calculation  is  done  for  all  other  units.  Once  all 
intervals  are  calculated,  convoy  number  plus  values  for 
number  of  track  and  wheeled  vehicles,  unit  lengths,  and 
unit  rate  characterize  each  convoy  entity. 

Convoy  order  on  each  route  is  established  with  an 
advanced  guard  element  at  the  head,  with  the  remainder  of 
each  regiment  and  other  units  behind  it.  The  division  main 
and  alternate  headquarters  are  placed  on  separate  routes 
for  redundancy.  Combat  service  support  elements  are  placed 
to  the  rear  of  columns,  and  the  main  combat  elements  and 
their  controlling  headquarters  are  placed  toward  the  head 
of  each  column. 

Figure  2.1  is  a  static  depiction  of  how  convoy 
movement  is  modeled  and  Figure  2.2  is  a  simplified  road 
network  depiction.  As  mentioned  earlier,  this  model  assumes 
a  hypothesized  three  route  move  for  a  MRD.  Figure  2.2  shows 
the  representation  of  the  main  route  for  each  convoy.  The 


numbers  between  nodes  (circles)  represent  the  distance  in 
kilometers  between  these  points.  Convoys  are  assumed  to 
have  been  moving  for  50  kilometers  prior  to  entering  the 
road  network  of  interest.  These  values  were  obtained  on  a 
visit  to  Sandia  National  Laboratory  and  at  random  verified 
using  manual  plotting  on  a  military  map  of  the  area  in 
Germany.  This  is  a  simplified  network  because  the  extremely 
complex  road  network  atailable  (other  roads  and  forest 
trails)  is  not  shown.  Since  this  model  assumes  a  hypothe¬ 
sized  travel  route  based  on  Warsaw  Pact  doctrine,  these 
routes  are  selected  to  approximate  where  convoys  might 
actually  travel.  Based  on  this  route  of  travel,  terrain- 
restricted  points  are  selected  for  interdiction.  In  most 
cases,  a  node  is  a  terrain  restricted  point. 

The  convoys  are  modeled  as  having  a  head  and  a  tail 
entity.  The  vehicles  in  between  represent  the  time  it  takes 
for  a  convoy  with  all  its  elements  to  pass  from  one  point 
to  another  point.  This  is  referred  to  as  passage  time. 

Since  the  length  of  the  convoy  is  stochastic,  the  vehicles' 
distance  in  between  the  head  and  tail  of  the  convoy  varies 
by  a  triangular  distribution.  The  head  of  the  lead  unit  in 
each  convoy  starts  at  "time  now"  (TNOW) ,  which  represents 
zero  time.  The  tail  of  the  convoy  is  computed  to  start 
behind  the  head  of  the  convoy  at  the  time  TNOW  plus  the  unit 
length  divided  by  the  rate.  Dividing  the  unit  length  by  the 
rate,  also  drawn  from  a  triangular  distribution,  converts 


the  distance  to  time  so  TNOW  can  be  added  to  it. 


All  t i me  s 


are  in  hours.  Since  the  interval  between  convoys  is  also 
stochastic,  a  sample  taken  from  a  triangular  distribution 
divided  by  the  rate  of  unit  one  determines  how  far  back  the 
second  unit  head  is  located.  Adding  this  time  (interval 
between  units  divided  by  rate)  to  the  time  for  the  tail 
element  of  the  first  convoy  gives  the  starting  time  for  the 
head  of  the  second  convoy.  As  with  convoy  one,  the  loca¬ 
tion  of  the  tail  element  for  convoy  two  is  computed  after 
sampling  frc m  a  triangular  distribution  with  parameters  for 
unit  two's  lengths.  This  procedure  continues  until  all 
convoys  are  on  the  road  with  a  head  and  tail  entity.  The 
march  rates  for  all  convoys  behind  the  first  unit  on  each 
route  are  adjusted  to  insure  that  the  units  finish  in  the 
same  order  in  which  they  start .  Rates  are  adjusted  by  decreasing 
or  increasing  the  rate,  depending  on  whether  it  is  greater 
or  less  than  the  lead  rate.  In  actuality,  the  rate  of  the 
lead  unit  determines  the  march  rate  for  units  following  it 
because  failure  to  maintain  a  rate  less  than  or  equal  to 
the  lead  rate  would  create  congestion  or  an  accordion  effect 
on  the  vehicles. 

Sortie  Attack 

The  movement  of  the  convoys  in  the  preceding  section 
represents  the  steady  state  condition  of  the  system  under 
study.  The  time  to  travel  the  portion  of  the  network  being 
examined  represents  the  time  that  the  MRD  would  take  without 
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any  impediments  such  as  interdiction.  This  section  will 
outline  how  air  interdiction  might  affect  the  steady  state 
time.  Examining  the  system  with  interdiction  will  give 
insights  into  how  additional  time  delays  will  affect  arrival 
times  of  convoys  at  the  battle  area.  The  difference  between 
the  steady  state  time  and  the  time  with  interdiction  will 
represent  the  delay  time.  Also,  the  destruction  of  vehicles 
computed  as  a  result  of  interdiction  represents  the  decrease 
in  combat  power  of  the  MRD. 

A  flight  of  two  aircraft,  each  loaded  with  12 
500-pound  MK-82  general  purpose  bombs,  represents  the  strike 
mission.  The  assumption  is  made  that  the  aircraft  are 
F4D/E,  with  a  6000-foot  slant  range,  in  poor  combat  condi¬ 
tions,  using  a  dive  toss  bombi  g  run.  Each  aircraft  is  also 
assumed  to  make  one  pass  and  drop  all  its  bombs  in  pairs. 
Figure  2.3  depicts  the  assumed  pattern  of  impact.  Consulting 
JE.MM ,  Weapons  Characteristics  (Ref  39:1-184G),  an  unclassified 


estimate  was  obtained  on  crater  diameter  for  a  MK-82  dropped 
on  12  inches  of  concrete  with  medium  soil  underlayer.  JEMM 
says  to  treat  roads  as  runways.  The  crater  estimate  for 
this  12  inches  of  concrete  was  about  35  feet.  Even  though 
the  road  being  modeled  is  not  12  inches  of  concrete,  the 
crater  diameter  of  the  bomb  is  assumed  to  be  30  feet.  In 
reality,  the  crater  diameters  would  probably  be  greater  due 
to  the  road  material  not  being  12  inches  of  concrete.  The 
int ravolometer  setting  is  assumed  to  be  set  so  that  the 


timing  of  release  between  bombs  would  give  the  pattern  in 
Figure  2.3.  A  25-mil  circular  error  probable  (CEP)  is  also 
assumed  to  arrive  at  a  150-foot  CEP  in  the  ground  plane. 

The  150-foot  CEP  was  computed  by  using  the  relationship 
that  1  mil  at  1000  feet  equals  1  foot.  So,  multiplying  the 
slant  range  by  the  25-mil  CEP  gives  150  feet. 

Since  the  intention  of  this  model  is  to  capture 
interdiction  effects,  and  not  aircraft  interdiction,  the 
interdiction  of  the  roads  is  modeled  by  determining  the 
probability  of  cutting  an  assumed  50-foot  road  width.  To 
gain  an  understanding  of  the  geometry  involved,  a  model  was 
built  to  examine  the  problem.  This  appears  in  Appendix  B. 
By  varying  different  parameters,  such  as  CEP,  range  error 
probable  (REP),  deflection  error  probable  (DEP),  attack 
angle,  target  width,  and  target  length,  an  estimate  was 
obtained  of  the  probability  of  cutting  the  assumed  road  for 
1000  bomb  drops.  Next,  a  group  of  fighter  pilots  was  con¬ 
sulted  about  the  reasonableness  of  the  estimate  obtained. 
Conditions  as  outlined  above  were  also  presented  to  them. 
Based  on  the  geometry  model  parametric  studies  and  these 
conversations,  the  probability  arrived  at  was  0.35.  Thus, 
based  on  the  given  conditions,  roads  are  assumed  to  be  cut 
35  percent  of  the  time  by  a  two-aircraft  sortie. 

Another  assumption  made  was  that  aircraft  penetrate 
the  FEBA  without  loss.  Also,  enemy  engagement  of  the 
attacking  aircraft  is  ignored,  but  partially  considered 


in  the  CEP  calculation  by  assuming  poor  combat  conditions. 

A  higher  CEP  allows  for  target  area  unfamiliarity  and  the 
threat  environment. 

While  a  specific  scenario  is  fixed  for  this  study, 
in  reality  the  aircraft  interdiction  process  is  a  complex  one 
that  is  full  of  variability  and  uncertainty.  For  instance, 
each  aircraft  has  its  own  system  weaknesses  and  strengths, 
pilots  have  different  levels  of  proficiency,  delivery  condi¬ 
tions  are  dependent  on  such  factors  as  air  speed,  delivery 
angle  and  target  type,  and  the  threat  environment  is  such 
that  the  probability  of  success  of  aircraft  penetration, 
target  attack,  and  egress  is  decreased. 

The  total  number  of  daily  sorties  available  for  the 
corps  sector  is  assumed  to  be  30.  Considered  was  that  a 
percentage  of  aircraft  would  be  nonoperational  and  that 
about  10  percent  of  the  total  aircraft  available  would  be 
committed  to  other  special  missions.  From  the  remaining 
number  of  aircraft,  an  allocation  has  to  be  made  of  aircraft 
to  interdict  this  three-route  move  of  the  MRD.  A  further 
complication  involves  knowing  that  the  MRD  under  study  is 
probably  one  of  several  forces  in  the  corps  sector  that 
require  interdiction  such  as  Figure  1.2  depicts.  Thus,  the 
allocation  of  sortie  resources  within  the  corps  sector  is 
among  competing  activities.  The  actual  basis  for  allocation 
will  involve  the  design  of  the  experiment,  to  be  discussed 
in  Chapter  IV. 


Damage  Due'  to  Sortie  Attack 

The  convoys  are  assumed  to  be  targeted  at  the  head 
element.  To  simplify  calculations  of  vehicles  damaged  due  to 
the  air  strike,  any  vehicles  within  the  radius  of  the  bomb 
impacts  are  considered  destroyed.  Because  the  majority  of 
targets  bombed  are  hard  targets  (tanks),  the  assumption  made 
is  that  bombing  would  kill  at  most  one  vehicle  per  pass. 

So,  a  total  of  two  vehicles  are  assumed  killed.  Once  the 
destroyed  vehicles  are  subtracted,  the  convoy  length  is 
adjusted  to  reflect  the  loss  in  vehicles. 

Delay  Due  to  Sortie  Attack 

An  Algorithm  for  Determining  Delays  Imposed  on  Ground 
Forces  I  ,e  To  Interdiction  Air  Strikes  Revisited,  Technical 
Paper  5-79  (Ref  3),  was  consulted  to  understand  more  about 
air  interdictions  and  the  effects  on  ground  forces.  Other 
papers  examined  include  Minutes  of  an  Exploratory  Meeting  on 
Interdiction  Study,  AC  243,  NATO  Panel  VII  (Ref  46,  NATO 
RESTRICTED)  and  An  Algorithm  for  Determining  Delays  Imposed 
on  Ground  Forces  Due  to  Interdiction  Air  Strikes,  Technical 
Paper  3-79  (Ref  25,  NATO  RESTRICTED).  These  two  references 
were  used  for  background  information. 

TP  5-79  explains  the  "four  mutually  exclusive, 
exhaustive  delay  events,  with  associated  delay  times  defined 
as  follows:" 

Heads  Down  Time  (HD)--the  initial  reaction  of  the 

Red  Force  to  Blue  air  attack. 


47 


Damage  Assessment  Time  (DA)--the  time  required 
for  a  commander  to  receive  reports  from  his 
subordinates,  assess  the  situation,  and  then 
report  to  his  higher  headquarters. 

Damage  Control  Time  (DC)--the  time  required 
to  treat  personnel  casualties  (first  aid)  and 
recover  damaged  and/or  slightly  damaged 
vehicles . 

Impaired  Movement  Time  (IM) — the  time  required 
to  bypass  craters  and  remove  burning  vehicles 
and  other  obstacles  from  the  road  (Ref  3:1). 

This  study  does  not  calculate  the  individual  times  for  the 
above  four  events.  Instead,  based  on  ideas  obtained  from  the 
NATO  Panel  study,  an  assumed  distribution  for  delay  time  is 
used.  If  the  convoy  route  is  cut  and  block,  sampling  is  made 
from  a  triangular  delay  distribution  with  minimum,  mode,  and 
maximum  values  of  30,  34,  and  38,  respectively.  This  time 
is  in  minutes.  If  the  road  is  not  cut,  then  a  delay  time  is 
still  imposed  to  capture  the  time  associated  with  getting 
the  convoy  reoriented,  reporting  to  higher  headquarters,  and 
negotiating  the  road  damaged  by  the  strike.  Sampling  is  made 
from  another  triangular  distribution  containing  parameters 
for  minimum,  mode  and  maximum  values  of  10,  12.5,  and  15, 
respectively.  Triangular  distributions  are  selected  for  the 
same  reason  given  in  development  of  convoy  movement  para¬ 
meters  . 

Since  the  lead  element  of  each  convoy  is  a  reconnais¬ 
sance  element  operating  in  front  of  the  convoy  at  greater 
distance  than  the  interval  between  the  other  units,  inter¬ 
diction  of  this  lead  element  allows  it  to  radio  back  the 


resulting  attack,  whether  successful  or  unsuccessful.  The 
rest  of  the  convoy  is  then  assumed  to  take  an  alternate 
route  to  bypass  the  obstacle  or  damaged  road  on  the  original 
route.  The  lead  element  is  assumed  to  be  delayed  by  the 
amount  required  to  repair  and  negotiate  the  obstacle  and 
continue  on  the  original  route.  The  idea  is  to  continue 
moving  the  rest  of  the  force  and  minimize  delay  or  bunching 
of  vehicles  as  an  obstacle  is  approached  and  negotiated. 

Tue  lead  element  is  not  allowed  to  turn  around  and  attempt 
to  resume  its  lead  place  in  the  pack.  Convoys  continue  in 
the  route  specified  without  attempts  to  remerge  them  because 
to  merge  remnants  with  the  original  convoy  requires  waiting 
time  for  some  part  of  the  convoy.  This  is  considered 
unrealistic  given  that  Warsaw  Pact  doctrine  requires  con¬ 
stant  movement  forward. 

If  the  lead  element  is  not  interdicted,  but  allowed 
to  pass  on,  then  the  next  element  of  the  convoy  will  be 
bombed.  By  allowing  the  reconnaissance  element  to  go  on 
uninterdicted,  the  main  column  presents  a  bigger  target  of 
opportunity  due  to  shorter  intervals  between  units.  The 
reaction  time  to  radio  back  the  bombing  and  the  issuance  of 
orders  to  change  routes  would  cause  the  following  units  to 
quickly  close  the  interval.  Assuming,  an  average  interval 
of  4  kilometers  and  an  average  rate  of  25  kilometers  per 
hour,  a  unit  following  the  interdicted  one  would  close  the 
distance  between  them  in  about  10  minutes.  Unless  immediate 


This  model  assumes 
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action  is  taken,  units  begin  to  pile  up. 
that  the  second  unit  will  not  be  able  to  find  a  bypass  and 
waits  behind  the  lead  unit  until  the  obstacle  is  repaired. 
The  third  unit  and  others  behind  it  are  assumed  to  have 
received  the  order  to  take  an  alternate  bypass.  As  before, 
the  blocked  units,  once  the  obstacle  is  negotiated,  do  not 
rejoin  the  rest  of  the  convoy  at  a  point  down  the  original 


route , 


A  few  observat ions  are  in  order  to  discuss  the 


implications  of  the  above  convoy  actions.  The  unit’s  elec¬ 
tronic  signature  surfaces  when  communication  silence  is 
broken  to  report  the  attack  to  higher  headquarters.  This 
action  may  facilitate  target  acquisition  or  detection. 

Also,  interdicting  the  reconnaissance  element  in  the  column 
and  successfully  cutting  the  road  will  leave  the  rest  of 
the  convoy  without  a  quick  obstacle  repair  capability  to  its 
front  because  the  lead  element  will  have  to  fix  or  negotiate 
the  obstacle.  Taking  away  this  capability  by  delaying  and 
separating  it  from  the  main  column  reduces  quick  obstacle 
re  air.  It  also  reduces  MRD  advance  warning  time  about 
future  activities.  At  any  future  obstacles,  the  convoy  will 
have  to  depend  on  what  is  at  hand  to  repair,  that  is,  man¬ 
power,  or  wait  for  some  type  of  engineer  equipment  to  make 
its  way  up  to  affect  repair.  On  the  other  hand,  interdicting 
this  element  will  give  the  main  body  of  the  convoy  more  time 
to  react  and  select  another  route.  This  would  allow  the  main 
body  to  move  on  without  slowing  it  down. 
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Interdicting  he  lead  element  of  the  main  body 
appears  to  offer  more  in  the  way  of  destruction  and  disrup¬ 
tion.  Not  only  is  there  less  reaction  time  available  to 
following  units  to  make  a  quick  route  change,  but  there  is 
more  bunching  of  units.  A  more  dense  and  lucrative  target 
is  created  for  wic  -area  anti-armor  munitions  (WAAM).  The 
possibility  then  exists  for  greater  disruption  and  destruc¬ 
tion.  By  modeling  the  process  of  moving  separated  units 
into  different  routes,  congestion  at  some  future  point  is 
eliminated.  Otherwise,  units  will  be  stopped  while  waiting 
for  another  unit  to  be  inserted.  Vulnerability  is  decreased, 
but  disruption  is  increased. 

All  the  above  assume  perfect  intelligence  and  good 

weather . 

Missile  Attack 

Once  the  road  is  interdicted,  a  determination  needs 
to  be  made  about  retargeting  the  same  area  with  a 
ground  launched  Army  missile  loaded  with  wide-area  anti¬ 
armor  munitions  (WAAM).  As  mentioned  previously,  knowing  that 
there  is  a  stationary  target  located  behind  a  road  blockage 
removes  the  uncertainty  of  target  location.  Cratering  or 
damaging  the  road  with  air  interdiction  should  create  a 
target  suitable  for  WAAM  engagement.  A  generic  missile 
capable  of  carrying  WAAM  is  assumed  to  have  an  approximate 
CEP  of  110  feet.  This  CEP  is  obtained  by  using  unclassified 
Eq  (1)  taken  from  a  letter  requesting  data  from  the  Corps 
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Support  Weapon  System  (CSWS)  office,  Fort  Sill,  Oklahoma 
(Ref  15:3)  (CONFIDENTIAL).  The  system  CEP  is  assumed  to  have 
a  circular  normal  distribution. 

CEP  System  =  /  CEPMPI2  +  CEPP2 

where  CEPMPI  =  0 . 5877(RSIGMPI  +  DSIG.MPI) 

CEPP  =  0 . 5887( RSIGP  +  DSIGP)  (1) 

Values  assumed  for  computing  system  CEP  follow: 

Range:  100  KM 

RSIGMPI :  20  M  (65  feet) 

(Range  deviation,  mean  point  of  impact) 

DSIGMPI :  20  M  (65  feet) 

(Deflection  deviation,  mean  point  of  impact) 

RSIGP:  20  x.l  (65  feet) 

(Range  deviation,  precision) 

DSIGP:  20  M  (65  feet) 

(Reflection  deviation,  precision) 

The  shape  of  the  dispersal  pattern  for  the  WAAM  is  assumed 
to  be  elliptical  with  a  minor  axis  value  of  46  meters 
(approximately  150  feet)  and  major  axis  value  of  240  meters 
(approximately  800  feet).  All  vehicles  within  this  pattern 
are  assumed  killed.  Assumptions  made  for  the  missile  include 
reliability  of  one  ( for  pre-launch  and  in-flight)  and  perfect 
communication  and  intelligence  to  allow  for  notification  of 
air  interdiction  mission.  Also,  another  assumption  made  is 
that  the  WAAM  always  impacts  with  the  major  axis  of  damage 
parallel  to  the  road.  As  with  the  air  sortie  interdiction, 


the  missile  attack  is  not  modeled  explicitly.  How  it  is 
done  follows  in  the  next  section. 

Damage  Due  to  Missile  Attack 

A  probability  density  function  is  computed  to  deter¬ 
mine  damage  due  to  WAAM  attack.  Knowing  that  the  minimum 
spacing  and  maximum  spacing  between  vehicles  are  25  and  50 
meters  (approximately  80  and  160  feet),  respectively,  a 
methodology  is  needed  to  determine  the  probability  of  killing 
vehicles  based  as  a  function  of  vehicle  spacing  and  munition 
impact.  Vehicles  are  assumed  to  be  on  the  roadway  center 
line.  An  average  of  both  extremes  of  spacing  is  determined 
to  be  120  feet  [(80  +  160)  f  2].  Vehicle  spacing  is  deter¬ 
mined  by  dividing  the  unit  length  by  the  number  of  vehicles. 
This  number  is  then  converted  to  feet.  Once  spacing  between 
vehicles  is  determined,  the  value  obtained  is  compared 
against  this  average.  If  the  value  obtained  is  less  than 
the  average,  then  the  assumed  interval  between  vehicles  is 
80  feet,  and  160  feet  otherwise. 

Using  the  minimum  interval  of  80  feet  and  an 
assumed  vehicle  length  of  20  feet,  an  average  of  8  targets 
are  hypothesized  as  being  within  a  road  cut  made  by  the 
ellipse  of  between  750  and  800  feet.  Figure  2.4  depicts 
the  elliptical  pattern  given  that  the  roadway  cut  is  ±400 
feet  from  center  of  missile  impact,  or  800  feet  for  the  total 
length.  This  represents  the  maximum  cut  possible  since  the 
major  axis  is  800  feet  long  and  8  vehicles  is  the  maximum 
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length  =  20  feet 
spacing  =  80  feet 


lmui 


number  that  could  fit  into  this  interval  given  their  length 
and  interval.  Using  an  interval  of  100  feet  (20  +  80  =  100), 
subsequent  calculations  are  made  to  determine  the  average 
number  of  vehicles  within  that  interval.  Another  assumption 
made  is  that  a  half  or  greater  portion  of  a  vehicle  within 
this  area  cut  is  also  considered  destroyed.  Table  2.2 
represents  the  interval  road  cuts  and  the  number  of  vehicles 
killed  within  that  distance. 

TABLE  2.2 


Road  Distance  Cut  and  Number  of  Kills — 
Minimum  Interval 


Distance  (feet) 

Kills 

Probabi 1 ity 

cut 

<  50 

0 

0.1074 

50  < 

cut  < 

150 

1 

0.0068 

150  < 

cut  < 

250 

2 

0.0118 

250  < 

cut  < 

350 

3 

0.0238 

350  < 

cut  < 

450 

4 

0.0338 

450  < 

cut  < 

550 

5 

0.0584 

550  < 

cut  < 

650 

6 

0.1052 

650  < 

cut  < 

750 

7 

0.2282 

750  < 

cut  < 

800 

8 

0.4246 

An  equation  for  an  ellipse  is  used  to  determine  the 
distance  of  the  y  cut  knowing  the  x  distance  cut. 

2  2 

Equation:  x  y_  _  - 

a2  b2  (2) 

where  x  =  value  of  cut  in  the  x  direction 
a  =  the  value  of  the  major  axis 
y  =  value  of  cut  in  y  direction 
b  =  the  value  of  the  minor  axis 
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This  equation  is  solved  for  y  by  taking  one  half  the  maximum 
x  value  in  each  interval.  This  y  value  is  then  normalized 
by  dividing  it  by  the  system  standard  deviation  of  93.43 
feet.  A  normal  distribution  table  is  then  used  to  integrate 
from  zero  to  this  normalized  value.  Doubling  this  value 
and  subtracting  it  from  the  value  obtained  for  the  previous 
y  cut  (explained  below),  gives  the  probability  of  killing 
that  many  vehicles. 

Two  exceptions  to  this  method  are  determining  the 
probability  of  zero  and  eight  kills.  Any  y  value  greater 
than  the  minor  axis  value  is  outside  effective  radius  of 
kill.  Dividing  the  minor  axis  by  the  standard  deviation 
gave  a  normalized  value.  Finding  the  probability  for  this 
value  in  the  normal  tables  and  doubling  it,  then  taking  its 
complement  gave  the  probability  of  zero  kills.  Once  the 
probability  of  zero  kills  was  known,  0.1074,  the  remaining 
probability  mass,  0.8926,  was  used  in  calculating  the  other 
probabilities.  For  example,  using  the  doubled  table  values 
of  one  kill,  0.8858,  and  subtracting  it  from  0.8926,  gave 
the  probability  of  killing  exactly  one  vehicle,  0.0068. 

Using  the  doubled  table  value  for  two  kills,  0.8740,  and 
subtracting  it  from  the  doubled  table  value  for  one,  0.8858. 
gives  a  probability  of  0.0118  for  killing  exactly  two  vehi¬ 
cles.  This  process  is  continued  to  determine  Table  2.2. 

The  complement  of  the  sum  of  the  number  of  zero  to  seven 
kills  then  gives  the  probability  of  eight  kills. 
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The  same  process  is  used  in  arriving  at  the  values 
in  the  bottom  half  of  Table  2.3,  the  case  for  the  maximum 


interval.  Instead  of  80  feet  between  vehicles,  a  160  feet 
interval  is  used.  A  random  number  drawn  from  a  uniform 
(0,1)  distribution  will  then  be  used  to  determine  destruc¬ 
tion  by  WAAM.  Vehicles  destroyed  are  then  subtracted  out 
of  the  rest  of  the  convoy.  A  new  convoy  length  is  then 
computed  and  the  delay  time  associated  with  this  attack  is 
determined. 


TABLE  2.3 

Road  Distance  Cut  and  Number  of  Kills — 
Maximum  Interval 


Distance  (feet) 

Kills 

Probability 

cut  <  170 

0 

0.1074 

170  <  cut  <  350 

1 

0.0424 

350  <  cut  <  530 

2 

0.0804 

530  <  cut  <  710 

3 

0.2290 

710  <  cut  <  800 

4 

0.5408 

Multiple  missile  attacks  are  treated  as  independent 
events  by  assuming  targeting  points  on  the  same  convoy  being 
at  least  four  standard  deviations  apart.  Treating  multiple 
missile  shots  this  way  maximizes  coverage  of  a  linear  target 
like  a  convoy. 

The  delay  time  associated  with  this  attack  is  dis¬ 


cussed  next . 


Delay  Due  to  Missile  Attack 

The  same  four  mutually  exclusive  delay  events 
associated  with  air  interdiction  also  apply  here.  Although 
the  road  is  not  cut  as  was  done  with  the  aircraft  interdic¬ 
tion,  damage  due  to  missiles  imposes  delay  due  to  more  likely 
destruction  of  vehicles  and  disruption  over  a  wider  area  due 
to  the  wide  area  coverage  of  the  WAAM . 

Delay  time  imposed  varies  as  a  function  of  number  of 
vehicles  killed  and  the  completed  actions  necessary  prior  to 
commencing  the  march.  This  model  assumes  that  the  missile 
attack  comes  at  the  end  of  the  impaired  movement  (IM)  time 
due  to  air  attack.  The  delay  associated  with  the  missile 
attack  is  added  to  the  old  delay  time. 

Imposed  delay  time  varies  as  a  function  of  number  of 
vehicle  kills.  If  five  or  less  vehicles  are  killed,  then 
delay  time  is  obtained  by  sampling  from  a  triangular  distri¬ 
bution  with  parameters  15,  18.5,  and  25.  For  kills  less 
than  or  equal  to  10  vehicles,  the  triangular  distribution 
parameters  are  22.5,  30.0,  and  37.5.  Finally,  if  the  number 
of  kills  is  greater  than  10,  the  distribution  parameters 
are  30,  41.5,  and  50.  These  parameters  were  obtained  by 
comparing  against  the  parameters  for  the  air  sortie  case. 
While  no  actual  distribution  exists  for  interdiction  with 
WAAM,  the  distributions  hypothesized  above  represent  the 
idea  that  as  more  targets  are  destroy;- d,  more  delay  time  is 
imposed  because  more  time  is  needed  for  clearing  roadways 
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of  destroyed  vehicles,  administering  first  aid,  and  restor¬ 
ing  order  to  continue  the  march.  Once  the  unit  accomplishes 
the  required  actions,  it  is  ready  to  continue  movement. 

Summary 

Chapter  II  has  described  how  the  problem  was  concep¬ 
tualized  so  that  computerization  can  begin.  An  MRD  was 
described  as  broken  down  into  44  units.  The  44  units 
traveling  on  the  hypothesized  three-route  network  represent 
the  steady  state  condition  of  the  system.  Route  interdic¬ 
tion  consisted  of  bombing  the  road  at  preselected  points  so 
that  the  convoys  are  blocked.  The  assumption  was  that  the 
aircraft  and  convoy  to  be  interdicted  arrived  simultaneously 
at  the  point  of  interdiction.  The  strike  mission  consisted 
of  a  flight  of  two  aircraft  each  making  one  pass  and 
dropping  all  12  500-pound  bombs  in  one  pass.  The  resulting 
stationary  convoys  would  then  increase  the  probability  of 
striking  the  convoy  with  one  or  more  missiles  loaded  with 
wide-area  anti-armor  munitions  (WAAM; . 

Chapter  III  follows  with  the  computerization  of  the 


Ill  Simulation  Model 

This  chapter  describes  how  the  system  structure  in 
Chapter  II  was  computerized.  The  description  includes  the 
convoy  movement  of  the  44  units  through  the  road  network,  the 
effects  of  interdiction  on  the  reconnaissance  element  on  each 
route,  and  the  subsequent  effects  on  the  arrival  time  for 
the  MRD.  The  model  structure  is  network  based  and  it  incor¬ 
porates  the  use  of  nodes  to  model  specific  geographic  fea¬ 
tures.  The  time  to  travel  between  nodes  represents  the  con¬ 
voy  travel  time  for  the  distances  represented  between  nodes. 

The  network  structure  of  STOPEM  and  its  interdiction 
process  were  both  modeled  using  Simulation  Language  for 
Alternative  Modeling  (SLAM).  SLAM  is  a  powerful  simulation 
language  that  provides  the  user  multiple  capabilities  to 
model  networks,  discrete  events,  and  continuous  events  or 
any  combination  of  these  three  processes.  Since  the  system 
structure  has  a  network  base  and  interdiction  processes, 
SLAM's  power  and  flexibility  provided  the  capability  to 
model  the  network  with  event  nodes.  The  interdiction  pro¬ 
cesses,  representing  changes  of  state  within  the  system, 
were  modeled  within  the  event  nodes  utilizing  user-written 
subroutines.  Another  feature  of  SLAM  allows  the  modeler  to 
assign  up  to  100  attributes  per  entity.  However,  STOPEM 
only  uses  16  attributes  per  entity  SLAM's  global  variables, 
XX(-)  with  their  values  known  to  all  subroutines  within  the 


program,  provided  the  capability  to  denote  such  events  as 
interdiction  and  to  act  as  indicator  variables.  One  capa¬ 
bility  not  exercised,  but  available  to  global  variables  is 
SLAM's  capability  to  provide  time  statistics  on  global 
variables  specified  by  the  user.  This  capability  is  above 
what  would  normally  be  associated  with  global  variables 
being  passed  within  user-written  subroutines.  Finally,  the 
waiting  time  for  interdicted  convoys  was  modeled  using 
AWAIT  nodes,  a  form  of  QUEUE  node.  The  few  SLAM  capabilities 
exercised  in  this  model  represent  but  a  small  portion  of 
SLAM's  overall  power.  The  reader  is  referred  to  Introduction 
to  Simulation  and  SLAM  (Ref  49)  for  a  more  complete  descrip¬ 
tion  of  the  language  and  SLAM's  capabilities. 

Network 

A  total  of  44  MRD  units  are  entered  into  the  road 
network  represented  by  the  Figures  3.1,  3.2,  and  3.3,  corres¬ 
ponding  to  Routes  BLUE,  GREEN,  and  RED,  respectively.  This 
represents  the  MRD  movement  on  the  hypothesized  three-route 
march  through  the  portion  of  the  corps  sector  being  examined. 
Route  BLUE,  Figure  3.1,  will  be  used  to  explain  the  computer¬ 
ization.  References  to  EVENT  and  ENTER  nodes,  and  other 
SLAM  structures  will  be  made  using  node  labels  associated 
with  each  structure.  This  label  is  located  immediately 


below  the  structure  in  the  referenced  figures.  Similar 
logic  was  used  on  Routes  GREEN  and  RED. 


Route  GREEN  (Continued) 


An  initialization  subroutine,  INTLC ,  is  used  to 
establish  the  units  for  the  model.  The  reader  is  referred 
to  Appendix  A,  the  SLAM  computer  code,  for  further  details. 
Fach  of  the  44  units  has  16  unique  attributes.  Attribute 
one,  the  unit  number  (same  number  that  appears  in  Table  2.1), 
is  used  to  rank  the  units  within  the  model  structure.  SLAM 
uses  a  system  of  files  to  maintain  order  as  specified  by 
the  user.  This  order  is  specified  as  a  ranking  within  the 
file.  The  entity  file  ranking  specified  for  file  one  was 
low  value  first  of  attribute  one,  the  unit  number.  This 
means  that  unit  one  with  all  its  associated  attributes  is 
first,  unit  two  is  second,  and  so  on  within  the  file.  This 
ranking  scheme  allowed  for  ease  of  accessing  the  units  out 
of  the  file  and  scheduling  them  for  movement  as  will  be 
explained  shortly.  Table  3.1  contains  a  description  of  the 
other  15  attributes.  Other  attributes  within  this  table 
will  be  explained  as  they  are  introduced. 

The  start  of  convoy  movement  on  the  road  network  is 
initiated  in  the  model  by  the  last  three  statements  in  INTLC. 
CALL  SCHDL  (1,  0.0,  ATRIB)  for  Route  BLUE  is  a  SLAM  state¬ 
ment  that  schedules  event  one  to  occur  at  time  of  0.0  or  at 
the  start  of  the  simulation.  The  last  argument,  ATRIB,  is  a 
vector  array  of  an  entity's  attributes  being  called  within 
that  event.  SLAM  dimensions  ATRIB  to  a  vector  (1,  100),  but 
the  model  only  has  a  (1,  16)  array  since  each  entity  has  16 


attributes.  All  40  events  used  in  STOPEM  are  found  within 


TABLE  3.1 


Entity  Attribute  Description 


Att  ribute 

Numbers 

Description 

1 

Unit  designator 

2 

Head/Tail  designation 
(1  =  Head,  2  =  Tail) 

3 

Number  of  track  vehicles 

4 

Number  of  wheeled  vehicles 

5 

Minimum  unit  length  (km) 

6 

Maximum  unit  length  (km) 

7 

Average  unit  length  (km) 

8 

Unit  length  (km) 

9 

Travel  time  for  first  node  (hrs) 

10 

Travel  time  for  subsequent  nodes  (hrs) 

11 

Day  travel  rate,  mixed  convoy  (km/hr) 

12 

Interval  between  convoys  (km) 

13 

Convoy  start  time  (hrs) 

14 

Convoy  finish  time  (hrs) 

15 

Event  code  (Assigned  by  SLAM) 

16 

Event  time  (Assigned  by  SLAM) 
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subroutine  EVENT.  The  schedule  statements  for  the  three 
routes  (1,  20,  and  32)  represent  the  event  number  that 
initiates  convoy  movement  for  that  route.  Within  each  event, 
the  first  unit  and  its  associated  attributes  are  referenced 
by  a  CALL  COPY  (entity,  file  number,  atrib)  statement  that 
copies  the  referenced  entity  with  associated  attributes 
from  the  appropriate  file.  For  Route  BLUE,  since  unit  one 
is  the  advanced  guard  for  that  regiment,  it  is  called  first 
from  file  one.  Since  the  time  keeping  mechanism  of  SLAM  is 
stopped  while  within  an  EVENT  node,  all  operations  within 
the  EVENT  node  are  of  zero  time  duration.  It  is  within  the 
EVENT  nodes  that  attributes  can  be  changed  to  model  system 
structure  changes. 

Evenc  one,  EVT1,  schedules  the  entry  into  the  net¬ 
work  for  the  head  and  tail  entities  for  convoy  one  and  the 
head  entity  for  convoy  two.  This  is  accomplished  by  the 
use  of  the  SCHDL  subroutine.  A  counter  variable,  II,  is 
used  within  event  one  to  keep  track  of  which  unit  out  of 
file  one  is  being  addressed.  Its  initial  value  is  assigned 
by  the  global  variable  XX(1),  initialized  to  one  and  incre¬ 
mented  by  one  before  the  event  termination.  Thus,  II  is 
initially  set  equal  to  one.  Within  GREEN  and  RED,  the 
counters  are  initialized  to  16  and  32,  respectively,  corres¬ 
ponding  to  the  lead  units  on  each  route.  To  regulate  convoy 
rates,  Subroutine  RESHUF1  is  called.  A  global  variable  XX(4) 
is  assigned  the  rate  value  for  unit  one.  XX(4)  is  then  used 


in  subsequent  calls  to  RESHUF1  to  compare  against  other  unit 
rates.  Based  on  whether  the  following  units  have  rates  that 
are  less  than  or  equal  to  or  greater  than  this  global  vari¬ 
able,  unit  rates  are  adjusted  up  or  down  to  match  the  lead 
unit  rate.  RESHUF2  and  RESHUF3  perform  the  same  function 
on  Routes  GREEN  and  RED.  This  assumes  that  all  units  within 
that  convoy  maintain  the  rate  set  by  the  lead  element.  The 
head  of  the  lead  unit  starts  the  convoy  movement  at  "time 
now"  (TNOW)  of  zero.  Its  entry  time  into  the  network  is 
placed  in  attribute  (13).  ENTER  node  ENT2  inserts  the  head 
entity  into  the  network.  As  its  name  implies,  the  ENTER 
nodes  functions  to  place  entities  into  the  network  at  the 
appointed  time.  The  tail  entity  for  the  first  convoy  is 
scheduled  to  enter  the  network  at  enter  node  ENT 3  at  the 
time  equal  to  the  passage  of  vehicles  between  the  head  and 
tail.  The  length  of  the  convoy,  attribute  (8),  is  divided 
by  the  rate,  attribute  (11),  to  obtain  the  travel  or  passage 
time.  Attribute  (8)  was  sampled  from  a  triangular  distribu¬ 
tion  composed  of  parameters  of  attributes  five,  six,  and 
seven.  This  calculated  time  is  added  to  the  time  of  entry 
for  the  head  entity,  which  is  zero  for  the  first  unit,  and 
used  to  schedule  entry  time  for  the  tail.  To  distinguish 
the  tail  from  the  head,  the  value  of  attribute  (2)  is  changed 
to  two.  The  time  it  enters  the  network  is  then  stored  in 
attribute  (13).  II  is  incremented  by  one  and  the  attributes 
of  unit  two  are  copied  into  the  array  ATRIB.  The  head 
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entity  of  the  unit  two  is  scheduled  to  arrive  at  the  time 
that  the  tail  of  one  entered  plus  the  time  between  the  tail 
and  itself.  The  value  of  attribute  (12)  for  the  first  unit, 
the  interval  between  the  two  units,  is  assigned  to  the  vari¬ 
able  CINT1  ( CINT2  and  CINT3  on  GREEN  and  RED,  respectively). 
This  distance  is  converted  to  time  by  dividing  CINT1  by  the 
adjusted  rate  of  attribute  (11)  for  unit  two.  RESHUF1  is 
called  again  before  this  operation  to  insure  the  rate  is 
adjusted  to  match  the  lead  rate.  The  time  is  then  added  to 
the  entry  time  for  the  tail  of  unit  one  to  obtain  the  entry 
time  for  head  two.  This  calculation  reflects  when  the  head 
of  unit  two  enters  the  road  network.  As  with  the  other  two 
entities,  its  entry  time  is  recorded  in  attribute  (13).  To 
continue  this  process,  the  scheduled  time  of  convoy  two 
triggers  event  one  again.  The  last  operation  before  exiting 
the  event  is  to  increment  XX(1)  by  one  so  that  the  next  time 
event  one  occurs,  the  counter  II  will  be  set  to  two  to 
schedule  the  tail  of  convoy  two  and  the  head  of  convoy 
three.  This  cyclic  process  continues  until  all  convoys  are 
moving  on  the  road  network. 

At  the  conclusion  of  each  event,  the  SLAM  event 
calendar,  which  keeps  record  of  when  to  schedule  entities, 
is  updated  with  the  unit  and  its  associated  attributes. 

When  TNOW  equals  the  value  of  that  entity  on  the  event 
calendar,  it  is  processed  at  the  next  SLAM  structure  if  the 
duration  activity  between  the  two  structures  is  zero.  The 


activities  between  SLAM  structures  represent  the  branches 
where  explicit  time  delays  are  modeled.  Otherwise,  if  the 
duration  activity  is  greater  than  zero  between  the  struc¬ 
tures,  the  entity  takes  that  long  to  get  to  the  next  struc¬ 
ture.  This  represents  the  convoy  movement  between  nodes. 

The  duration  between  structures  mentioned  above 
represents  the  time  in  the  model  associated  with  the  network 
travel.  This  aspect  is  modeled  by  dividing  the  distance  to 
travel  between  event  nodes  by  the  rate.  For  the  first  node, 
this  value  was  assigned  to  attribute  (9).  For  subsequent 
travel  distances,  the  travel  time  is  assigned  to  attribute 
(10).  Both  values  for  these  attributes  are  initialized  in 
INTLC,  but  recomputed  within  EVENTS  1,  20,  and  32  because 
convoy  march  rates  were  adjusted  in  these  events. 

Other  attribute  values  that  characterize  each  unit 
are  read  in  by  INTLC  for  appropriate  attribute  array  for 
all  units.  These  values  include  designation  for  the  head, 
and  the  number  of  wheeled  and  track  vehicles. 

Discrete  Events 

As  was  mentioned  previously,  SLAM  provides  the 
capability  of  modeling  the  system  changes  of  state  within 
each  event  node.  For  example,  within  selected  event  nodes 
user-written  Fortran  subroutines  are  used  to  model  system 
changes  such  as  interdiction  due  to  aircraft  and  missiles. 
Before  interdiction  is  described,  a  network  description  is 

i 
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given  of  the  system’s  steady  state.  The  next  section  then 
discusses  how  the  interdiction  processes  are  modeled. 

Steady  State.  Figure  3.1  is  the  description  of  the 
system  steady  state.  An  indicator  variable,  XX(10),  is 
used  to  denote  when  there  is  interdiction.  A  value  of  zero 
denotes  steady  state,  while  a  value  of  one  denotes  inter¬ 
diction.  As  shown  on  Figure  3.1,  entities  leaving  the  two 
enter  nodes  encounter  a  GOON  node  G01 .  This  GOON  node 
functions  as  a  continuation  node  between  EVT1  and  EV11. 

The  branch  or  activity  time  from  the  enter  nodes  to  G01  is 
zero  because  there  is  zero  travel  time  involved  between  the 
structures.  The  activity  time  between  G01  and  EV11  repre¬ 
sents  the  travel  time  for  convoys.  The  number  two  in  the 
box  beneath  the  branch  between  the  structures  is  the  activity 
number  that  uniquely  identifies  that  branch.  Once  at  EV11, 
representing  event  eleven,  a  conditional  statement  is  used 
within  the  event  node  to  recompute  attribute  (10)  since  it 
changed  from  the  value  computed  in  INTLC.  An  exiting  entity 
EV11  encounters  conditional  branching.  The  value  of  one  in 
each  EVENT,  ENTER  and  GOON  nodes  denotes  that  only  one  branch 
can  be  taken  on  the  exiting  end  of  the  structure.  This 
logic  models  the  convoys  taking  only  one  route.  Since  the 
system  is  at  steady  state,  that  is  XX(10)  equals  zero,  and 
all  units  on  Route  BLUE  are  to  go  on  this  branch,  that  is 
attribute  (1)  is  greater  or  equal  to  one,  only  one  branch 
can  be  taken  and  that  is  the  top  branch.  The  value  for 
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attribute  (10)  computed  in  EV11  is  used  for  the  duration 
from  EV11  to  EV14.  This  time  reflects  the  passage  time  for 
the  head  and  tail  entities  for  each  convoy.  The  same  pro¬ 
cess  continues  through  the  remainder  of  the  network.  Values 
for  the  convoy  passage  for  the  next  branch  are  computed  in 
the  preceding  event  node.  The  conditional  statement  at  EV16 
routes  the  entities  to  EV18  since  it  meets  the  top  condi¬ 
tion  first  and  only  one  branch  can  be  taken.  EV19  is  the 
last  event  for  Route  BLUE.  A  conditional  statement  within 
EV19  permits  only  the  tail  entity  of  each  convoy  to  enter 
the  enclosed  subroutines.  The  subroutine  calculates  total 
time  in  the  system  for  each  entity  by  subtracting  from  TNOW 
the  value  of  attribute  (13),  when  the  entity  entered  the 
road  network,  and  assigning  the  value  to  attribute  (14). 
Selected  values  are  sent  to  a  tape,  which  was  used  in  model 
verification,  discussed  in  Chapter  IV.  A  sample  of  these 
tapes  14  and  17  appear  in  Appendix  E.  Also,  event  19  main¬ 
tains  a  count  of  armored  vehicles  by  incrementing  XX(23) 
every  time  an  entity  enters  EV19.  When  the  last  unit 
finishes  (its  tail  entity),  its  finish  time  is  assigned  to 
XX(19).  This  represents  the  total  travel  time  for  all  units 
on  BLUE.  Subroutine  GAFMGF  is  then  called  to  increment 
XX(22)  by  XX(19)  so  that  XX(22)  represents  the  total  time  for 
all  three  routes  when  the  other  two  route  finish  times  are 
added  to  it  also.  For  the  other  two  routes,  this  occurs  at 
EV31  and  EV40 .  XX(28)  is  incremented  by  one  within  EV19 


and  used  as  a  test  within  GAFMGF  to  determine  when  the  last 
unit  of  the  MRD  has  finished.  Once  the  very  last  unit  has 
finished,  an  average  finish  time  for  the  MRD  is  obtained  and 
divided  into  the  total  number  of  armored  vehicles  to  arrive 
at  the  measure  of  effectiveness,  number  of  armored  vehicles 
per  unit  time.  This  value  and  other  values  (discussed  in 
Chapter  IV)  are  sent  to  an  output  tape.  A  sample  of  this 
output  tape  appears  in  Appendix  C. 

Upon  exiting  EV19,  conditional  branching  occurs  again. 
The  head  entities  take  the  bottom  branch  and  are  routed 
through  another  GOON  node,  GOE,  to  the  terminate  point  where 
the  entity  disappears  from  the  network.  This  action  corres¬ 
ponds  to  the  head  of  the  unit  arriving  at  its  release  point. 
The  tail  entity  takes  the  top  branch  and  is  routed  through 
two  COLLECT  nodes.  These  nodes  collect  statistics  on  the 
entities  passing  through  them.  The  collection  information 
is  a  function  of  the  argument  in  the  left  portion  of  the 
node.  For  example,  C0L1  collects  interval  statistics  on 
the  time  between  arrival  of  entities  based  on  attribute  (13). 
C01  collects  the  value  of  attribute  (3),  armored  vehicles, 
and  records  it.  At  the  end  of  each  simulation  run,  the 
statistics  collected  are  printed  in  table  format.  The  tail 
entity  travels  through  GOE  on  its  way  to  the  terminate  node, 
where  it  is  deleted  from  the  system.  The  tail  entity  arriv¬ 
ing  at  this  point  represents  the  completed  movement  of  each 
convoy . 
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Interdiction  Processes.  By  changing  the  value  of 
XX(10)  to  one,  the  system  structure  now  has  interdiction. 
Another  global  variable,  XX(ll),  denotes  the  number  of  mis¬ 
siles  targeted  against  convoys.  The  combination  of  these 
variables  denote  interdiction  and  the  levels  of  interdiction 
within  the  model.  Level  of  interdiction  is  defined  as  the 
combination  of  one  sortie  and  either  one,  two,  or  three 
missiles.  These  levels  will  be  discussed  further  in  the 
design  of  the  experiment. 

All  processes  discussed  in  the  steady  state  section 
pertaining  to  convoy  movement  and  calculation  of  travel 
times  apply  within  this  section  also.  The  difference  now  is 
that  the  system  structure  has  been  subjected  to  interdic¬ 
tion.  System  interdiction  is  modeled  only  for  the  recon¬ 
naissance  element.  STOPEM  as  it  is  now  does  not  incorporate 
targeting  the  main  body.  The  model  interdiction  occurs  at 
EVENT  nodes  12,  23,  and  34  for  Routes  BLUE,  GREEN,  and  RED, 
respectively.  The  basic  process  within  each  event  is  similar 
One  difference  in  EVT1  of  BLUE  is  that  attribute  (9)  is 
recomputed  when  therg  is  interdiction  because,  when  the 
reconnaissance  element  is  interdicted  and  the  main  body 
bypasses  the  point  of  interdiction,  the  distance  used  in 
INTLC  has  changed  from  64.784  kilometers  to  50  kilometers. 
Since  convoy  one  continues  on  its  original  route  when  it  is 
interdicted,  EV11  computes  the  remaining  distance  time 
originally  computed  in  INTLC  for  attribute  (9) 


(64.784  -  50  =  14.784).  This  value  divided  by  the  rate  is 
then  assigned  to  attribute  (10).  This  attribute  (10)  is 
the  time  duration  from  EV11  to  EV12.  EV11  also  calculates 
the  new  adjusted  attribute  (10)  for  the  main  body.  This 
new  time  represents  the  time  to  bypass  the  new  distance  and 
is  the  time  duration  from  EV11  to  EV13. 

Conditional  statements  within  EVENT  11  insure  that 
the  appropriate  values  for  attribute  (10)  are  assigned  to 
the  correct  units.  Conditional  branching  on  the  exit  end 
of  EV11  only  allows  unit  one,  both  head  and  tail  entity, 
to  be  routed  to  EV12.  The  main  body  of  the  convoy  is  routed 
to  EV13  using  the  bottom  branch.  This  structure  simulates 


main  body  bypass  of  the  interdicted  unit  at  EV12.  This 
process  simulates  Warsaw  Pact  doctrine  of  pressing  forward 
with  the  movement.  The  main  body  continues  on  through  EV14 , 
EV15,  EV16 ,  and  EV17.  This  process  portrays  the  main  body 
moving  through  the  road  network  without  interdiction  and  is 
similar  to  the  steady  state  description.  The  main  body  will 
be  discussed  again  at  the  conclusion  of  the  discussion  of 
interdiction  at  EV12. 

As  the  middle  branch  exiting  from  EV11  shows,  only 
unit  one  can  take  this  branch  when  there  is  interdiction. 

The  activity  duration  models  the  time  until  the  head  entity 
is  interdicted  at  EV12.  A  conditional  statement  at  the 
beginning  of  EV12  insures  that  only  the  head  entity  is 
allowed  to  enter  the  interdiction  subroutine.  Once  within 


v  n 


r 


EV12,  subroutine  SORTIE  is  called  and  a  two-ship  strike 
mission  is  simulated  against  the  head  of  the  convoy.  In 
SORTIE,  the  variables  DTIME  and  DTIMEP  represent  a  sample 
from  a  triangular  distribution  denoting  full  or  partial  road 
blockage,  respectively.  SORTIE  then  draws  a  sample  from  a 
uniform  (0,1)  distribution  and  assigns  this  random  variate 
to  the  variable  PASS.  If  PASS  is  less  than  or  equal  to 
0.35,  then  the  road  is  considered  cut.  If  PASS  is  greater 
than  0.35,  then  partial  road  blockage  is  obtained.  A 
counter  variable,  XX(9),  is  incremented  by  one  if  the  road 
is  cut  and  zero  if  the  road  is  partially  blocked.  Another 
sample  is  taken  from  the  uniform  distribution  to  denote  the 
second  aircraft  pass.  The  same  test  is  performed  on  PASS  to 
determine  if  it  is  less  than  or  equal  or  greater  than  0.35. 
SORTIE  returns  a  value  for  XX(9)  of  zero,  denoting  partial 
road  blockage,  or  a  value  greater  than  zero  denoting  that 
the  road  was  cut.  This  value  of  xx(9)  is  then  tested  to 
determine  how  much  delay  time  was  imposed  as  a  result  of 
interdiction.  Attribute  (3)  is  decreased  by  two  to  denote 
two  armored  vehicle  kills.  The  event  time  is  increased  by 
the  delay  time  imposed,  which  is  either  DTIME  or  DTIMEP. 

The  event  time  represents  when  the  entity  will  be  scheduled 
on  the  network  again.  SLAM  places  the  entity  event  time 
in  the  attribute  that  is  two  greater  than  the  maximum  number 
of  attributes  for  the  model.  Since  the  maximum  number  of 
attributes  is  14,  the  event  time  is  stored  in  attribute  16. 


Since  entity  one  has  suffered  damage  and  delay,  the  number 
of  armored  vehicles  is  adjusted  and  its  event  time  is 
increased  by  an  appropriate  amount  as  sampled  from  SORTIE. 

Since  the  assumption  made  was  that  the  delay  time 
due  to  missile  attack  would  be  added  to  the  delay  time  due 
to  missile  attack  would  be  added  to  the  delay  due  to  SORTIE, 
subroutine  MISSILE  is  called  next.  A  delay  time  correspond¬ 
ing  to  the  number  of  kills  is  sampled  from  a  triangular  dis¬ 
tribution  and  assigned  to  the  variables  DTMIS1,  DTMIS2,  and 
DTMIS3,  where  the  order  of  variables  indicate  an  increasing 
time  delay.  The  number  of  missiles  shot  is  assigned  to 
MSHOT  by  XX(ll).  Either  one,  two,  or  three  missile  salvos 
can  be  fired  at  the  target.  The  number  of  missiles  fired 
will  be  a  function  of  interdiction  policy,  resource  con¬ 
straints,  and  target  priority.  In  order  to  compute  target 
kills,  the  spacing  within  vehicles  is  required  to  be  computed 
for  each  unit.  The  value  of  attribute  (8)  is  decreased  by 
the  value  represented  by  the  number  of  companies  within  each 
unit  times  the  average  interval  (37.5  meters)  between  units. 
This  value  is  then  divided  by  the  total  number  of  vehicles 
in  that  unit  and  converted  to  feet.  MISSILE  then  uses  a 
looping  structure  to  simulate  firing  from  one  to  three 
missiles.  For  each  loop  passage,  a  random  sample  drawn  from 
the  uniform  (0,1)  distribution  is  assigned  to  SHOT.  Then 
based  on  the  vehicle  interval  computed  outside  the  loop,  the 
subroutine  tests  to  see  if  the  vehicle  interval  was  less 


than  or  equal  to  120  feet  or  greater  than  120  feet .  The 
variable  SHOT  is  then  tested  within  the  appropriate  branch 
and  the  number  of  kills  is  assigned  to  the  variable  NUMKIL. 
NUMKIL  is  incremented  for  each  missile  shot.  MISSILE 
returns  to  event  12  the  number  of  kills  and  the  delay  time 
imposed.  The  number  of  vehicle  kills  is  subtracted  from 
attribute  (3)  and  the  event  time  is  incremented  by  the  delay 
time.  If  the  number  of  armored  vehicles  is  less  than  zero, 
attribute  (4)  (wheeled  vehicles)  is  decreased  by  this 
amount  to  model  the  correct  number  of  kills.  The  head 
entity  exits  the  event  and  is  filed  on  the  event  calendar. 
Once  the  event  calendar  schedules  it  again,  the  delay  time 
has  transpired  and  it  travels  to  AWAIT  node  MAF,  with  dura¬ 
tion  equal  to  attribute  (10). 

The  lower  portion  of  event  12  contains  the  condi¬ 
tional  statenv"  ~  for  the  tail  entity.  Its  values  for  attri¬ 
butes  three  and  16  must  be  adjusted  to  correspond  to  the 
head  entity  values.  Two  variables  in  the  head  entity  block, 
VAL1  and  DELCH1,  are  assigned  the  head's  values  for  attri¬ 
butes  16  and  three.  These  variables  are  then  assigned  to 
the  tail  entity's  attributes  three  and  16.  Additionally, 
attribute  (16)  is  increased  by  the  separation  distance 
between  the  head  and  tail.  This  is  accomplished  by  dividing 
attribute  (8)  by  the  rate.  Once  this  is  accomplished,  the 
tail  entity  exits  the  event  and  is  filed  on  the  event 
calendar.  When  it  appears  again  on  the  event  calendar,  it 


travels  to  AWAIT  node  MAF,  with  duration  equal  to  attri¬ 
bute  (10). 

AWAIT  node  MAF  stores  convoy  one  in  file  two  while 
unit  one  waits  for  GATE  BL0K1  to  open.  SLAM's  concept  for 
a  gate  is  as  its  name  implies.  A  GATE  is  either  opened 
or  closed  and  either  allows  or  denies  passage,  respectively. 
In  this  model,  the  initial  gate  status  for  BL0K1  is  closed 
as  denoted  in  the  GATE  block  to  the  left  of  EVT1 .  This 
SLAM  structure  shows  that  BL0K1  is  initially  closed  to 
model  that  the  main  body  will  not  stop  or  slow  up  to  allow 
interdicted  units  to  re-enter  between  units  of  the  main  body. 
This  action  is  modeled  in  this  fashion  because  allowing  the 
interdicted  unit  to  re-enter  would  cause  congestion  on  the 
main  body  or  a  slowing  down  that  might  present  a  better 
target  owing  to  units'  interval  decrease.  The  ranking  of 
file  two  is  based  on  low  value  first  of  attribute  two. 

This  implies  that  the  head  entity  is  filed  before  the  tail 
entity.  Opening  BLOK1  is  modeled  by  the  tail  entity  of  the 
last  unit,  unit  15.  The  exit  end  of  EV17  shows  the  top 
conditional  branch  being  reserved  only  for  the  tail  entity 
of  unit  15,  the  last  unit  in  the  main  body.  The  structure 
at  STR1  is  an  OPEN  node  that  opens  BLOK1.  This  action  models 
the  tail  entity  of  15  clearing  the  road  junction  ior  the 
blocked  unit  to  enter.  No  interval  distance  is  modeled 
between  units  15  and  one  to  denote  that  unit  one  enters 
immediately  behind  15.  The  bottom  branch  exiting  EV17  is 


85 


AD-fll39  495  STOPEM  fl  SIMULATION  INTERDICTION  MODEL  OF  A  MOTORIZED  2/3 
RIFLE  DIVISIONS  AIR  FORCE  INST  OF  TECH 
HR  I GHT -PATTERSON  AFB  OH  SCHOOL  OF  ENGINEERING 
UNCLASSIFIED  G  fl  FULTON  NflR  83  AFIT/GST/0S/83H-2  F/G  15/7  NL 


taken  by  all  units,  including  the  head  of  15  and  unit  one. 
The  units  then  maintain  this  order  and  finish  with  unit  one 


last.  The  remainder  of  the  network  for  BLUE  is  as  described 
in  the  steady  state  case. 

Route  GREEN  is  similar  to  BLUE  because  an  AWAIT  node 
EKF  is  also  modeled  there.  Both  routes,  even  if  just  inter¬ 
dicted  once,  cause  convoys  to  meet  at  some  future  point  or 
take  a  much  longer  bypass  to  prevent  this  meeting.  This 
particular  portion  of  the  road  network  creates  disruption 
and  bunching  of  units.  Even  though  the  road  network  is 
seemingly  infinite,  interdicting  these  two  routes  shows  an 
instance  where  convoy  disruption  can  be  maximized. 

STOPEM  only  does  one  interdiction  of  the  head 
element  of  the  reconnaissance  unit.  An  examination  of  the 
effects  of  interdiction  for  this  road  network  on  the  main 
body  would  create  more  lucrative  targets  for  WAAM  targeting, 
with  increased  probabilities  of  destruction. 

On  the  other  hand,  Route  RED  is  an  autobahn  that  if 
interdicted  offers  the  possibility  for  quick  exiting  and 
entering  again  at  some  future  point.  Route  RED  models  that 
the  reconnaissance  element  is  interdicted  and  the  main  body 
exits  and  bypasses  the  obstacle.  As  on  the  other  two 
routes,  interdiction  occurs  and  damage  and  delay  are  assessed 
However,  instead  of  waiting  while  the  main  body  gets  back  on 
the  main  route,  unit  one  exits  the  route  prior  to  the  point 
where  the  main  body  enters.  Event  node  36  represents  unit 

86 


one’s  exit  and  event  node  38  models  the  main  body’s  entry 
back  onto  the  original  route.  The  difference  between 
Routes  BLUE  and  GREEN,  and  Route  RED  is  that  there  are  more 
likely  exit  points  on  RED  heading  in  the  original  direction 
of  travel  that  do  not  cause  bunching  as  opposed  to  the 
restrictions  on  BLIT]  and  GREEN. 

Summary 

The  entire  MRD  broken  down  into  44  units  with 
associated  attributes  is  initialized  into  file  one.  Once 
in  the  file,  SLAM  enter  and  event  nodes  schedule  the  entry 
of  convoys  through  the  road  network.  Each  convoy  is  modeled 
as  two  entities  with  unique  attributes  except  for  attribute 
two.  This  attribute  denotes  whether  the  entity  is  the  head 
or  tail  element.  The  indicator  variable  XX(10)  denotes 
steady  state  when  set  equal  to  zero  and  interdiction  when 
set  equal  to  one.  The  flexibility  of  the  model  is  given  by 
event  nodes  where  user-written  subroutines  model  changes  of 
state  for  the  system.  When  XX(10)  is  set  equal  to  one, 
interdiction  occurs  in  selected  event  nodes  on  each  route. 
Within  this  event,  subroutines  SORTIE  and  MISSILE  are  called 
to  determine  interdiction  effects.  At  the  conclusion  of 
each  convoy  arrival  at  the  event  node  for  that  route,  data 
is  collected  for  analysis  by  calling  subroutine  GAFMGF. 

The  above  concludes  Chapter  III  and  how  model  com¬ 


puterization  was  accomplished.  Chapter  IV  is  next  with  an 
explanation  of  validation  and  verification. 


IV 


Model  Verification  and  Validation 


The  model  verification  and  validation  process  is  an 
important  and  necessary  area  to  which  the  modeler  must  sub¬ 
ject  the  model  if  he  expects  to  have  an  acceptable  level  of 
confidence  in  the  inferences  drawn  about  the  model  outputs 
of  the  system.  This  chapter  will  discuss  the  process  of 
verification  and  validation  used  for  STOPEM. 

Shannon  in  Systems  Simulation  divides  the  process  of 
evaluating  simulations  into  three  categories: 

1.  verification  -  insuring  that  the  model  behaves 
the  way  an  experimenter  intended. 

2.  validation  -  tescing  the  agreement  between  the 
behavior  of  the  model  and  that  of  the  real 
system. 

3.  problem  analysis  -  the  drawing  of  statistically 
significant  inferences  from  the  data  generated 
by  the  computer  simulation  (51:30). 

Chapter  VI  will  discuss  the  problem  analysis,  while  the 

remainder  of  this  chapter  will  discuss  the  first  two  points. 

Verification 

Two  categories  of  tests  were  performed  to  verify  the 
internal  consistency  of  the  model.  The  first  involved 
statistical  testing  of  distributions  to  determine  their 
behavior  in  the  model.  The  second  category  of  testing  con¬ 
sisted  of  using  the  SLAM  trace  option,  print  statements,  and 
data  sent  to  tapes  to  verify  that  the  model  and  its  activi¬ 
ties  were  performing  as  desired. 


Sortie  and  .Missile  Distributions.  Appendix  E  con¬ 
tains  the  SPSS  output  on  testing  of  the  Uniform  (0,1)  dis¬ 
tribution  used  in  SORTIE  and  MISSILE.  The  null  hypothesis 
tested  using  a  Kolmogorov-Smirnov  test  was  that  the  samples 
generated  came  from  the  distribution  indicated.  The  null 
hypothesis  failed  to  be  rejected. 

Triangular  Distributions.  No  statistical  testing 
was  conducted  on  the  triangular  distributions.  The  idea 
used  to  develop  the  distribution  was  based  on  the  low  and 
high  values  given  in  the  references  cited.  A  mode  value 
was  computed  from  these  values  and  sampling  conducted  during 
the  simulations.  Checks  made  throughout  model  development 
using  the  SLAM  trace  option,  explained  next,  showed  that 
the  values  obtained  were  within  the  specified  ranges. 

Trace  Monitoring.  Utilizing  SLAM's  trace  option 
and  print  statements  throughout  the  program  allowed  each 
entity  and  four  other  attributes  to  be  followed  from  event 
to  event.  Since  the  trace  option  also  showed  TNOW  and  dura¬ 
tion  to  next  event,  checks  were  made  to  determine  if  enti¬ 
ties  were  arriving  at  their  next  point  at  the  appointed 
time.  This  option  also  gave  the  terminate  time  for  each 
entity.  By  knowing  when  a  unit  entered  the  system  and  when 
it  arrived  at  its  destination,  finish  time  data  was  checked 
for  accuracy.  Print  statements  such  as  shown  in  Appendix  E 
also  gave  feedback  if  interdiction  subroutines  were  working 


as  modeled.  In  all  cases  checked,  the  model  was  functioning 
as  intended. 

Validation 

Three  possible  tests  exist  for  examining  model 
validity.  The  first  is  face  validity.  This  test  consists 
of  checking  the  model  at  extreme  values  for  parameters  and 
examining  for  the  reasonableness  of  results.  The  latter 
method  is  referred  to  as  a  Turing  test.  This  test  consists 
of  using  experts  to  examine  model  outputs  to  compare  these 
outputs  against  the  actual  system.  The  other  two  methods 
consist  of  'testing  of  assumptions"  and  "the  testing  of 
input-output  transformations"  (Ref  51:29).  These  last  two 
methods  are  concerned  with  statistically  related  testing  for 
such  things  as  tests  of  means  and  variances  and  analysis  of 
variance  (51:29).  STOPEM’s  validation  process  consisted  of 
a  combination  of  the  first  and  the  last  of  the  above  tests. 
The  usage  of  these  tests  in  STOPEM's  validation  is  one  of 
degree,  as  will  be  discussed  after  an  explanation  for  the 
purpose  of  the  model  is  made. 

STOPEM  represents  an  abstraction  of  the  real  system. 
To  balance  the  reality  of  the  system  under  study  with  the 
limitations  of  time,  resources  such  as  computer  availability, 
and  tractability  of  the  model,  assumptions  discussed  in 
previous  chapters  were  made  to  formulate  STOPEM  as  it  is 
now.  To  simulate  the  complete  system  under  study  would 
require  a  model  of  such  complexity  that  the  required  time 


to  exercise  it  would  consume  incalculable  resources  and  time 


As  an  abstraction,  STOPEM  represents  a  portion  of  the  system 
with  those  variables  perceived  to  be  more  important,  that 
when  acted  upon  by  external  forces,  present  insights  about 
the  system  and  its  interactions,  if  any.  The  validation  of 
STOPEM  must  be  evaluated  with  this  in  mind. 

Variation  of  the  parameters  for  the  triangular  dis¬ 
tributions  was  not  conducted  at  the  extreme  values.  How¬ 
ever,  a  ten  percent  reduction  in  the  mode  parameters  for 
the  convoy  rate,  length,  and  interval  was  made.  These 
variations  will  be  discussed  further  in  Chapter  V.  The 
results  given  based  on  these  variations  proved  to  be  consis¬ 
tent.  Results  will  be  discussed  in  Chapter  VI.  The  sub¬ 
routines  SORTIE  and  MISSILE  contain  concepts  of  the  contri¬ 
butions  of  air  sorties  and  missile  strikes.  The  modeling 
of  air  sortie  results  in  the  program  in  Appendix  B  and  the 
discussion  of  the  reasonableness  of  these  results  with 
informed  experts  gave  credibility  to  the  parameter  used  in 
determining  when  the  road  was  cut  or  partially  blocked. 

The  effects  of  a  missile  attack  was  hypothesized  based  on 
information  received  from  Fort  Sill,  Oklahoma.  The  generic 
missile  postulated  and  the  effects  gave  consistent  results. 
Finally,  the  uniform  (0,1)  distribution  used  in  MISSILE  and 
SORTIE  was  tested  as  discussed  in  the  verification  portion 
of  this  chapter.  All  the  above  combined  give  STOPEM 
validity . 
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Once  the  model  verification  and  validation  is  com¬ 
plete,  testing  of  the  system  can  begin.  This  chapter  has 
explained  how  the  model  was  verified  and  validated. 

Chapter  V  discusses  the  data  collection  and  how  the  experi 
ment  was  conducted. 


In  order  to  draw  meaningful  inferences  about  the 
system  under  study  once  the  model  has  been  verified  and 
validated,  an  experimental  design  is  needed  to  run  the  model 
to  obtain  the  data  necessary  to  make  an  analysis  (51:145). 
Designs  of  experiments  have  two  purposes:  "(1)  they  are 
economical  in  terms  of  reducing  the  number  of  experimental 
trials  required;  and  (2)  they  provide  a  structure  for  the 
investigator's  learning  process"  (Ref  51:145).  This  chapter 
discusses  how  the  experiment  was  designed  to  obtain  the 
required  data. 

Measure  of  Effectiveness  (MOE) 

In  this  thesis,  a  MRD  consisting  of  44  units  was  put 
through  the  system  in  each  simulation  run.  A  measure  of  effec 
tiveness  (MOE)  was  developed  to  use  as  a  standard  to  compare 
model  outputs.  The  MOE  developed  is  the  arrival  rate  per 
hour  of  armored  vehicles  of  the  MRD  at  its  destination. 

This  MOE  was  selected  as  a  way  to  measure  the  arriving  com¬ 
bat  power,  represented  by  armored  vehicles,  because  one 
objective  of  this  thesis  was  to  determine  how  the  arrival  of 
combat  systems  can  be  affected  by  interdiction.  The  more 
the  MRD  can  be  delayed,  disrupted,  or  destroyed,  the  less 
the  impact  on  the  friendly  force  at  the  MRD  point  of  commit¬ 
ment.  The  design  of  the  experiment,  which  follows  next, 


relates  how  the  model  experimentation  will  be  conducted  to 
determine  the  MOE  or  response  variable. 
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Exprimental  Design 

Shannon  in  Systems  Simulation  describes  experimental 
design  as  a  systematic  three-way  process.  The  three  steps 
are:  "(1)  design  of  the  structural  model;  (2)  design  of 

the  functional  model;  and  (3)  design  of  the  experimental 
model"  (Ref  51:151).  The  structural  model  and  the  functional 
differ  in  that  the  former  presents  what  should  be  done  in 
terms  of  factors  and  factor  levels  to  be  analyzed,  while 
the  latter  says  what  can  be  done.  This  difference  is  attrib¬ 
uted  to  constraints  such  as  computer  time  or  funds  that  might 
require  a  reduction  of  factors  or  levels  in  the  experimental 
design.  This  redesign  becomes  the  functional  model.  A 
definition  follows  for  factors,  levels,  and  cells  in  the 
context  used  above  and  in  the  remainder  of  this  research. 

A  factor  is  an  independent  variable  being  examined,  while 
the  factor  level  refers  to  a  factor  value.  Table  5.1  pre¬ 
sents  the  factors  and  associated  levels  analyzed.  These 
will  be  discussed  further  in  subsequent  paragraphs.  A  cel' 
refers  to  the  "basic  structure  or  building  block  of  an 


experiment"  (Ref  51:153).  j 

Equation  5.1  was  used  to  determine  the  structural  \ 


model : 


Ns  =  ((li^cl2^q3)  •  ’  * 
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«  ‘  V*  N  .S  .N  .*»  ^ 


(5.1) 


where 


N  =  number  of  cells  in  the  experiment 
s 

k  =  number  of  factors  in  the  experiment 

=  number  of  levels  of  ith  factor,  i=l,2,3,...k 

(Ref  51:153) 

Based  on  this  formula,  the  structural  model  for  this  experi¬ 
ment  is  (4)(2)(2)(2)  =  32  cells.  Because  computer  time  was 
not  a  constraint  and  the  model  run  \me  was  not  excessive, 
each  cell  had  a  response.  This  m;  »  the  functional  model  a 
complete  model  as  opposed  to  an  i  >mplete  model  that  has 
fewer  responses  than  the  number  oi  ells. 

Sample  size  determination  for  each  cell  or  the  number 
of  replications  per  cell  was  based  on  what  is  considered 
adequate.  An  adequate  number  of  replications  per  cell  is 
eight  (Ref  51:163).  To  determine  where  the  response  vari¬ 
able  mean  would  lie  based  on  this  number  of  replications, 

Eq  (5.2)  was  used  to  make  a  calculation  to  determine  this. 

(0Za)2 
n  *  g— " 

a  (5.2) 


where 

n  =  number  of  runs 

a  =  population  standard  deviation 

Z  =  two-tailed  standardized  normal  statistic 
a 

2 
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ci  =  alpha  level 
d  =  allowable  error 

Using  a  =  0 . 05  ,  n  =  8  ,  and  d  =  —r  ,  the  above  equation 

X 

gave  that  the  response  for  the  sample  mean  would  lie  between 
U  ±  0.70  ,  where  u  is  the  unknown  true  population  mean. 

This  implies  that  using  eight  replications  per  cell  gives  a 
95  percent  probability  that  the  interval  y  ±  0.70  contains 
the  sample  mean  for  the  response  variable. 

In  order  to  capture  the  model  response  due  to  the 
variation  of  these  four  factors  and  their  associated  levels, 
a  full  factorial  design  was  run.  A  full  factorial  design 
implies  that  "all  levels  of  a  given  factor  are  combined  with 
all  levels  of  every  other  factor  in  the  experiment"  (51:164). 
This  design  allows  for  the  study  of  the  factor  main  effects 
as  well  as  interactions  between  factors  if  they  exist. 
Interactions  refer  to  both  the  individual  influence  of  a 
factor  and  the  combined  influence  of  two  or  more  factors  on 
the  response  variable.  Also,  since  the  effect  of  each  fac¬ 
tor  is  estimated  at  several  levels  of  the  other  factors, 

"the  conclusions  reached  hold  over  a  wide  range  of  condi¬ 
tions"  (51:165). 

The  factors  and  the  levels  used  in  this  experiment 
were  structured  into  logical  patterns  so  that  the  necessary 
observations  could  be  made.  As  previously  mentioned, 

Table  5.1  contains  the  factors  and  levels.  The  subsequent 
paragraphs  explain  in  more  detail  the  factors  and  levels. 
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Factor  one,  the  interdiction  policy,  was  set  at 


four  levels  to  denote  no  interdiction  and  the  combination 
of  one  sortie  plus  either  one,  two,  or  three  missiles. 

These  four  levels  represent  four  possible  interdiction 
policies,  with  increasing  number  of  missiles  denoting 
increased  need  for  destruction  and  delay  or  a  higher  priority 
target . 

Convoy  rate,  the  second  factor,  was  set  at  two 
levels:  sample  from  a  triangular  distribution  (20,25,30) 

and  sample  from  another  triangular  distribution  (20,22.5,30). 
The  first  rate  represents  the  distribution  with  the  mode  as 
the  average  of  the  minimum  and  maximum  rates,  whereas  the 
second  represents  a  -ten  percent  shift  in  the  mode  to  repre¬ 
sent  a  partial  decrease  due  to  such  effects  as  road  damage 
and  weather. 

The  third  factor,  convoy  length,  was  also  considered 
at  two  levels:  the  first  level  represents  the  length  as 
drawn  from  another  triangular  distribution  with  the  param¬ 
eters  left  as  computed  in  the  initialization  subroutine. 

The  second  level  represents  a  -ten  percent  decrease  in  the 
mode  parameter  for  each  entity.  This  represents  the 
decrease  in  convoy  length  due  to  mechanical  vehicle  losses 
or  prior  interdiction. 


The  fourth  factor,  convoy  interval,  represents  the 
distance  in  kilometers  between  convoys.  As  with  the  other 
factors,  the  first  level  represents  the  initial  computation 


of  parameters,  while  the  second  level  represents  decreased 
convoy  interval. 

Table  5.2  is  the  design  matrix  for  the  blocks  of 
runs.  As  this  table  shows,  a  total  of  256  simulation  runs 
were  made.  These  256  runs  w^ere  divided  into  blocks  of  64 
runs,  where  there  were  8  cells  within  each  block,  and  8 
replications  per  block,  for  a  total  of  64. 


TABLE  5.2 

Design  Matrix  for  Block  of  Simulation  Runs 


Run  Number 

Policy 

1-64 

One  (No  interdiction) 

65-128 

Two 

129-192 

Three 

193-256 

Four 

Table  5.3  shows  how  each  block  from  Table  5.2  was 
designed.  This  table  shows  that  runs  one  through  eight 
represent  the  first  cell  with  all  factors  set  at  level  one. 
The  second  set  of  runs,  nine  through  16,  show  that  all  the 
factor  levels  remain  the  same  except  that  the  interval  now 
is  at  level  two.  This  same  process  continues  for  the  other 
six  cells.  For  runs  65  through  128,  129  through  192,  and 
193  through  256,  Table  5.3  would  be  the  same,  except  that 
the  levels  for  the  policy  would  be  two,  three,  and  four, 
respectively . 


TABLE  5.3 


r 

t 


Desi gn 


Matrix 


Within  Each 


Block 


Run 


Factors 

Policy  Rate  Length  Interval 


1-8  1 

9-16  1 

17-24  1 

25-32  1 

33-40  1 

41-48  1 

49-56  1 

57-64  1 


1 

1 

1 

1 

2 

2 

2 

2 


1 

1 

2 

2 

1 

1 

2 

2 


1 

2 

1 

2 

1 

2 

1 

2 


The  experiment  was  run  once  the  measure  of  effective 
ness,  the  appropriate  sample  size,  and  the  experimental  de¬ 
sign  were  determined.  The  analysis  of  this  data  is  pre¬ 
sented  in  Chapter  VI. 


VI 


Data  Analysis 


The  next  step  after  gathering  data  from  the  experi¬ 
mental  design  is  to  conduct  an  analysis  of  the  data.  This 
part  of  the  analysis  process  allows  the  analyst  to  make 
valid  inferences  relative  to  the  model's  output.  Analysis 
of  variance  (ANOVA)  using  the  Statistical  Package  for  the 
Social  Sciences  (SPSS)  is  an  available  tool  to  test  hypoth¬ 
eses  to  determine  if  several  populations'  means  are  equal. 
ANOVA  provides  the  analyst  with  a  capability  to  determine 
if  response  variable  means  are  equal  based  on  the  factors 
and  levels.  The  analyst  can  then  make  inferences  about  the 
model  outputs. 

Both  four-way  and  three-way  ANOVAs  were  conducted 
on  the  data  generated  by  the  design  of  experiment  discussed 
in  Chapter  V.  The  terms  three-way  and  four-way  refer  to 
three  and  four  factors  in  the  experiment.  Both  ANOVAs  were 
also  conducted  with  all  second  and  higher  order  interactions 
suppressed.  This  implies  that  these  ANOVA  tables  presented 
in  Appendix  C  do  not  have  interaction  effects  present.  The 
interactions  were  suppressed  to  gain  further  insights  about 
the  main  effect  contributions  since  all  second  and  higher 
order  interactions  proved  to  be  insignificant.  Appendix  C 
contains  all  SPSS  outputs  and  tape  1,  which  is  the  actual 
data.  A  discussion  follows  of  the  results. 
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Three-Way  and  Four-Way  ANOVAs 


Both  ANOVAs  with  interactions  showed  that  only  two 
of  the  main  effects  were  significant  using  an  alpha  level  of 
0.05.  All  two-way  and  higher  interactions  proved  to  be 
insignificant  at  this  same  alpha  level.  The  factor  length 
was  deleted  to  make  the  three-way  ANOVA .  The  reason  for 
its  deletion  was  based  upon  it  having  the  lowest  value  of 
significance  based  on  the  four-wav  ANOVA  with  interactions. 

The  lower  F  value  pointed  to  its  contributing  less  to  explain¬ 
ing  the  main  effects.  The  ANOVAs  without  interactions  show 
the  same  two  main  effects  being  significant.  In  order  to 
gain  a  more  robust  representation  of  the  statistical  contri¬ 
butions  of  the  main  effects,  a  selection  was  made  from  the 
two  ANOVAs  without  interactions  of  the  one  with  the  lowest 
mean  square  error.  The  model  thus  selected  would  tend  to 
explain  more  about  the  main  effects.  The  four-wav  ANOVA 
proved  to  have  the  lowest  mean  square  error.  This  is  reason¬ 
able  considering  that  all  factors,  in  spite  of  two  being 
insignificant,  make  for  a  more  complete  system  representation 
since  the  convoy  length  and  convoy  interval  represent  key 
parts  of  the  system. 

Main  Effects.  Main  effects  refers  to  each  indepen¬ 
dent  variable  or  factor.  The  SPSS  output  in  Appendix  C  shows 
that  only  two  of  the  four  main  effects  were  considered  sta¬ 
tistically  significant  at  the  alpha  level  of  0.05.  This 
says  that,  at  the  stated  alpha  level,  the  mean  value  for  the 
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response  variable  is  statistically  different  for  different 
levels  of  the  factors  of  policy  and  rate.  To  graphically 
illustrate  that  this  statistical  difference  does  exist,  a 
plot  is  made  with  the  main  effects  and  their  associated 
levels  along  the  horizontal  axis  and  the  response  variable 
along  the  vertical  axis.  The  lines  drawn  between  points 
only  serve  to  emphasize  changes  in  the  response  variables 
due  to  factor  levels  and  do  not  represent  linear  relation¬ 
ships.  As  the  plots  in  Figure  6.1  show,  the  change  in  slope 
in  policy  and  rate  indicate  that  a  difference  does  exist 
between  levels  of  these  two  factors.  In  contrast,  the  other 
two  factors  show  relatively  no  change  in  slope  and  hence 
show  no  statistical  differences.  This  implies  that,  for  the 
given  levels,  the  two  factors  do  not  have  any  significant 
impact  on  the  response  variable. 

Main  effects  behave  as  expected.  As  policy  level 
one  shows,  no  interdiction  allows  a  rate  of  arrival  of 
approximately  80  armored  vehicles  per  hour  at  the  destination 
point.  The  interdiction  effects  show  up  as  the  steep  change 
in  slope  of  Figure  6.1  when  going  from  no  interdiction  to 
one  sortie  and  one  missile.  Policies  two,  three,  and  four 
show  no  appreciable  difference.  As  would  be  expected, 
policy  four  was  a  combination  of  one  sortie  and  three  mis¬ 
siles  and  it  has  the  lowest  vehicle  arrival  rate.  A  bigger 
difference  between  policies  two,  three,  and  four  might  have 
been  expected  if  the  main  body  of  the  convoy  had  been 
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interdicted  also.  The  other  factor,  convoy  rate,  shows  a 
drop  in  vehicle  arrival  rate  as  the  rate  is  decreased  ten 
percent.  This  is  again  consistent  with  what  would  be 
expected  because  a  reduced  convoy  rate  means  units  arrive 
at  the  destination  later.  The  other  two  main  effects,  while 
not  significant,  also  show  consistency.  As  the  length  and 
interval  go  from  level  one  to  two,  the  arrival  rate  increases. 
Since  level  two  is  a  ten  percent  decrease  in  the  convoy 
length  and  interval,  the  units  are  neither  as  long  nor  is 
the  distance  to  travel  between  them  as  far.  This  implies  that 
a  lesser  distance  takes  less  time  to  travel,  so  the  arrival 
rate  increases. 
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Fig  6.1.  Influence  of  Main  Effects 

Two-Way  and  Higher  Interactions.  As  the  SPSS  print¬ 
out  in  Appendix  C  shows,  there  are  no  significant  two-way  or 
higher  interactions.  An  explanation  offered  for  this  lack 
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of  interactions  is  the  measure  of  effectiveness  (MOE) 
variable.  Since  the  MOE  has  time  as  part  of  its  measure, 
it  is  likely  that  the  rate  and  policy  factors  dominate  the 
other  two  or  that  the  amount  of  variation  (ten  percent 
reduction)  was  not  sufficient  to  get  a  bigger  change  in  the 
response  variable  for  length  and  interval.  A  sensitivity 
analysis  conducted  at  the  extreme  values  for  these  factors 
might  provide  some  additional  insights  into  how  the  model 
is  working.  Another  possible  explanation  for  this  lack  of 
two-way  interactions  and  that  an  almost  significant  three- 
way  interaction  of  length,  rate,  and  interval  was  close  to 
the  0.05  cutoff,  is  that  the  equal  variance  assumption 
between  cells  might  have  been  violated.  A  test  for  equality 
of  variance,  which  was  not  conducted,  might  show  that  this 
is  the  case.  Since  the  number  of  replications  per  cell, 
eight,  was  not  large  enough  to  invoke  the  central  limit 
theorem  (CLT),  the  replications  per  cell  would  have  to 
increase  to  meet  CLT  criteria  of  30  or  greater  per  cell. 

The  CLT  says  that  as  the  number  of  replications  approaches 
infinity  the  random  variable  response  approaches  normality 
(Ref  51:187).  Since  the  above  are  explanations  that  were 
not  tested,  they  are  offered  as  suggestions  for  follow-on 
work  or  checking  of  the  model  outputs. 

This  chapter  has  explained  the  data  analysis  of  the 
experiment.  Chapter  VII  presents  the  Results  and  Conclu¬ 
sions  . 
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VII  Conclusions  and  Recommendations  and 
Recommendations  for  Follow-on  Studv 


This  thesis  modeled  a  motorized  rifle  division 
(MRD),  represented  by  44  units,  moving  through  a  hypothesized 
three-route  network.  The  experimentation  consisted  of  com¬ 
paring  the  steady  state  of  the  system,  that  is  no  interdic¬ 
tion,  against  the  three  policies  of  interdiction,  consisting 
of  combinations  of  sorties  and  missiles.  The  measure  of 
effectiveness  was  the  arrival  rate  of  armored  vehicles. 

The  analysis  consisted  of  a  four-way  ANOVA  conducted  on  the 
data  generated  by  the  levels  of  the  four  factors.  This 
chapter  presents  the  conclusions  and  recommendations  and 
recommendations  for  follow-on  study. 


Conclusions  and  Recommendations 


The  objective  of  this  thesis,  as  stated  in  Chapter  I, 
was  to  develop  a  model  of  a  MRD  moving  through  a  road  net¬ 
work  to  determine  how  interdiction  would  delay  the  MRD  com¬ 
bat  power  arrival  rate  into  a  portion  of  a  corps  sector. 

The  conclusions  are  as  follows: 

1.  Interdiction  of  just  one  element  on  each  route 
does  cause  delav,  disruption,  and  destruction 
of  the  MRD. 

2.  The  difference  in  increased  levels  of  missile 
usage  did  not  appreciably  show  a  significant 
decrease  in  the  combat  arrival  rate  as  Figure  6.1 
shows.  This  result  is  attributed  to  the  target¬ 
ing  of  only  reconnaissance  units  since  the 
rerouting  of  the  main  body  onto  a  different 
route  to  bypass  the  interdicted  point  did  not 
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add  appreciable  delay  time  to  the  MRD.  This  is 
due  in  part  to  a  road  network  that  permits 
relative  ease  of  bypassing  the  interdicted 
point . 

3.  Modeling  the  interdiction  of  roads  and  convoys 
by  conventional  munitions  is  a  difficult  task 
that  requires  close  coordination  between  the 
Army  and  the  Air  Force. 

Based  on  the  above  conclusions,  the  following 
recommendations  are  made: 

1.  That  the  model  be  enhanced  further  to  include 
interdiction  of  the  main  body  to  gain  insights 
into  further  decreasing  the  combat  power 
arrival  rate. 

2.  That  efforts  to  improve  coordination  between 
the  Air  Force  and  the  Army  continue  to  insure 
the  greatest  effect  of  employment  of  limited 
resources  against  the  second-echelon  threat. 


Recommendations  for  Follow-on  Study 

Since  this  research  effort  could  not  address  all 
aspects  of  the  system  studied  or  provide  insights  to  all 
the  questions  that  need  to  be  asked,  some  areas  are  recom¬ 
mended  for  thought  or  follow-on  study  concerning  second- 
echelon  interdiction. 

As  was  explained,  STOPEM  assumed  a  hypothesized 
route  for  the  MRD.  Developing  a  model  that  would  determine 
shortest  route  distance  based  on  time  to  travel  would  pro¬ 
vide  ideas  on  which  routes  would  be  more  likely  travel 

* 

routes.  Such  a  model  would  necessitate  a  decision  structure 
that  would  incorporate  type  of  road  (autobahn,  primary, 
secondary,  or  tertiary);  the  load  carrying  capacity  of  any 
bridges  along  that  route  that  would  limit  route  to  vehicles 


less  than  a  certain  weight  classification;  a  route  that 
would  facilitate  convoy  passage  (whether  a  winding  road, 
hilly,  numerous  defiles);  road  interdiction  potential;  and 
the  number  of  built-up  areas  to  traverse. 

The  effect  of  interdiction  and  convoy  passage  will 
degrade  the  road  quality  over  time.  As  time  into  the  con¬ 
flict  increases,  convoy  rates  would  decrease  as  a  result. 

The  hypothesized  triangular  distribution  for  the  rate  then 
might  better  be  represented  by  a  mode  that  is  closer  to  the 
low  value  parameter  instead  of  the  average  of  the  low  and 
high  values. 

This  model  assumes  that  the  head  element  of  the 
convoy  is  always  interdicted.  This  process  needs  to  be 
modified  to  include  the  stochasticity  of  when  and  where  the 
sortie  and  convoy  would  actually  meet  at  the  point  of  inter¬ 
diction.  The  model  assumes  perfect  intelligence  to  say  that 
existing  technology  and  new  technology  will  facilitate  tar¬ 
get  acquisition.  Also,  sortie  generation  is  considered 
constant,  with  no  modeling  of  the  FEBA  penetration. 

One  aspect  not  incorporated  is  the  refugee  problem  and 
how  it  would  affect  the  movement  of  both  combat  forces  and 
the  logistical  effort.  Until  the  refugee  problem  subsides, 
it  can  be  expected  to  have  a  lowering  effect  on  convoy 
movement  rates. 

Development  of  an  interactive  program  to  model  real 
time  movement  of  convoys  would  provide  a  valuable  training 


108 


tool  for  training  of  corps  and  division  staff  elements  in 
the  process  of  getting  a  feel  for  how  rapidly  second-echelon 
forces  will  close  the  gap  from  detection  to  commitment. 

Such  a  program  would  show  the  importance  of  getting  timely 
intelligence,  processing  it,  deciding  whether  or  not  to 
target,  ordering  the  mission,  executing  the  mission,  and 
analyzing  the  result.  Such  a  model  might  show  that  incoming 
data  could  easily  overwhelm  the  information  processors,  and 
thus  delay  the  order  to  interdict,  resulting  in  a  missed 
opportunity.  This  overwhelming  of  the  processors  would  be 
due  to  the  amount  of  incoming  information  and  the  arrival 
rate  of  this  data. 

Since  the  Warsaw  Pact  Armies  are  expected  to  consume 
prodigious  quantities  of  both  petroleum  products  and  ammuni¬ 
tion  with  their  rates  of  advance  and  reliance  on  firepower, 
an  indirect  approach  to  getting  at  this  one  weakness  might 
be  to  examine  what  the  required  rates  are  to  maintain  this 
momentum  and  how  interdiction  of  these  factors  would  affect 
the  capability  to  fight  or  sustain  an  offensive. 

STOPEM  assumes  that  echelonment  of  forces  is  taking 
place.  While  the  technique  of  echelonment  is  not  charac¬ 
teristic  in  all  attack  scenarios,  attack  without  echelonment 
would  mean  a  dispersal  of  forces  along  a  wider  frontage. 
Intelligence  nets,  processing  of  information,  and  decision 
making  would  be  deluged  with  inputs.  This  would  make  it 
difficult  to  stop  the  attacking  forces  conventionally 
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because  the  enemy  force  density  would  be  fairly  uniform, 
arriving  at  the  FEBA  simultaneously.  However,  this  might 
raise  the  U.S.  nuclear  threshold  because  of  presented  target 
densities.  This  might  also  increase  the  case  for  more  and 
more  lethal  vide  area  weapons  to  stop  this  wide  attack.  One 
key  point  here  is:  do  not  assume  away  the  possibility  of 
an  attack  on  a  wide  frontage  if  Warsaw  Pact  intention  is  to 
overwhelm  a  dispersed  defense  and  prevent  lateral  shifting 
of  forces.  Also,  if  the  Warsaw  Pact  perceives  NATO's  reluc¬ 
tance  to  use  nuclear  weapons,  the  Warsaw  Pact  could  use  high 
speed  roads  such  as  autobahns  to  double  up  convoys  abreast 
and  get  units  to  the  front  faster.  Without  nuclear  weapons, 
conventional  forces  would  be  strained  to  stop  such  rapid 
movement.  As  before,  such  a  movement  tactic  could  be 
countered  with  WAAM  type  weapons.  The  point  of  all  this  is 
that  any  studies  done  need  +  remain  cognizant  that  echelon- 
ment  and  dispersal  a~e  not  fixed  doctrine  because  circum¬ 
stance  could  dictate  exceptions  to  what  is  considered  normal 
operating  procedure. 

STOPEM  assumes  that  the  interval  between  vehicles 
and  units  varies  between  25  and  50  meters.  A  study  to  deter¬ 
mine  Warsaw  Pact  equipment  operator  driving  habits  and  how 
they  affect  intervals  and  march  discipline  might  provide 
some  interesting  insights  of  how  accurate  estimates  are.  One 
would  expect  that  little  or  no  previous  driving  experience 
of  the  average  Warsaw  Pact  soldier  would  mean  stricter 


110 


compliance  with  intervals.  The  contrast  with  the  western 
counterpart  who  has  more  driving  experience  and  thus 
carries  bad  habits  (failure  to  maintain  interval)  with  him 
when  operating  military  vehicles  while  in  convoy  would 
probably  show  that  the  Warsaw  Pact  armies  have  better  march 
discipline.  This  characteristic  might  make  it  harder  to 
create  vehicle  bunching. 

The  areas  mentioned  in  this  chapter  are  considera¬ 
tions  for  further  study  with  STOPEM  or  as  separate  efforts. 
Some  are  related  directly  and  others  indirectly  since  second- 
echelon  interdiction  is  a  small  subsystem  of  a  much  larger 
and  complex  system.  The  ideas  presented  in  this  chapter 
were  meant  as  much  to  stimulate  thought  as  well  as  possibly 
provide  ideas  for  enhancement  of  STOPEM  or  creation  of  new 


models  or  studies. 
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|  APPENDIX  A 


SLAM  COMPUTER  MODEL 


The  computer  model  starts  on  the  next  page  and  is 
divided  into  three  parts:  (1)  job  control  cards;  (2)  the 
SLAM  coded  portion;  and  (3)  the  FORTRAN  coded  portion.  The 
SLAM  portion  has  an  explanation  of  the  global  variables  and 
the  attributes.  The  FORTRAN  coded  portion  has  comment  cards 
at  strategic  points  to  further  clarify.  The  user  is  cau¬ 
tioned  about  numbering  of  enter  nodes.  SLAM  is  limited  to 
25  enter  nodes.  Thus,  any  enter  nodes  numbered  greater 
than  25  will  result  in  a  SLAM  input  error.  A  way  to  get 
around  this,  as  this  model  shows,  is  to  number  any  enter 
nodes  with  the  beginning  numbers  of  the  network,  that  is 
all  network  enter  nodes  are  numbered  one  through  eight. 


CAF,C‘<  165000.T400, 1 0400.  T«2 038° ,  FULTON ,  42 SR 

ATTACH , PR  OC  FI L , I n=A » 10 ] 7 1 , S N=A SDAB . 

BEGIN, NOFRILE. 

ATTACi: ,  PROCFIL,  SI.AMFROC ,  I  n=AFIT. 

GET,THF22 , ID=FULTON. 

GE T ,  TUS  2  2 , 1  D»FU  LTON . 

FTN5(I=TUS22,ANSI=0). 

PEG  IN' ,  SLAMI I , ,  I=THF,22 ,  M=LGO . 

FYIT,U. 

REWIND, TAPE  11. 

REWIND, TAPE  15. 

REWIND, TAPE  ]R. 

REWIND.TAPE23. 

PE  PLACE, TAPE  1 1 , ID=FULTON. 

RE PLACE, TAPE  15, ID-FULTON. 

REPLACE, TAPE  19 , ID*FULTON. 

RF.  PLACE ,  T A  PF.  2  3 ,  ID-FULTON . 

* 

* 

* 

* 

* 

* 

* 

GEN.C.A.  FULTON, STOPF.M, 03/ 10/83,256; 

LIMITS, 3,  14,300; 

niIORITY/l,LVF(  1); (INITIAL  FILE  WITH  44  ENTITIES) 
PRIORITY/2, LVF(2); (AWAIT  NODE  PLOKl) 

PRIORITY/3, LVF(2); (AWAIT  NODE  BL0K2 ) 

NETWORK; 


********************************************************************* 

DESCRIPTION  of  KEY  GLOBAL  VARIABLES  AND  ATTRIBUTES  FOLLOW.  EACH  OF 
THE  THREE  ROAD  NETWORKS  IS  OUTLINED  BELOW.  INITIALLY, ATTRIBUTES 
9  AND  10  ARE  READ  INTO  FILE  ONE.  SUBSEQUENT  VALUES  FOR  ATTRIBUTE  9  &  10  ARE 
RECOMPUTED  DUE  TO  RATE  CHANGE  IN  SUBROUTINE  EVENT.  GATES  BLOK  1  6  BLOK 2 
ARE  INITIALLY  CLOSED  TO  DENOTE  THAT  THE  MAIN  BODY  HAS  PRIORITY  AT  A  ROAD 
JUNCTION  &  THAT  THE  UNIT  ENTERING  HAS  TO  WAIT. 

********************************************************************* 


XX(  1)«C0UNTER  TO  ACCESS  ATTRIBUTES  FILE  KENTITIF.S  1-15) 
XX(2)»C0l’!:rER  TO  ACCESS  ATTRIBUTES  FILE  1(F.NTITIES  16-3  1) 
XX( 3 ) “COUNTER  TO  ACCESS  ATTRIBUTES  FILE  1(ENTITIES  32-44) 
XX(  4  )*VAP.TARLF.  FOR  TESTING  RATE  ,  RESHUF  1 
XX( 5 )“VARI ABLE  FOR  TESTING  RATE , RESHUF2 
XX(6)-VARIABLE  FOR  TESTING  RATE , RESHUF3 


XX(7)=NUMRER  OF  SORTIES 
XX(8)=NUMRER  OF  MISSILES 

XXr  0)=COUNTF.R  IN’  SORT  IF.  FOR  NUMBER  OF  KILLS 

XX( 10)=VARIABLE  TO  INDICATE  V.’HETHEP  OR  LOT  INTERDICTION^  —  NO;  1--YFF) 
XX(  1  !)=COUNTER  FOR  NUMBER  OF  MISSILES  SI’.OT(  l,2,OR,3) 

XX( I2)=VARIABLE  TO  INDICATE  THAT  EVENT  12  HAS  OCCURRED 
(0— NO;  I— YES), INITIALIZED  TO  0.0  IN  INTLC 
XX(  13)=VAP.IARLE  TO  INDICATE  THAT  EVENT  23  HAS  OCCURRED 
(0— NO;  1— YES), INITIALIZED  TO  0.0  IN  INTI.C 
XX(  ]4)*VARIARLF.  TO  INDICATE  THAT  EVENT  34  HAS  OCCURRED 
(0 — NO;  1 — YF.S)  .INITIALIZED  TO  0.0  IN  INTLC 
XX( 15)»VARIABLE  TO  ADJUST  MODE  PARAMETER  IN  CONVOY  LENGTH 
DISTRIBUTIONS — NONE ; 0.  3—10%  REDUCTION) 

XX(  16  )=VARIARl,E  TO  ADJUST  MODE  PARAMETER  IN  CONVOY  RATE 
DISTRIBUTION (2 5.0 — REGULAR;  22.5 —  3  0%  REDUCTION) 

XX(  17)=VARIABLE  TO  ADJUST  MODE  PARAMETER  IN  CONVOY  INTERVAL 
DISTRIBUTION— RECON  UNITS  ONLY(L’NITS  1,16,32) 

( 2 5 .0—  REGULAR; 22.5—10%  REDUCTION) 

XX(  18)=VARIABLE  TO  ADJUST  MODE  PARAMETER  IN  CONVOY  INTERVAL 
DISTRIBUTIONS— ALL  OTHER  UNITS 
(4.0— REGULAR; 3. 6— 10%  REDUCTION) 

XX(  19)=TIME  OF  FINISH  FOR  CONVOYS  ROUTE  #  1 

XX(20)=TIME  OF  FINISH  FOR  CONVOYS  ROUTE  il  2 

XX(2  1)«TIME  OF  FINISH  FOR  CONVOYS  ROUTE  if  3 

XX(22)*(XX(  19)+XX(20)+XX(2  1) )/3 — TIME  OF  TRAVEL  MRD 
XX(23)=ARMORF.D  VEHICLES  FINISHED  ROUTE  #  1 
XX(24)=ARMORED  VEHICLES  FINISHED  ROUTE  ft  2 
XX(25)=ARMORED  VEHICLES  FINISHED  ROUTE  *  3 

XX(26)=XX(23)+XX(24)+XX(25)— TOTAL  ARMORED  VEHICLES  FINISHED 
XX( 2 7 )=MEASUEE  OF  MERIT— NUMBER  OF  ARMORED  VEHICLES  ARRIVED  AT 
TERMINATION  POINT  rER  HOUR— XX(26)/XX(22) 

XX(28)«COUNTEF  USED  IN  OUTPUT  SUBROUTINE  GAFMCF 
ATRIB( 1)=UNIT  TYPE 

ATRIB(  2 ) :  1=HEAD  OF  CONVOY,  2=TAIL  OF  CONVOY 
ATRIB(3)=N0.  OF  TRACK  VEHICLES 
ATRIB( 4 )=N0.  OF  WHEEL  VEHICLES 
ATR IB( 5 )=MIN  UNIT  LENCTH  (KM) 

ATRIB( 6 )*MAX  UNIT  LENGTH  (KM) 

ATRIB( 7 )=AVG  UNIT  LENGTH  (KM) — ( ATRIB( 5 )+ATRIB( 6 )/2 ) 

ATRIB( 8 )*UNIT  LENCTH, SAMPLE  FROM  TRIANGULAR  DIST  (KM) 

ATRIR(9 )=TRAVEL  TIME  FOR  1ST  NODE  (DISTANCE/ATRIB( 1 1)— HRS) 

ATRIB(  10)=  TRAVEL  TIME  SUBSEQUENT  NODES  (DISTANCE/ATRIB(  1  1)— HRS) 
ATRIB(  11)=DAY  RATE, MIXED  CO N\’0Y, SAMPLE  FROM  TRIANGULAR  DIST  (KM/HR) 
ATRIR( 12)=INTERVAL  LENGTH, SAMPLE  FROM  TRIANGULAR  DIST  (KM) 

ATRIB(  13)=TIMF.  ENTITY  STARTS  ( T NON ) 

ATRIB( 14)=TIME  ENTITY  FINISHES  (TNOU-ATRIB(  13) ) 

ATRIB(  15)-ASSIGNED  BY  SLAM  AS  **AX  NUMBER  OF  ATTRIRl’TES(  14)+  ] , 

AS  EVENT  CODE 

ATRTB(  16)=ASSICNF.D  BY  SLAM  AS  MAX  NUMBER  OF  ATTRIBUTES (  14)+2, 

EVENT  TIME 
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;roete  blue (entities  1-15) 

» 

CATF/BLOK 1, CLOSE, 2; 

CATF/BL0K2 , CLOSE, 3 ; 

» 

EVT 1  EVENT,],]; 

ACT2  ACT,,, CO]; 

FN'Tl  ENTER,  2,]; 

ACT, , ,C01; 

ENT  2  ENTER,  3,1; 

ACT  1  ACT /  1 , , ,GO 1 ; 

FNTR  ENTER, R, 1; (NOT  USED) 

ACT, ATRIB( 10) , , MAF ; 

C0  1  GOON,  1; 

ACT 3  ACT/2 ,ATRIB( 9) ; 

FV11  EVENT,  11,1; 

ACT4  ACT/3, ATRIB( 10),XX( 10) .EO.O.O.AND.ATRIR( 1).CE.  1.0.EV14; 
ACT, ATRIB(  lo),XX(  10). EC.  1. O.AND. ATRIB(  1).E0.  1.0.EV32; 
ACT, ATRI B(  10),XX(  10). EO.  1.0.  AN'D.ATRIIU  1).CT.  1.0.EV13; 

EV  12  EVENT,  12,1; 

ACT/4,ATRIR( 10), , MAF; 

MAF  AWAIT(2),BL0K1, 1; 

ACT/5, ,,EV 17; 

EV 13  EVENT,  13,1; 

ACT/6; 

FV 14  EVENT,  14,1; 

ACT/7, ATRIR(  10); 

FV 15  EVENT, 15,1; 

ACT, ATRIB(  10); 

FV  16  EVENT,  16,1; 

ACT, ATRI B( 10),XX(  10) .EO.0.0.AND.ATRIR(  1).CE.  1.0, EV  IF; 
ACT , ATRIB(  10)  ,XX(  10). EO.  1 .0 . AMD. ATRI B(  1).CT.  1.0, F.V  17; 

EV  1 7  EVENT,  17,1; 

ACT, ATRTB( 10),XX( 10). EQ. 1.0. AND.ATRIB( 1).E0. 15.0. AND. 
ATRIR(2).EC.2.0,STR1; 

ACT,ATRIB(  10)  ,XX(  10). EO.  1. 0.  AND.  ATRI  B(  1)  .LE.  15.0,  F.V  IP; 
STR  1  OPEN,  BLOK  1,  1; 

ACT,  ,  ,EV  18; 

F.V18  EVENT,  IP.,  1; 

ACT, ATRI R(  10); 


********************************************************************* 

HEAD  OF  THE  CONVOY  RYFASSFS  COLLECT  NODE .  THIS  NODE 

ONLY  INTERESTED  IN  TATL  TIME  FINISH  OF  EACH  ENTITY  AND  THEN 

LAST  I'NIT  FINISH  TIME. 
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EV  19  EVENT,  1°,]; 

ACT 5  ACT, ,ATRIB(2).EQ. l.O.GOE; 

ACT 6  ACT; 

COE  1  COLCT,IN'T(  13)  .UNIT  TIME  RT"  1 ; 
ACT; 

COl  COLCT, ATRIB( 3 ) .ARMORED  VEHICLES; 

ACT, , ,COE; 


•ROUTE  CREEK (ENTITIES  16-31) 

9 

EV20  EVENT, 20,1; 

ACT, , ,G0 1 1; 

ENT4  ENTER, 4,1; 

ACT, , , GO  1 1 ; 

ENT 5  ENTER, 5,1; 

ACT/8, ,,0011; 

ENT9  ENTER, 9 , 1 ; (NOT  USED) 

ACT, ,,EKF; 

C0  11  GOON,  1; 

ACT/<1,ATRIB(9); 

EV21  EVENT, 2 1,1; 

ACT /  10, ATRIR( 10); 

EV22  EVENT, 22,1; 

ACT/  1 1,  ATP.IB(  10),XX(  10) .FO.O.O.AND.ATRIR(  1).CE.  16.0.EV25; 
ACT,ATRIB( 10),XX( 10). EO. 1 .O.AND. ATRIB( 1).E0. 16.0.EV23; 
ACT,ATRIB( 1"),XX(  10). EQ. 1 .O.AND. ATRIB( 1).GT.  16.0.EV24; 
EV23  EVE NT, 23,  1; 

ACT/ 12 , ATRIR( 10),,EKF; 

EKF  AWAIT( 3 ) , BL0K2 , 1 ; 

ACT/ 13, ,,EV26; 

EV24  EVENT, 24,1; 

ACT/  14, , ,EV25; 

EV25  EVENT, 2 5,1; 

ACT/  15; 

EV26  EVENT, 26,1; 

ACT,  ATRIB(  10),XX(  10)  .EO. 0.0. AND.  ATRIB(  1).CE.  16.0,F.V2R; 
ACT, ATRIB( 10),XX( 10). EQ. 1.0.AND.ATRIR(  1) .EC. 3 1.0. AND. 
ATRIB(2).E0.2.0,STR2; 

ACT, ATRIB( 10),XX( 10). EQ. 1 .O.AND. ATRIB( 1).GE. 16.0.EV27; 
RTR2  OFEN.BLOK2, 1; 

ACT, , ,EV28; 

FV27  EVENT, 27, 1 ; 

ACT, ATRIB( 10),XX( 10).E0. 1.0.ANP.ATRIB( 1).CE. 16.0.EV2R; 
EV28  EVENT, 28,1; 

ACT/  16,ATRIB( 10); 


FV2°  EVENT, 29,  1; 

AFT , ATR I R ( 10); 
r.v3o  eve't,3o,i; 

ACT, ATR IP ( 30); 


HEAD  0^  THF.  CONVOY  BY  TASSOS  COLLECT  NODE.  THIS  NODE 

ON'l.Y  INTERESTED  IN  TAIL  TIME  FINISH  OF  EACH  ENTITY  AND  THEN 

THE  LAST  UNIT  FINISH  TIME. 

********************************************************************** 


FV31  EVENT, 3 3,1; 

ACT,  , ATRIB(  2).E0.3.0,  COF, ; 

ACT; 

C0L2  COLCT,INT(  13)  ,l!NIT  TIME  RT/-‘2  ; 
ACT; 

CO 2  COLCT, ATR IR( 3), ARMORED  VEHICLES; 

ACT,,, COE; 


RO!TE  RED(ENTITIES  32-44) 


EV32  EVE  NT, 32,1; 
ACT, , ,C02 1; 


ENT 6  ENTER, 6,  3; 

ACT,  ,  ,C02  1; 

ENT 7  ENTER, 7, 3; 

ACT/  17, , ,C02  1; 

EN10  ENTER, 10, 1; (NOT  USED) 

ACT, , , G03Q; 

CO 2 3  COON ,  1; 

ACT /  IP , ATRIR( 9 ) ; 

F.V33  EVENT, 33,1; 

ACT, A7RIR( 10),XX( 10) .EO. 0.0. AND. ATRIB( 1 ) .CE . 32 .0, EV36; 
ACT,  , ATRIR(  1)  .EO. 32.0. AND. XX(  30).F.Q.  1.0.EV34; 

ACT,  ATRIR(  10)  ,  ATRIR(  1)  .CT.  32 .0.  AND.XXC  30). EO.  1.0.F.V35; 
EV34  EVE  NT,  34,1; 

ACT, , , 0030; 

CO30  COON,  1; 

ACT, ATRIR(  10), ,EV36; 

EV35  EVENT, 35,1; 

ACT; 

EV3A  EVE  NT, 36,1; 

ACT,  ATRTR(  10),XX(  10)  .F.O  .0 .0 .  AND.  ATRIP  (  3 )  .CE.  32 .0  ,F.V37 ; 
ACT,  ATRIB(  10),XX(  10). EO.  1 .0.AN’D.ATRIR(  1 )  .CT .  32 . 0 ,  EV3R ; 


ArT,  ATRTP-(  10),XX(  ]0).r.0.  l.P.AND.ATPTR(  ]  )  .EC .  32 .0 ,  F.V39 ; 
FV37  EVENT, 37,1; 

AFT/  19,  ATRIP.f  10); 

FV3R  EVENT,  38,1; 

ACT/20, ATRIB( 10) , ,EV40; 

FV  39  EVE NT, 39,1; 

ACT/2  1,ATRTR( 10) , ,FV40; 


HEAD  OF  THE  CONVOY  BYPASSES  COLLECT  NODE.  THIS  NODE 

ONLY  INTERESTED  IN  TAIL  TIME  FINISH  OF  EACH  ENTITY  AND  THEN 

THE  LAST  UNIT  FINISH  TIME. 

*********************************************************************** 


FV40  EVENT, 40,1; 

ACT, ,ATRIB(2).E0. 1.0, COE; 

ACT; 

C0L3  COLCT,  I  NT/  13) ,  UNIT  TIME  PT/(3  ; 
ACT; 

C03  COLCT, ATRIR( 3), ARMORED  VEHICLES; 

ACT,,, COE; 


COE  COON,  1 ; 

TE  TERM ; 

END; 

TNIT , 0 ,  150; 

;  START  OF  256  SIMULATION  RUNS,  DONE  IN  BLOCKS  OF 
;  8,  RANDOM  NUMBER  SEEDS  REINITIALIZED  AT  START  OF 

;  EACH  BLOCK  OF  64  RUNS. 

» 

INTLC  ,XX(  10)=0.0,xxm)=0.0,xx(  15)=0.0,XX(  16)=2  5.0,XX(  17)=25.n, 
XX(  18)-4.0; 

SEEDS; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

SIMTLATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX(  10)-0.0,X\7  1 1)«0.0,XX(  15)=0.0,XX(  16)=25.0,XX(  17)=22.5, 
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Vv( 

SEEPS; 

S  IMULATE ; 

SIMULATE ; 

SIMULATE; 

SIM I  LATE ; 

SIMULATE ; 

SIMULATE ; 

SIMULATE; 

SIMULATE; 

INTLC ,  XXI  10)=0.0,XX(  n)=0.0,XX(  IS)=0.n,XX(  Ifi  )=2  2  .  A ,  XX(  ]7)  =  2  5.F 
XX( 1 8 ) =4 .0; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

IUTLC.XXI  10)-n.0,XX(  n)=0.n,XXM5)=0.0,XX(  16)=22.5,XXI  37)=2 2.5 
XX(  1R)  =  3 .ft; 

SEEPS ; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XXI 10)=P.o>xx(  1  l)=n.o,XX(  15)=0. 1,XX( 16)=2  5.0,XX(  17)«25.P 
XX(  1R)  =  4.0; 

SEEPS ; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

I!!TLC,XX(  10)=0.0,XX(  1 1  )*0 .0 ,  XXII 5  )=P .  1 ,  XX  (  lf>)=2  5 .0,  XX(  17  )  =  22  .  S 
XXI  IB)-3.I>; 

SEEDS ; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 


SIMULATE; 

SIMULATE ; 

SIMULATE ; 

SIMULATE; 

INTLC ,  XX  (  10)=0.0,XX(  1!)=0.0,XX(  1  5)=0 .  1 ,  XX(  16 )=2 2 . 5 ,  XX(  17 )=2  S .0 f 
XX( 1R)=4 .0; 

SEEDS ; 

S  IMPUTE; 

SIMULATE; 

SIMULATE; 

SIMTLATE ; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC,XX( 10)=0.0,XX( 11)*0.0,XX( 15)=0.  1,XX( 16)=22. 5 ,XX( 17)=22 . 5 , 
XXT 18)*3 .6; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

* 

;  2ND  RLOCK  OF  64  RUNS  STARTS  HERE. 

INTLC ,XX( 10)=1.0,XX( 1 1)«1.0,XX( 15)=0.0,XXf 16)=2 S.O,XX( 17)=25.0, 
XX( IP )*4 . 0 ; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX( 10)«=3.0,XX(  1I)= 1.0, XX( I5)=0.0,XX(  16)=25.0,XX(  17)=22.5, 
XX(  lR)-3 .6 ; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 
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S  TMULATF. ; 

I  NTLC  ,  XX(  10)=].0,vv(  ]])-1.0,XX(  3  S  )=0 . 0 , XX(  If.  )=22 . 5  ,  V  v(  17)  =  2r..O 

xx c  ]«)=4.n; 

53  F.E  ns ; 

SIMULATE; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

SIMULATE ; 

INTLC.XXf  30)=1.0,XX(  1  1)*  1  .0  ,XX(  15)=0.O,XX(  If. )=22 .  S , XX(  17  )=2  2 . 5 
XX( 18)=3 .6; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC ,XX(  10)=  1.0, XX(  11)=).0,XX(  15)=0. 1,XX(  If>)*2  5.0,  XX(  27 )=2 5 . 0 
XX( lfi)=4 .0; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

SIMULATE; 

INTLC ,XX(  10)=1.0,XX( 1 1)=I.0,XX( 1S)=0. ],XX( 16)=25.0,XX( 17)=2  2.5 
XX( 18)=3 .6; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX(  10)=1.0,XX( 1 1)-].0,XX( 15)=0. 1,XX( I6)=22.5,XX(  17)=2  5.0 
XX(  18)  =  A.O; 

SEEDS; 

SIMULATE; 

SIMULATE; 
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SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC ,XX(  10)=1.0,XX(  ll)=1.0,XXf]5)=n. 1,XX( 1A)=22.S,XX( 17)=22.5 
XXf  1P)=3.6; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

! 

;  3RD  BLOCK  OF  64  RUNS  STARTS  HERE. 

INTLC, XX(  10)=  1.0, XXf H)=2 .0,XX(  1 5)=0 .0 ,XX(  16)=2  5.0,XX(  17)-2  5.0 
XX(  1R)=4 .0; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX( 10)=1.0,XX(  ll)  =  2.0,vv( ]5)=0.0,XX( J6)=25.0,XX(  17 )*2 2.5 
XX(  IP )=3 . 6 ; 

SEEDS ; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX( 10)=].P,XX( 1 1)=2.0,XX( 15)=0.0,XX(  1 6 )=2 2 . 5 , XX( 17 )«2 5 . 0 
XXf  1R)=4.0; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 


simulate ; 

SIMULATE ; 

SIMULATE; 

i rnrLC , xx(  io>=i.d,xx(  ii)=2.o,xx(  35)=n.n,vy(  if>)=22.5»y>:(  ]7)«22.5, 
XX(  18  5  =  3 .6; 

SEEDS; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

SIM I  LATE; 

SIMULATE; 

SIMULATE ; 

SIMULATE; 

SIMULATE; 

INTLC , XX( 10)-].0,XX(11)-2.0,XX( 15)-0. 1,XX(  16)*25.D,XX( ]7)-25.0, 
XX( 18)=4 .0; 

SEEDS; 

S IMULATE ; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE ' 

INTLC , XX(  10)=1.0,XX( ll)-2 .0,XX( 15)=0. 1,XX( 16)»25.0,XX(  17 >-22.5, 
XX (  135-3.6; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX(  10)=].o,XX( 11)=2.0,yv( i5)=o.  1,XX(  16)-22.5,XX(  37 )  =  2  5.0, 
XX(  lR)-4.0; 

SEEDS ; 

Sim.’ LATE ; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC, XX ( 10)-].0,XX( I1)=2  .0,XX(  15)-0 .  1 ,XX(  16)=22 . 5 ,XX(  17 )=2  2 . 5, 
XX(  13)— 3 .6; 

SEEDS; 
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SIMULATE 
SI  Ml' LATE 
SIMl'LATF. 
SIMILATE 
SIMULATE 
SIMULATF 
SIMULATE 
SIMULATF. 


4TH  RLOCK  OF  A4  RUN'S  STARTS  HERE. 


IUTLO , XX (  10)=I.0,XX(  U)=3.o,xy.(  15)*0.0,XX(  16)=25.0,XX(  17)-2  5.0 
XX(  1R)=4 .0; 

SEEDS; 

SIMULATE; 

SIMULATE 
SIMULATE 
SIMILATE 
SIMULATE 
SIMULATE 
SIMULATE 
SIMULATE; 

IN'TLC ,XX(  10)=1.0,XXni)  =  3.0,XXf  15)=0.0,XX(  16)=25.0,XX(  17)=22.5 
XXf 18 )=3 .6; 

SEEDS; 

SIMULATE; 

SIMULATE 
SIMULATE 
SIMILATE 
SIMULATE 
SIMULATE 
SIMULATE 
SIMULATE; 

INTLO ,XX( 10)-1.0,XX( 1  1)  =  3.0,XX( 1S)-A.0,XX( 16)=22.5,XX(  ]7)=2S.O 
XX( 18)«4.0; 

SEEDS; 

SIMULATE; 

SIMULATE 
SIMULATE 
SIMILATE 
SIMULATE 
SIMULATE 
SIMULATE 
SIMULATE; 

INTL  r,vy(  10)-].otxxn  J)-3.0,XX(  I  S  )-D .  n  t  XX(  I  b  )«22  .  S ,  X  XM  7  )  =  2  2 .  S 
XX(  18). 3.^; 

STEDS; 

S IMULATE ; 


SIMULATE 

SIMULATE 


SIMILATE ; 

SIMULATE ; 

SIMULATE ; 

SIMULATE ; 

SIMULATE; 

IN'TLC  ,XX(  m)=i.O,XX(  11)  =  1.0>XX(  ]S')=n.  1,XX(  16)=2  5.0  XX  (  ]  7  )  =  ■>  A  .0 

xxf  iR)=A.n;  ’ 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

IKTLC ,  XX(  10)=  ]  .0,XX(  1  1)=3 .0,  vy(  15)=o.  1,XX(  lf>)=25.0,XX(  1 7  )=2 2.5 
XX( IP )=3 • 6 ; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

INTLC,XX(  10)=1.^,XX(  1 1)=3 .0,XX(  15)=0.  1,XX(  16 )=2 2 . 3,XX(  1 7 )=*2  5.0 
XX(lfl)-4.0;  ’ 

SEEDS ; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

TNTLO ,XX( 10)=l.O,XX(ll)=3.0,XX( 15)-0. 1,XX( 1A)«22.S,XX( 17)-22.5 
XX( 1P)=3 .6; 

SEEDS; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMILATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

SIMULATE; 

FIN; 


* 

* 


* 

* 

* 

PR  OCR  AM  MAIN  ( I N  rJT ,  OUT  PUT ,  TA  PF.  5=1  N  PUT ,  TA  TE  f>=OUT  PUT ,  T  A  PF.  7 ,  T  A  PF,  ]  5 ) 
DIMENSION  NSET( lOOOO) 

COMMON/SCOM  1/ATRIP>(  100), DD(  100)  ,  DDL(  100)  ,DTNOW  ,  II , MFA, MSTOP, NCLNR 
+  , NC  RDR , N  PENT , NNRUN , UNSET , KTA  PE , S  S ( 1 00 ) , R  SL ( 100) ,TNEXT ,TN0V , 

+  XX(  100) 

COMMON  0SET( 10000) 

EQUIVALENCE  (NSET(  1)  ,OSF.T(  1)  ) 

NCRDR-5 

N PR NT-6 

NTAPE=7 

KNSF.T- 10000 

CALL  SLAM 

STOP 

END 

* 

* 

************************************************************************ 

* 

*  VALUES  USED  THROUGH  OUT  THE  PROGRAM  ARE  INITIALIZED  AND  READ  INTO 

*  FILE  ONE, WHERE  THE  RANKING  IS  BASED  ON  THE  LOW  VALUE 

*  OF  ATRIBUTE  I, THE  UNIT  DESIGNATION. 

* 

************************************************************************ 

* 

* 

SUBROUTINE  INTLC 

COMMON/SCOM  1/ATRIB( 100) , DD( 100 ) , DDL( 100 ) , DTNOW , 1 1 ,MFA, MS TO  P, NCLNR 
+  , NCR DR , N  PR NT , NNRUN , NNS ET,NTAPE,SS( 100) , SSL ( 100) ,TNEXT, TNOW , 

+  XX( 100) 

COMMPN/UCOM  1/TIME, TIME  1  ,TIME2,TI*1E3,TIME4, TIMES,  POLICY, ROUTE 
+  ,12,13,  J  1,J2,  J3,DIF  1,DIF2,DIF3,DTIME,DTIMF.P,NUMKIL,DTMIS  1, 

+  DTMIS2 .DTMIS3 , I , J 
SAVE/UCOM 1/ 

DIMENSION  A(20),B(20),C(20) 

REAL  TIME .TIME  1  ,TIME2 ,TIME3 .TIME4 ,TIME5 , DIF  1 , DIF 2 , DIF 3 , 

+  DT  IMF. ,  DT IM  E  T ,  DTM I S  1 ,  DTM I S  2  ,  DTMI S  3 ,  T IMEM 

INTEGER  I, I 1,I2,I3,J,J1,J2,J3,K,L,M,N, POLICY, ROUTE, 

+  KUMKI L ,  NU  PLAN ,  MS110T 

DATA  A, B,C/20*0. 0,20*0. 0,20*0.0/ 

DIF  1-0.0 
DIF2-0.0 
DIF3-0.0 
XX( 1)-1 
XX(2)-16 
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.sV- ‘  V 


XX ( 3  )  = 
XX(4)= 
>:x(5)= 
xx(f.)= 

XX(  7 )  = 
XX(9)= 
XX(  12) 
XX(  13) 
XX(  14) 
XX( 19) 
XX(20) 
XX(2  1) 
XX(22) 
XX(23) 
XX(24) 
XX(25) 
XX(26) 
XX( 27) 
XX(  28  ); 


32 
0.0 
0.0 
O.o 
30.0 
‘  100.0 
=  0.0 
1=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 
=0.0 


t*********************************************************************** 

k  ROUTE  BLUE  (ENTITIES  1-15).  ATRIBUTE5  1-12  FOR  ENTITIES  1-15  ARE 

k  READ  INTO  FILE  ONE.  ATRIBUTE  ONE  IS  THE  UNIT  NUMBER.  OTHER 

k  A  TRIBUTES  ARF,  EXPLAINED  IN  THE  SLAM  CODED  PORTION  0^  THE  MODEL  ABOY 

* 

************************************************************************ 


A(  1)=1 
A(2  )=  1 .0 
A(3)=57 
A(4 )=32 
A( 5 )=2 . 4 
A( 6 )=4  .  B 

A(  7  )=3 .  A-XX( 15)*3 .6 
A(8)=TRIAC(A(5),A(7),A(6),4) 

A(  1  l)=TRIA0(20.ntxX( 16), 30.0, 3) 
A(  12)=TRIAG(20.0,XX(  17), 30. 0,3) 
A(9)=64.784/A(  11) 

A(  10)=14.864/A(  11) 

CALL  BILFM( 1 , A) 

A(  1  )=2 
A( 2 )=  1 .0 
A( 3 )*5 
A( 4 )  =  2  5 
A( 5 )=0 .75 
A(6)=  1.5 

A(7)= 1.  125-XXf  15)*1.  125 
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M P  )-TRI .V* ( A< 5 ) , Af  7  > , A( b ) , 4 ) 

A(  !  l)=TP.IAC(20.ritvy(  ]ft)  ,30.0,  ]) 
AC  12)=TfvIAC(3.o,y>;(  IP)  ,5.0,3) 
A(0)=»64.7P,4/A(  1 1) 

A( in)»14.RA4/A( 11) 

CALL  F I I.F.M C  1 ,  A  ) 

A(  1)  =  3 

A(2)*J.O 

A(3)=2» 

A( 4 )*2 
A( 5 )=0 . 9 
A(6)»1.R 

A(7)*1.35-XX( IS)*!. 35 
A<<3)=TF.IAG(A(5),A(7),A(A),4) 

A( 1 1)=TRIAG(20.0,XX( 16), 30.0, ]) 
AC  12)*TRIAC(3.0,XXC IS), 5. 0,3) 
A(9)*64.784/A(  1 1) 

A(  10)«14.R6A/A( 11) 

CALL  F1LEMC 1 , A) 

AC  1  )*4 
AC2)*3.0 
AC  3 )*33 
A<4)=9 
AC  5 )=  1 .  1 
A(  A )=2 .2 

AC  7  )■=  1 .  A  5-XXC  15)*1.65 
A(R)=TRIACCA(5),A(7),ACA),4) 

A( 11)=TRIACC20.0,XXC 16), 30.0, 1) 
AC  12)-TRIACC3.0,XXC 18), 5. 0,3) 
AC9)=R4.7R4/AC 11) 

AC  10)*  14 .  R64  /AC  11) 

CALL  FI  LEM ( 1,A) 

A( 1)*5 
AC2)=1.0 
AC  3 )=30 
A( 4 )=9 
A(5)-).0 
AC  6)=2 . 0 

AC  7 )= 1 . 5-XXC  15)*  1 . 5 
A(R)*TRIAC(A(5) ,A(7) ,  A(6)  ,4) 

A(  11)-TR1ACC20.0,XXC 16), 30.0,  ]) 
AC  12)-TRTACC3.0fXXC  18), 5. 0,3) 
A(9)=64.784/A( 1 1) 

AC  I0)*14. 864/aC 11) 

CALL  FTLEMC 1 , A) 

AC 1)*6 
AC  2 )■  1 .0 
AC  3 )*7 

AC4)-R9 

AC5)-2.5 
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AC  A  )=5 .0 

Af7)=3.7S-XX( 35)*3.75 
A(8)=TRIAC(A( 5 ) ,AC  7  ) ,  A(  A  )  ,  4  ) 

AC  ll)=TRIAC(20.0,yv( 16),30.n  j> 
A(  12)=TRTAC(3.0,XX( 1ft), 5. 0,3) 
AC9)=fS4.784/A(  11) 

A( 10)=  14.P64/A(  11) 

OALL  FILEM( 1 , A) 

A(  1  )=7 
A( 2)»  1 .0 
A( 3  )  =  5  1 
A(  4 )=2 
A( R )=  1 ,4 
A(  6)=2 .8 

A( 7 )=2 .  1-vy(  15)*2.  1 

A( ft )=TRI AG( A( 5 ) , A( 7 ) , A( 6 ) , 4  ) 

Ar  1  D“TRlAO(20.0,yx(  1  A),  30.0,  1) 
A(  12 )=TRIACC3 . 0, XX(  1ft), 5. 0,3) 
AC°)*A4.7R4/AC 11) 

A(  10)=14.R64/A(  13) 

CALL  FILEM(l.A) 

AC  l)-8 
A(2)=1.0 
A( 3  )  =  5 
A(4)=30 
A( 5 )*0 .875 

A( A )=  ] . 7  5 

A( 7 )* 1 . 3  125-XXC 15)* 1 . 312ft 
AC8)*TRIACCAC5),A(7),ACA),4) 

A( 1 1)-TRIAC(20.0,XX( 36),3O.0, 3) 

A( 1 2 )=TRT ACC  3 .0 , vv( 3ft)  5,d  3) 

A(9)»64.7«4/A(  11) 

AC  10)= ]4.ftA4/A(  1 1) 

OALL  FI LF.MC  1 , A) 

AC  l)=9 
A(2)-  1.0 
AC3)=0 

A(  4  )=  15 
AC  5)=0 .4 
A( 6 )*0 . ft 

A(7)=0.6-XX( 15)*0,f> 

A( 8 )*TRI ACC A( 5) , A( 7 ) , A( 6) ,4 ) 

AC  1 1)«TRIACC20.0,XXC 1 A), 30.0,  1) 

A( 12)=TRIAG(3.0,XX( 1ft), 5. 0,3) 
AC9)»A4.7R4/AC  31) 

AC  10)= 14.R64/AC  11) 

CALI  FILFM(l.A) 

AC  1)-10 
AC  2 )»  1 .0 
A(3)=0 
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A  f'O  =  50 
A(  5  )  =  1 . 6 
A(6>)«3.2 

A(7)=2.4-XX( 1S)*2.4 
A(8)=TRIAC(A(5),A(7),A(6),4) 

A(  ll)-TRIAC(2n.n,xx(  16), 30.0,  1) 
AC  12)=TRIAC(3.0,XX(  18), 5. 0,3) 
A(9)=64.784/A<  11) 

A(  10)=  14.864 /A( 11) 

CALL  FI  LEM (1, A) 

A( 1 )=  1  1 

A(2)-I.O 

A(3)=0 

A(4)=59 

A(5)=1.6 

A( 6)  =  3 .2 

A(7)=2.4-XX( 3  5 ) *2 . 4 
A(8)=TRIAC(A(5),A(7),A(6),4) 

A(  11)=TRIAC(20.0,XX(  16), 30.0,  1) 
A(  12)«TFIAC(3.0,XX( 18), 5. 0,3) 
A(9)=64.784/A(  11) 

A(  10)= 14. B64/AC  11) 

GALL  FI  LEM (  1,A) 

A( 1)=  12 
A( 2 )=  1 .0 
A( 3 )=0 
A( 4 )=  18 
A(5)«=0.5 
A(  6 )=  1 .0 

A(7)=0.75-XX(  15) *0 .75 
A(8)=TRIAC(A(5),A(7),A(6),4) 

A( 1  l)*TRIAC(20.0,vx( 16), 30.0, 1) 
A( 12)«TRTAC(3.0,XX(  18), 5. 0,3) 
A(9)=64.784/A(  1 1) 

A(  10)= 14.864/A(  11) 

CALL  FILEM(l.A) 

A(  1  )=  13 
A(2)-1.0 
A(3)-0 
A( 4 )=7  1 
A(  5  )=  1 . 9 
A(6)=3 .8 

A(7)=2.85-XX(  15) *2 .85 
A(8)=TRIAC(A(5),A(7),A(6),4) 

A( 1  1)*TRIAC(20.0,XX(  16), 30.0, 1) 
A(  12)=TRIA0(3.0,XX( 18), 5. 0,3) 

A( 9)=64 . 784/A(  11) 

A(  10)=  14 . 864  /  A(  1 1) 

CALL  FI  LEM ( 1,A) 

A(  1)=  14 
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r/.  -• 

t- n  ; 


* 

★ 


4(2)=  1.0 
a(3)=o 
AC  4  )=4  R 
A( 5)=  1 . 3 
AC  A  )=2 .6 

A(7)=  1.Q5-XXC 15)* ] .05 
A(P)=TRIAC.(AC5),A(7),A(ft),4) 

A( 1  1)=TR1AG(20.0,vv( lfi) ,30.0, 3) 
A(  12  )=TRI ACC  3 .0 , XX( IP), 5. 0,3) 

AC  0 )=64 . 7P4 /A(  11) 

AC  10)= 14.R64/A(  11) 

CALL  FI  LEM (  1,A) 

A(  1)-15 
A(  2 )=  1.0 
A( 3)=0 
A( 4 )=202 
A(5)=5.  1 
A(  6 )=  10. 2 

A( 7 )=7 . ft5-XX( 1 5)*7 .65 
A(B)=TRIACCA(5),A(7),A(6),4) 

AC  1 1)=TRIAC(20.0,XX( lft) ,30.0, 1) 
A(  12)=TRIACC  3.0.XXC IP), 5. 0,3) 
A(9)*64.7B4/A( 11) 

A(  10)= 14.864 /AC  1 1) 

CALL  FILEM(  1,A) 


********************************************************************** 

* 

*  ROUTE  GREEN  (ENTITIES  16-31).  ATRIRUTES  1-12  FOR  ENTITIES  16-3  1  ARE 

*  READ  INTO  FILE  ONE.  ATRIRUTE  ONE  THE  UNIT  NUMBER.  OTHER 

*  ATRIBUTES  ARE  EXPLAINED  IN  THE  SLAM  CODED  PORTION  OF  THE  MODEL  ABOVE 

* 

********************************************************************** 

* 


B( 1)=  16 
B( 2 )=  1 .0 
B(3)»57 
B( 4 )=32 
B( 5)*2 .4 
B(ft)=4.R 

F(7)=3.ft-XXC  15)*3.ft 

B( R )=TRIAG(B( 5) ,B(7),B(6) ,4) 

B(  11)-TRIAC(20.0,XX(  16), 30.0,  1) 
P,C  12)=>TRIACC20.0,XX(  17), 30. 0,3) 
P(9)=64.972/B(  11) 

PC  10)=7.62R/B(  11) 

CALL  FILFMC 1,B) 

R(  1)=17 


m 

Sv 


RC  2  )=  1 . 0 
PC  3)  =  5 
R(4)=25 
P(  5  )*0 . 7  5 
B(6)=1.5 

B(7)-l.  125-XXC 15)*J.  125 
B(  °  )-TRI  ACC  P(  5  )  ,  B(  7  )  ,  P,(  6)  ,4  ) 

B( n)=TRIAC(20.n,yv( 16), 30.0, 1) 

B(  12)*TRIAG(3.0,XX(18),5.0,3) 

PC  9  )=*64 . 972 /R  (  1 1) 

B(  30)=7.628/PC  1 1) 

CALL  F I  LEM (  1 ,  P. ) 

R(  1 )=  18 
P(  2 )*  1.0 
B( 3)  =  28 
R(  4  )*=2 
R( 5 )=0 .9 
R(  6)*=  1 .8 

R( 7 )=  1 . 35-XXC  15)*]. 3S 
PC8)-TRIAC(P.(5),RC7),B(6),4) 

PC  1 1)*TRI AO( 20.0,XXC 16), 30.0, 1) 
RC  12)=TRIACC3.0,XXC  18), 5. 0,3) 
RC9)=64.972/RC  11) 

BC  10)*7 .628/BC  11) 

CALL  FILEMC 1 , B) 

RC  1>= IP 

BC  2 )=  1 .0 
BC3)«33 
RC4)-9 
R(  5 )■  1 . 1 
RC  6 )=2 .2 

RC  7  )=  1.65-XX( 15)* 1 .65 
B(8)-TRIAGCPC5),B(7),BC6),4) 

PC  1  1)*TRIACC20.0,XXC 16), 30.0, ]) 
PC  12)=TR1ACC3.0,XX( 18), 5. 0,3) 
P.C9)=64.972/RC  11) 

BC  10)*7. 628/BC 11) 

CALL  FILEMC 1,R) 

RC  l)-20 
EC  2 )■  1.0 
RC  3 )*30 
BC4)-9 
BC  5)»  1 .0 
RC6)-2.0 

PC  7  )*=  1. 5-XXC  15)*  1 . 5 
BC*)=TRIACCBC5),BC7),BC6),4) 

RC  11)-TRIACC20.0,XXC 16), 30.0, 1) 
BC  12)-TRIACC3.0,XXC 18), 5. 0,3) 
BC9)-64.972/RC  11) 

RC  10)=7. 628/BC  11) 
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CALL  FI  LEM (  1,P) 

B(  ])-2  1 

R(2)=1.0 

B( 3 )=7 
P(4)=89 
B( 5 )=2 . 5 
B(6)=5.0 

B( 7 )  =  3 . 75-XX( 1 5 ) *3 . 75 
B(8)*TRIAC(B(5),B(7),R(6),4) 

B(  11)=TRIAC(20.0,XX(  16), 30. Cl,  I) 
B(  12 )=TRI AC( 3 . 0 , XX( IB), 5. 0,3) 

B( 9  )=64 . 972 /B( 1  1) 

R(  10)=7 . 628/R(  11) 

CALL  FILF.M(  1,B) 

P.(  1  )=22 
R(  2 )®  1 .0 
B(  3 )=0 
B(4)=59 
B( 5)-  1.6 
B(6)*3.2 

B(7)=2.4-XX( 15)*2.4 
B(8)*=TRIA0(F(5),P(7),B(6),4) 

B(  11)=TRIAC(?0.0,XX( 16), 3O.0, 3) 
B(  12  )*TP.IAG(  3 .0,  XX(  18),  5. 0,3) 
B(9)=64.972/B( 1  1) 

B(  10)=»7.62P/B(  11) 

CALL  FILEM(l.B) 

B(  1)=23 
B(2)-1.0 
P( 3 )=0 
B( 4 )=57 
R( 5 )=  1 . 5 
B( 6 )  =  3 .0 

B( 7 )  =  2 .25-XXC 15) *2 .25 

B(  P )kTRTA0( B( 5 ) , B( 7 ) , B( 6 ) , 4  ) 

B( 11)-TRIAC(20.0,XX(  16), 30.0,  1) 
P(  12)=TRIA0(3.0,XX(  18), 5. 0,3) 
B(9)-64.972/B(  11) 

B( 10)*7.628/B( 11) 

CALL  FI  LEM ( 1,B) 

B(  l)-24 
B(2)-1.0 
B(3)-0 
B(  4  )=56 
B( 5 )=  1 . 5 
R(6)*3 .0 

B(7)=2.25-XX(  15) *2 .25 
B(8)«TRIA0(B(5),B(7),R(6),4) 

F( 11)“TRIAO(20.0,XX(  16), 30.0,  1) 
B(  12)*TRI  AC.(  3 .0,  XX(  18),  5. 0,3) 
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prp)=<si.r'72/r’(  i i) 

F.r  io)=7.62P/r,(  id 
CALL  FT  LEV'(  1,B) 

Ti(  1  )  =  2  5 
Br2)«1.0 
p.n>=2  2 

R(4)=76 
P( 5)=2 . 5 
P(6)*5.n 

B(7)  =  3.75-XYC  1 S ) *3  -  75 
P(8)-TLIAC(R(5),RC7),R(6),4) 

R(  n)-Tt'IAC(20.0,XX(  16),3n.O,  ]) 
B(  12)=TP.AG(3.0,XX( 1R) , 5.0,3) 
B(9)«64.972/B(  11) 

P(  10)=7 .62R/PC 1 1) 

CALL  FILF.M(1,R) 

B(  1)=26 
B( 2  )=  1 .0 
B(3)-0 
P(4)»56 
B(  5 )*»  1 . 5 
R(6)*3.0 

B( 7 )=2 . 2  5-XX(  15)*2 . 2  5 
B(R)=TRIAC(B(5),B(7),B(6),4) 

B( ll)-TRIAC(20.0,XXf 16), 30.0, ]) 
B( 12)-TRIAC(3.0,XX( 18), 5. 0,3) 
R(9)~64.972/P.(  11) 

P(  10)«7.628/B( 11) 

CALL  F I LEM (  1,B) 

B(  1  )=27 
B(  2 )=  1.0 
B(3)=0 
R( 4 )=66 
B(  5  )=  1 . 8 
B(6)=3.6 

RC7)=2.7-XX(15)*2.7 

B(8)=TRIAC(P(5),R(7),R(6),4) 

B( 11)-TRIAC(20.0,XX( 16), 30.0, 1) 
R(  12)-TRIAG(3.0,XXf 18) ,5.0,3) 
R(9)«64.972/B(  11) 

R(  10)-7.628/B( 11) 

CALL  FI  LEM ( 1,B) 

R(  1)“2R 
B(2 )*  1.0 
B(3)-0 
R(4)«40 
B( 5)*1 . 1 
P(6)-2.2 

R(7)»1.65-XX(  15)*  1 .65 
P(8)«TRIAG(B(5),B(7),B(6),4) 


Bf  1  ])=TRTAr(20.n>Yx(  if.), 30.0,  i) 

?•(  12)=TPIACf  3 .0,  yy.f  1 8 ) ,  5 . n ,  3 ) 

P(9)=64.Q72/Rf  1 1) 

Rf  10)=7.628/Bf 11) 

('ALL  FILEMf  1,B) 

Rf  1)=29 
Rf 2)=  1.0 
P(3)=5 
Bf  4  )=60 
R(  5 )=  1 . 7 
B(  6  )=3 . 4 

B( 7 )=2 . 55-XX( 15)*2 . 55 

Bf  8  )*=TRIAGf  R(5),Bf7),Bf6),4) 

P.(  ll)*TRIAGf20.0,XXf  16)  ,30.0,  1) 

R(  12)=TRIACf3.0,XXf  18) ,5.0,3) 

B(9)=64.972/Rf  11) 

B(  10)=7 . 62R/B( 11) 

CALL  FILEMf  1,P.) 

B(  1)=30 
B( 2)=  1 .0 
Rf3)=0 
R(  4 )=  18 
Bf 5 )=0 . 5 
Rf  6)*  1 .0 

R( 7 )=0 . 75-XXf 15) *0 .75 
R(8)-TRIAC(B(5),B(7),B(6),4) 

Rf  1  l)-TRIACf 20.0 , XXf  16), 30.0,  1) 

B(  12)=TRIACf3.0,XXf 18), 5. 0,3) 

B(9)*64.972/B( 11) 

R(  10)«7.628/Rf 11) 

CALL  FILEMf  1, P.) 

R( 1 )=3  1 
Rf 2)=  1.0 
B(3)«0 
Rf  4 )=3  1 
Pf5)=0.9 
Rf  6)  =  1 . 8 

Rf  7 )■ 1 . 35-XXf  15)*  1 . 35 

Rffl )=TRI ACf Bf  5) ,B(7 )  ,R(  6)  ,4) 

Rf 1 l)«TRIAGf20.0,XXf 16), 30.0,  1) 

Rf  12)-TRIACf3.0,XXf  18), 5. 0,3) 

Rf9)-64.972/Rf  11) 

Rf  10)»7 .628/Bf 1 1) 

CALL  FILEMf 1,R) 

* 

* 

******************************************** ********************* 

* 

*  ROUTE  RED  fF.NTITIES  32-44).  ATRIBUTES  1-12  FOR  ENTITIES  32-44  ARE 

*  REAn  INTO  PILE  ONE.  ATPIBUTE  ONE  IS  THE  UNIT  NUMBER.  OTHER 
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“DPF.L  AR^yr 


*  ATRIP. IT rr>  ARE  EX TLAI  NTEr>  IN'  TrC  ELAM  CODED  PRIORI I oy  OF  TDD 

★ 

*  ****************************************************************** 


c(  1)  =  37 
C(2  )=  1 .0 
C(3)=5 
0(4)= 30 
C( 5)=0.875 
0(6)=1.75 

C( 7 )= 1 . 3 125-XX(  15)*1. 3125 
C(8)=TRIAG(C(5),C(7),C(6),4) 

C(  11)*TRIAC(20.0,XX( 16), 30.0, 1) 
0(  12 )=TRIAG( 3 . 0, XX( 18) ,5.0,3) 
0(9 )=56.628/0( 11) 

0(  10)=  11.6 16/0( 11) 

CALL  FILEM( 1,C) 

C(  1)=32 
C( 2 )=  1 .0 
C(  3 )=39 
0( 4 )=8 
0(  5 )=  1 . 5 
C(6)=3 .0 

C(7)=2.25-XX( 15) *2 .25 

0(8 )=TRIAC(C( 5),C(7),C(6),4) 

C( 1 1)=TRIAG(20.0,XX(  16), 30.0,  1) 
C(  12)=TRIAC(20.0,XX(  17) ,30.0,3) 
G(9)=56.628/C(  11) 

0(  10)=  1 1.6 16/C(  11) 

CALL  FILEM(  1,C) 

C(  1)=33 
C( 2 )=  1 .0 
0(3 )=6 
0(4 )=  10 
C( 5)=0 . 5 
0(6)=  1.0 

0(7 )*0.75-XX(  1 5 ) *0 . 7  5 
0( 8)«TRTAO( C( 5), 0(7), 0(6), 4) 

C(  1 1)=TRIAG( 20.0,XX(  lrb 
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•  .N  .*• 


c(  i?)=triao(3.o,xx(  l°),s.n,3) 

r<,9)=56.62P/r(  n) 

r.(  10)=  1  ].ftlft/C(  1  ]) 

tall  filem(i,o 

C(  1)-3A 

C( 2 )=  1 .0 

C( 3 )=33 

C(4)=0 

C(5)»J.O 

C(6)=2 .0 

C(7)-1.5-XX(  15)*  1 . 5 
C(R)»TRIAC(C(5),C(7),C(6),4) 

C(  1  1)=TRIAC(20.0,XX( 1ft) ,30.0, 1) 
re  12)=TRIAC(3.0,XX(  18) ,5.0,3) 
C(9)=56.ft2£/C( 1 1) 

C(  10)= 11.6 16/C( 1 1) 

CALL  FI LEM ( 1,C) 

C(  1  )=3  5 
C(  2  )=  1 .0 
C(3)=23 
C(4)-0 
C( 5)=0 . 7 
C(ft)=  1.7 

C( 7 )=  1 . 2-XX(  1 5 )*  1 .2 

C(  8 )=TRIAC(0(  5),C(7),C(ft),4) 

C(  11)=TRIAC(20.0,XX(  16),30.0,  1) 
C(  12)=TRIAC(3 .0 ,XX( 18), 5. 0,3) 
C(9)=56.628/C( 11) 

C(  10)= 11.6 lft/C( 11) 

CALL  FI LEM (1,C) 

C(  1)=36 
C( 2 )=  1 .0 
C( 3)=  15 
C( 4 )=89 
C(5)=2.ft 
C( 6)*5 .2 

C( 7 )=3 . 9-XX(  15)*3 .9 
C(8)=TRIAC(C(5),C(7),C(ft),4) 

C(  11)*TRIAC(20.0,XX(  16), 30.0, 1) 
C(  12)=TRIAC(3 .0 ,XX( 18), 5. 0,3) 
C(9)-5ft.ft2R/C( 1 1) 

C( 10)=ll.ftlft/C( 11) 

CALL  FI LEM (1,C) 

C(  1  )=3P 
C(2)=1.0 
C(3)=0 
C( 4 )=  IP 
C( 5 )=0  .  ft 
C(ft)=1.0 

C(7)=0.75-XX(  1 5 ) *0 . 7  5 


r(F)*TPIAC(r(!>),cm,C('.),4) 

r.(  m=TP.i  v;(2^.o,yy(  i6),3o.nt  n 
0  (  1 2  )=T  FIT  AO  ( 3 .  ° ,  X  V  (  1 q  ) ,  *■ .  0 , 3 ) 
rf  o  )  =  56- 62<--/Cf  11) 
r.(  lo)=i  i.a ia/c(  i  ]) 

CALI.  FI  LE“(  1,C) 

C(  1)=3« 

C( 2 )=  1 .0 
C(  3 )  =  5  7 
C(4  )=32 
C( 5 )=2 . 4 
C( 6 )=4  .  P 

C(7)=3.6-XXf  15)*3.6 
C(8)*TRIACfC(5),C(7),C(6),4) 

C(  11)=TP.IAC(20.0,XX(  16), 30.0,  1) 
C(  12)«TRTAcn.0,xx(  18),  5. 0,3) 
C(9)=56.62F/C(  11) 

C(  10)= 11,6 16/C(  11) 

CALL  FILFMf 1,C) 

C(  1)=40 
C(2)=  1.0 
0O)= 5 
C( 4  )=2  5 
C( 5 )=0 .75 
C(6  )=  1 . 5 

C( 7 )= 1 .  125-XXf 1 S)*  1 .  125 
C(R)=TPIAC(C(5),C(7),C(6),4) 
cm)=TRIA0(20.0,XX(  16), 30.0,  1) 
C(  12)=TRTAC(3.0,XX( 18), 5. 0,3) 
C(9)=56.628/C( 11) 

C( 10)=11.616/C( 11) 

CALL  FILEM( 1,C) 

C( 1)=4  1 
C( 2 )=  1 .0 
C(3)=2R 
C( 4 )=2 
r(3)=o.9 
C(6 )=  1 .8 

C(7)=1.35-XX( 1 S)* 1 . 35 
C(8)=TRIAC(C(S),C(7),C(6),4) 

C( 1 1 )=TR1 AC( 20.0 , XX(  16), 30.0,  I) 
C( 12)=TRIAC(3.0,XX( 18), 5. 0,3) 
C(9)=56.628/C( 11) 

C( 10)= 1 1.6  !6/C( 1 1) 

CALL  FILF.M(l.C) 

C(  l)-42 
C( 2 )=  1 ,0 
C( 3 )=33 
C( 4 )*9 
C(S)=l.l 


lil 


r. 6)=2.? 

C(7)=  1.65-XX(  15)*1.A5 
C(P)=TRIAO(C(5),0(7),C(A),4) 

C( ll)=TRIACr20.0,XX( 16), 30.0, ]) 
C(  12)=TRIAG(3.0,XX(  ip.)  ,5.0,3) 
C(9)=5A.A2R/C(  1 1) 

c(  io)*ii.ai f>/c.(  ii) 

GALL  FILEM( 1,G) 


C(  l)-43 
C( 2 )=  1 .0 
G( 3 )=23 
C(4)=9 
C(  5 )=  1 .  P 
C(  A  )=  1 . 6 

C( 7 )=  1 . 2-XX(  15)* 1 . 2 

G(R  )=TRIAG(C( 5) ,C(7)  ,C(6) ,4) 

C( 11)-TRIAG(20.0,XX( 16), 30.0,  1) 
C(  12  )=TRIAC( 3 . 0 , XX(  18)  ,5.0,3) 
C(9)=56.62P/C<  11) 

C(  10)=  11.6 16/G( 11) 

CALL  FILEM( 1,G) 

C( 1 )=4 A 
C(2 )=  1 .0 
C(  3 )*  14 
C( 4 )=  124 
C(5)=3.7 
C( A )=7 .4 

C(  7 )=5 . 55-XX(  1 5 ) *  5 . 5  5 

C( 8 )*TRIA0(C( 5 ) ,C(7 ) ,C( A ) ,4 ) 

C(  1 1)=TRIAG(20.0,XX( 1A) ,30.0, 1) 
C(  12)=TR1AC(3.0,XX( IP), 5. 0,3) 
C(0)=56.A2P/C( 11) 

C(  10)=11.A16/C(  11) 

CALL  FILEM( 1,C) 


****************************************************************** 

* 

*  NETWORK  PASSAGE  IS  INITIATED  BY  CALLING  THE  EVENT  NODE 

*  AT  THE  START  OF  EACH  ROAR  NETWORK. 

* 

****************************************************************** 


CALL  SCEDL(  1,0.0, ATRIB) 
CALL  SCHDLC20,0.0,ATRIR) 
CALL  RCHDL( 32 ,0.0, ATRIP) 
RETURN 
END 


.s_.s  >Vv 


subroutine  evf.nt(  in) 

common/ech" 1/ atrip (  inn) .one  ioo),ddl(  ioo) ,dtnow,ii  ,mfa,mstdp,nclnr 
+  ,NCRPR,Nm:T,NMEVN,NNSET,NTAPE,ss(  ioo) ,ssi.(  inn) ,tnext,  tngv, 

+  XX(  100) 

C0MMON/UC0M1 /TIME, TIME  3 ,TIME2 , TIME  3 , TIME4 .TIME  5 ,  POLICY, ROUTE 
+  ,  T2,T3,  J  1,  J2,  J3.DIF  1 , DIF2 , DIF3, DTIME , PTIME P.NUMKIL, nP'IE  1 , 

+  P?MIS2,DTMIS3,I, J 
!  VE/UCOM  1  / 

R: AL  TIME, TIME  1  ,TIME2 ,TIME3,TIME4 , TIME 5, DIF  1,DIF2,DIF3, 

+  DT IME , DTTMFF, DTMTS 1 , DTMIS2 , PTMIS3 ,TIMEM 

INTEGER  I , I  1 , 1 2 , 1 3 , J , J 1 , J2 , J3 , F , L , M, N ,  mLICY, ROUTE, 

+  NUMKIL, NUFLAN.MSHOT 
POLICY=XX( 10) 

CO  T0( 100,200,300,400,500,600,700,800,900, 1000, 3100,1200,  3300, 

+  1400,  1500,  1600,  1700,  1800,  1900,2000,2100,2200,2300,2400, 

+  2 500, 2600, 2700, 2800, 2000, 3000, 3  100, 3200, 3300, 3400, 3500, 36no, 3700, 

+  3800, 3900, 4000), IN 

* 

* 

********************************************************************* 

* 

*  ROUTE  BLUE  (ENTITIES  1-15).  FIRST  FAT. NT  MOOR  STARTS  THE  SCHEDULING 

*  OF  THE  CONVOYS .  THE  HEAD  &  TAIL  OF  THE  CONVOY  ARE  SCHEDULED  TO  COME 

*  INTO  THF.  NETWORF.  WITH  THE  ENTER  MODES  2  &  3  AT  THE  TIME  THAT  ALLOWS  FOR 

*  PASSAGE  OF  THE  VEHICLES  IN  RETWEEN  THE  HEAD  AND  TAIL.  SUBSEQUENT  EVENTS 

*  RECOMPUTE  ATRIRUTE  10, THE  TIME  THROUGH  THE  NEXT  NODES. 

*  END  TIME  TS  COMPUTED  AND  SENT  TO  TATE  IN  EVENT  19. 

*  WHEN  ATRIRUTE  2  VALUE  IS  2, THIS  REPRESENTS  THE  TAIL  OF  THE 

*  CONVOY.  WHEN  THE  VALUE  OF  ATRIRUTE  2  IS  1(HEAD),IT  BYPASSES 

*  THE  LAST  EVENT  IN  THE  APPROPRIATE  ROAD  SEGMENT. 

* 

********************************************************************* 

* 

* 

100  TF(XX(  1) .LE.  15.0)THEN 
1 1  =  INT(XX(  1)) 

CALL  COPY(I  1,  l,ATRIRi) 

* 

*  CONVOY  RATE  IS  ADJUSTED. 

* 

CALL  RESHUF1 

IF(XX( 10). EO. 1.0)THEN 

ATRIB(9)=50.0/ATRIR( 13) 

ENDIF 

TF( I  1  .EO. DTHEM 
T IME -T NOW 
ATRIR( 1 3 )=TNOW 


I  ~  ' 


•  -r 


xf, 


*  schedule  ufad  °f  convoy  i 


CALL  E CURL (2  ,TIME, ATRIP) 
ENDIF 

ATRIR( 2 )=2 .0 


*  REINITIALISE  ATTRIBUTE(R)  DUE  TO  RATE  CHANCE 

* 

if(xx(  ]0).r.o.n.o)THnr: 

ATR£S(<l)=64.7R4/ATRIR(  11) 

ENDIF 

ctnti=atrir(  12) 

ATRIM  13)=T.N0V 
TIME=ATRTB(8)/ATRIP(  11) 


*  SCHEDULE  TAIL  ENTITY 


CALL  SCHDL(3 ,TIME, ATRIB) 
IF(I  l.EO. 15)00  TO  150 
1 1=1 1+1 

CALL  COTYd  1,  l.ATRIB) 
CALL  RESULT  1 


*  RATE  ADJUSTMENT 


IF( XX(  10). EO. 1.0)THEN 

ATRIB(9 )=50.0/ATRIB( 11) 
ENDIF 

TIME=(TIME+(CINT1/ATRIR( 1 1))) 
IF(XX(  10).E0.0.n)THEN 

ATRIR(  9  )=f>4 .  784  /  ATRIP  (  1  1) 
ENDIF 


*  SCHEDULE  HEAD  ENTITY  OF  NEXT  CONVOY 


* 

CALL  SCHDL(2, TIME, ATRIP.) 

ft 

ATRIR(  13)=TIMF. 

£ 

*  RESCHEDULE  EVENT  ft  1  TO  CONTIN 

1 » 

* 

i  k  * 

■ 

CALL  SCHDL( 1 ,TIME,ATRIB) 

id 

150 

XX(  1)=YY( 1)+1 

ENDIF 

.* 

RETURN 

* 

200 

CALL  F.NTER(2,ATRIB) 

y!; 

RETURN 

■y. 

300 

CALL  ENTER (3, ATRIP) 

>• 

RETURN 

r.j 

BOO 

CALL  ENTER (B, ATRIP.) 
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•  •  •  • 


v,..*.. 


RETURN 

1  100  IF!(XX(  ln)  .E^.O.rv)  .A\':i.(ATRir.(  n.CE.  l.o))THE'* 
ATRIP!  ]0)=  1 4 .  R  A4  /AT" I  P(  ]  l) 

end  if 

IF(!XX(  10). F,o.  1.0). AND.  (ATRIP-!  l).F.O.  1.0))THEN 
ATRIP!  10)  =  1 4 . 7 B4 /ATRI R(  11) 

F.  nr  T  F 

IF!  (XX(  10)  .F.o.  1.0)  .AND.(  ATRIP !  1)  .OT.  1.0))THEN 
ATRI R(  10)=2Q.Q0/aTF.TP,(  11) 

FNDIF 

RETURN 

* 

*  INTERDICTION  OF  CONV OY  #1,!!F.AP  ELEMENT 

* 

3200  IF! (ATRIR( l).EO. 1 .0) .AND. ( ATRIP! 2 ) .EO . 1.0). AND. 

+  (XX(  10)  .  F,0.  1).AN0.(XX(  12)  .EO.O.O)  )TI!EN 

CALI.  SORTIE 
IF(XX(0)  .EQ.O.O)THF.N 

* 

*  2  KILLS  SUBTRACTED  OUT  &  EVENT  TIME  ADJUSTED 

* 

ATRIB(3)=ATRIB(3)-2 
ATRIB(  If  )=ATRIB(  1 6 )+DT IMF, P/60.0 
ELSE 

* 

*  2  KILLS  SUBTRACTED  OUT  &  EVENT  TIME  ADJUSTED 

* 

ATRIR( 3 )=ATRIR! 3 )-2 
ATRTR(  lf>)=ATRIB(  16)+DTIME/60.0 
ENDIF 

CALL  MISSILE 
IF(KUMKIL.LE.  S)THF.V 

DTMIS  1=DTMIS 1/00,0  +  0.0 

* 

*  ADJUST  ATTRIBUTE( 3 )  CONVOY  LENCTH  CORRECTION  INCLUDES 

*  KILLS  FROM  SORTIE .AVERAGE  INTERVAL  BETWEEN  VEHICLES 

*  ASSUMED! 2 5+50/2*3 7.5).  CONVERSION  MADE  TO  KM. 

*  SAME  APPLIES  FOR  TWO  EI.SF.  STATEMENTS  BELOW. 

*  EVENT  TIME  IS  ADJUSTED  DUE  DELAY. 

* 


ATR  T  R  (  3  )  =  ATR I B  (  3  ) -  KUMK 1 1, 

DELCH  1»ATRIB( 3 ) 

TF( ATRIB( 3 ) .LT . 0)ATRIB! 4 )=ATRIB( 4 )+ATRIB( 3 ) 
ATRIB( B )*ATRIR( B )-! (MUMKT L*3  7 . 3+2*37 . 5) / 1000.0) 
ATRI R(  lf>)=ATRIB(  lf>)+DTMIS  1 
ATRTB(  10)=2  3.SR7/ATRIB( 1  1) 

XX!  12)*  1.0 
VAL1«ATRIB( 1A) 

ELSE  IF!NUMKIL.LE.  10)THEN 
DTMIS2-DTMIS2/A0. 0+0.0 
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ATRIR(  3)=ATRIB(  3  )-NUUl'I  I, 

nr.i.C!  ]=atp.  i  p(  3 ) 

T F ( ATPTB (  3  )  .  LT.  o  ) ATRI P. ( 4  )=  VTR I P. (  4  )+ATR T  r, (  3  1 
ATRlP(R)=ATRTB(R)-(  (NtMKIL* 3 7. 5+2*17.5)/  1000. n) 
ATRIR(  16)=ATRI?f  16)+DTMIS? 

ATRIP( 10)=23.R97/ATRIR( 11) 

XX(  12)=  1.0 
VAL1=ATRIB( 16) 

ELS  r. 

DT'tir.3=DTMIS3/60. 0+0.0 
ATRIP. ( 3  )=ATRIR (  3 )-NUMKI L 
PFLOH  1=ATRIR( 3 ) 

IF  (ATRI  P(  3  )  .LT .  0)ATRIP,  (4  )=ATP.  IB(  4  )+ATRIP(  3  ) 
ATRIP.(P>)*ATRTB(R)-(  (NUMKI  I.*3  7 . 5+2*37 .5)  /  1000.0) 
ATRIR( 16)=ATRIR( 16)+PTMIS3 
ATP,IB( 10)=23.P97/ATRIB( 31) 

XX(  12)=  3 .0 
VAL3*ATIUB(  16) 

ENDIF 


*  TAIL  OF  CONVOY  A]  RESOH DOLED  WITH  NEW  EVENT  TIME. 


ELSE  IF((ATRIR(  l).F.O.  1.0)  .AND.  (  ATRIB(2  )  .EO.2.0)  .AND. 

+  <XX<  10).F.0.1.0).AND.(XX(  12).E0.1.0))T11EN 
ATRIP.  (  16)  =  (ATRIP(R)/ATRIR(  13))+VAL1 
ATRIR(3)=DELC111 
ATRIR(  10)=23. 987 /ATRIP, (  11) 

ENDIF 

RETURN 

1300  ATRIB(  10)=4.0/ATRIB( 1 1) 

RETURN 

1400  ATRIR(  10)=3.RR/ATRIB( 11) 

RETURN 

1500  ATR1R(  10)=1.P24/ATRIB( 3 1) 

RETURN 

1600  IF((yX(  10). EO. 0.0). AND. (ATRIP, (  3).0E.  1.0))TUF,N 
ATPIR(  10)=7 .24/ ATP.IB(  1 1) 

ENDIF 

IF((XX(10).EO.  1.0)  .  AND.  (  ATRIBf  1).GT.  1.0))TUEN 
ATRIB( 10)=4.74/ATRIB( 11) 

ENDIF 

RETURN 

1700  I F(  (XX(  lo).EO.  1.0).AN0.(ATRIB(  1).0T.  1.0))TI!EN 
ATRIB( 10)=2.5/ATPIB( 11) 

ENDIF 

TF((XX(  10). FO. 1.0) .AND. ( ATEIB(  1).FQ.  l.n).AND. 

+  ( ATRIP ( 2 ) .EO. 1 .0) )TUEN 
TVAI  T2*=F  F  AWT  (  2  ) 

ATRIR(  10)=2.5/ATP,IB(  11) 


fi.se  if((xx(  20).ro.  i.n)  .and.(  atrtb(  i).fo.  i.o).and. 

+  ( ATP.  I P.  (  2  )  .  T.o .  2  .  o )  )T  RE  N 
TV AIT2=FFANT(?  ) 

ATRIP (  10)=( 2 . 5+ATRIR(  PI) /ATi’.I P(  1  2) 

END  IF 
RETURN 

1800  ATRIP (  10)=]2.f'2/ATRIB(  11) 

RETURN 

★ 

*  UPDATE  ARMOR  VEHICLES  TOTAL, ONLY  TAIL  OF  CONVOY 

*  EXITS  I!E»E.  THIS  EVENT  TAIL  TIMF  ELEMENT  IS 

*  RECORDED  AND  GAFMCF  CALL  IF  LAST  ELEMENT. 

*  THIS  IS  LAST  OF  "IF"  FLOCKS  THAT  DENOTES  ORDER  OF 

*  FINIS!!  FOR  LAST  UNIT  IS  DEPENDENT  ON 

*  WHETHER  INTERDICTION  OR  NOT. 

* 

]O00  IF(ATRIP(  2)  .EQ.2.0)THEN 

ATRIB( 14 )=TNON-ATRIB( 23) 

IF(NKPl’N.EO.  DTHF.M 

WF.IT(5,'  (3I4.3FR.2)' )  l.N.NRUN,  POLICY,  ATRTB(  3)  ,ATRIB(  2), 
+  ATRIB( 24) 

F.NDIF 

XX(23)-XX(23)+ATRIR(3) 

TF((ATRIB(  1).E0.  25.0)  .AND.  (XX(  20)  .EO.O.O)  )THF.N 
XX(2B)=XX(2R)+1 


XX( 2R)=XX( 29)+TNOV 
CALL  GAFMCF 
F.NDIF 

IF( ( ATRIR ( 2).EQ. 1.0) ,AND.(XX( 10). EO. 1.0))THEN 
XX(2P)=XX(28)+2 
XX(  ID )=XX(  1 0 )+TNOM 
CALL  GAFMCF 
F.NDIF 
F.NDIF 
RETURN 

* 

* 

******************************************************************* 

* 

*  ROUTE  GREEN  (ENTITIES  16-31).  EXPLANATION  IS  THE  SAME  AS  GIVEN 

*  ABOVE .  HOWEVER, EVENT  20, AND  ENTER  NODES  4  &  5  FULFILL  THE  ABOVE  ROLES 

*  OF  1,2,  AMD  3, RESPECTIVELY.  COMMENTS  ABOVE  APILY  TO  SIMILIAB  AREAS  IN 

*  IN  GREEN  BELOW . 

* 

******************************************************************* 

* 

* 

2000  IF(XX(2)  .LF..3  1.0)THEN 
J1«INT(XX(2)) 

CALL  COrY(Jl,  1, ATRIR) 
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CAT.!,  ’’ESNUF2 
IF(  J  l.EO.  ]6)T!!EN 
TIME2=TNOW 
ATKIB(  13  )=TNOL' 

CALL  SCHDL(4,TIME2, ATRIP) 

F.  N'T!  IF 

ATRIB(2)=2.0 
TIME2=ATRIR( 8  )  / ATRIB(  11) 

Cl  NT2*ATRIB(  12) 

ATRI8(  R  )“64 . 972  /ATRIP.(  11) 

ATRIP  (  13)=TNOV? 

CALL  SCHDLC  5 ,TIME2 , ATRI B) 
irCJ  KF.0.3DGO  TO  2050 
J1*J  1+1 

CAT.L  C0PY(J1,  l.ATRIB) 

CALL  RESHUF2 

TIMF.2*(TIME2+(  CINT2/ATRI  R(  1 1))) 

ATRIB(R )=64 . 972 /ATRIP ( 11) 

CALL  SCHDL(4,TIME2, ATRIP.) 

ATRIP.(  13)-TIME2 

CALL  SCHDL(20,TIME2 .ATRIP) 

2050  XX(2)«XX(2)+1 

ENDIF 
RETURN 

400  CALL  ENTER(4, ATRIR) 

RETURN 

500  CALL  ENTER (5, ATRIP) 

RETURN 

900  CALL  ENTER(9 ,ATRIE) 

RETURN 

2100  ATRir(  10)»7.628/ATRIB(  1 1) 

PETl’RN 

2200  IF(  (XX(  10).EO.O.0).AMn.(ATRir.(  1).CE.  18. 0)  )THE.N 
ATRIB( 10)=2 . 7/ATRIP( 1  1) 

END  IF 

IF(  (XX(  10). EO.  1.0)  .AND.  (  ATRIR  (  1).F,0.  16.0))THEM 
ATRIP  ( 10)=2 .7/ ATRIB( 1 1) 

ENDIF 

IF((XX(  10). EO.  1.0).AND.(ATRIB(  1).CT.  ]6.0))TUEN 
ATRIR( 10)-17.42/ATRIR( 11) 

ENDIF 

RETURN 

* 

*  INTERDICTION  OF  UNIT  *16.  SEE  EVENT  12  APOVE 

*  SAME  THOUGHT  PROCESS  APPLIES  HERE  AS  ABOVE, 

*  WITH  NUMRER  OF  CONVOY  AS  EXCEPTION. 

* 

2300  IF( (ATRIB( 1).E0. 16.0) . AND. ( ATRIBf 2 ) .EQ. 1.0). AND. 
+  (XX( 10). EO. 1 ) . AND. ( XX ( 13).EO.O.O))THEN 
CALL  SORTIE 
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mxx(O)  .f/'.o.o)t;iej: 

ATRIP.  n)=ATRIn(T)-2 

.AT.R7P (  16)=ATRTB(  16)4DTI.N.rr/60.o 

else 

ATPIB(  3)=ATRIr>(  3)-2 
ATRIP (  16)=ATRIB(  16)4r>Tr<r./60.n 
ENDTF 

CALL  MI SPILT 
IF('UTX,KIL.LF.5)THEV 

OTMIS 1=DTMTS 1/60.0  +  o.n 
ATR I B  (  3  )= ATRI  B  (  3  )  -N IJMKI L 

IF( ATRIR  f  3 ) .LT . 0)ATRIB( 4 )=ATRTB( 4 )+ATPIR( 3 ) 
ATRIB(  B  )=ATRIB(  p.  )-(  (NUMXI  L*37 . 5+2*37 . 5)  /  1O00. 0) 
ATRIB(  16)=ATRTRf  ]6)4-DTMIS  3 
ATP  T  P.  (  10)=10.56B/ATRIR(  11) 

XY( J3>=  1.0 
VAL2=ATRIB( 16) 

ELSE  I F( NUMKI L • LE .  10)TECN 
DTMIS2-DTMIS2/60. 040.0 
ATRIB(3)=ATRIB(3)-NU>1KIL 
DF.LCH2=ATRIB(  3  ) 

IF (  ATRIP  (  3  )  .  LT  .  0  ) ATRIP (  4  )= ATRIP  (  4  )4ATP.I B  (  3  ) 

ATRI R(  8  )=ATRIB( 8 )-( (NUMKIL*37 . 542*3  7 . 5 ) / 1000 . 0) 
ATRIP ( 16 )=ATRIR( 16)4DTMIS2 
ATRIB(  10)= 10. 5 63/ ATRIP ( 1  1) 

XX( 13)=  1.0 
VAL2=ATRIP( 16) 

ELBE 

TVTM I S  3=r>TMIS  3  /60 . 0+0 . 0 
ATRIP ( 3 ) =ATR I B ( 3 )-NTMKI L 
0ELC02=ATR IP ( 3 ) 

TF(  ATP.IB(  3 )  .LT.0)ATRIB(4)=ATRIP(4)4ATRIB(3) 
ATRIB(B)=ATRIB(R)-((ML’MKIL*37. 542*37. 5)/  3000.0) 
ATRIB(  16)=ATRIB(  16)4DTMIS3 
ATRIB( 10)=10.56R/ATRIB( 1  1) 

XX( 13)=  1.0 
VAL2=ATRIB( 16) 

E'TDIF 

ELSE  IF( (ATPIB( 1).E0. 16.0) . AND. ( ATRIB( 2 ) .E°. 2 .0) . AMD. 
4  (XX(  10). EO.  1.0)  .AND.  (XX(  13).EO.  1.0))TUEN 

ATR I B (  16 )  =  ( ATP.TR(  R ) /  ATRIP (  3  1))4VAL2 

ATRIB(  3  )»DEI.C(I2 

ATRIB(  10)=  10. 568/ATRIM  1  1) 

F.NDIF 
RETURN 
2400  RETURN 

2500  ATRIP (  10)= 10. 568/ ATRI B(  11) 

RETURN 

2600  IF(  (XX(  10).E0.0.0).AND.(ATRIB(  1).CF..  16 .0)  )TI!EN' 

ATRIPT  10)«5.9/ATRIB( 1  1) 


ir(  cxxc  io). rc.  i.n)  .A-:n.r  atp.ifc  i).r/\  36.o) .and. 

+  fATh,Tp(2>.rf>.  1.0))TI!FN 
TUATT3=FFAMTC3) 

ATP.TR (  ]0)=3.5/ATRir.(  11) 

ELSE  IF(  CXXC  10).F.O.  1.0). AND. C ATRIRC  3).EQ.  16.0). AMD. 

+  (ATRIRC?). EO. 2. 0))THEN 
TWATT3=FFA,.,,T(  3 ) 

ATRIRC  10 )=(  3 . 5+ATP.IP.  (  8  )  ) /ATR1 R(  11) 

ENDIF 

I F  (  (  X  X  (  10). EO.  1.0)  .AND. (ATRIRC  1).CT.  16.0))TIJEN 
ATRIRC  1 0 ) =2 .  4/ATRIR( 1  1) 

ENT  IF 
RETURN 

2700  IFCCXXC  10)  .F.O.  1.0).AND.(ATRIP(  1 )  .  CT .  16.0))TH"M 
ATP.IB(  10)  =  2.016/ATRIF.(  11) 

EMDIF 

IF( (XX(  10). EO. 1.0).  AND.  (ATP.IRC l).EO. 16.0))T!IEM 
ATP.IRC  10)=2 .4/ATRIRC  11) 

EMDIF 

RETURN 

2800  ATP.IRC  10)=  1.344 /ATP.IRC  13) 

RETURN 

?°00  ATP.IRC  10)=12.292/AT°.IRC  1  1) 

RETURN 

30°O  ATRIRC 10)= 13 . 628 /ATRIRC 1  1) 

RETURN 

* 

*  SEE  EVENT  19  ABOVE, SAME  THOUGHT  APPLIES  HERE. 

* 

3100  IFC  ATRIRC  2  )  .F.O. 2 .0)TUEN 

ATP.IRC  14)=TN0W-ATPTR(  13) 

IFCMNRl'N.EO.  1)THEN 

'TRITE C 15, ' C3I4.3FP.2) ' )2,NNRUN, POLICY, ATRIB( 3) , ATRIRC  1) , 
+  ATRIRC 14) 

ENDIF 

XX( 2  4 )=XX( 2  4 HATP.IR ( 3 ) 

TF( C ATRIRC  D.E0.31.0) .AND.CXXC 10) .EO.0.0) )THEN 
XX(28)=XX(28)+1 

XXC20)=XXC20)+TNO'j 
GALL  CAFMCF 
ENDIF 

IFC  (ATRIp( 1 ) .EO.  16.0). AND.CXXC 10). EO.  1.0))THEN 
XX(28 )=XX( 28 )+ 1 
XXC2O)=XXC20)+TN0v' 

CALL  CAFMCF 
ENDIF 
ENDIF 
RETURN 


★  •A************************************:*:***************-********************** 

★ 


*  ROUTE  RFn  (ENTITIES  32-44).  EXPLAN  'TION  IS  TEE  SA’T  AS  riV-'F  A.T'''VF. 

*  • » OW  E  V  E  R ,  F  Y  F.  IT  32  AND  E1TEH  NODES  At  7  FI'LFI  1 1.  T."F  ABOVE  pn^rf.  or 

*  1,2, AN"1  3,  PEE  FEOTIVELY.  COMMENTS  AROVF  IF.  RUT.  APPLY  IF  REP  AS  T.U.. 

* 

* 

* 

3200  IF(XX(3).LE.44.0)THEN 

v  1  =  I  NT  (  X  V  (  3 ) ) 

CALL  C0PY(K1, 1, ATRIR) 

CALL  RF.SFFF3 
TF(K  l.F0.32)THFF 
TIME4=T  VO'.' 

A TRIP ( 13)=TFOV 

CALL  SC!'DL(  A  ,TIME4 , ATRIR) 

F.MPIF 

ATRIP.(2)=2.0 
CIFT3=ATRIR( 12) 

ATRIR(9)=56.62R/ATRIR(  11) 

ATRIP ( 33)=TMOW 
TIMF.4=ATRTB(R)/ATRIB(  1  1) 

CALL  SCI!PL(7  ,TIMF,4, ATRIP) 

IF(K1.F,0.44)CO  TO  3250 
k l-r 1+1 

CALL  C0PY(K1, 1, ATRIR) 

CALL  RESHUF3 

TIME4=(TIME4+(CIMT3/ATRIB(  1  1))) 

ATRIP (°  )=56.628/ATRIF>(  11) 

CALL  SCI!DL(6,TIME4  ATRIP) 

ATRIR ( 13)=TIME4 

CALL  SC!!PL(32,TIMF,4, ATRIR) 

3250  XXI 3 )=XX(3 )+ 1 
ENOIF 
RETURN 

600  CALL  ENTER(6, ATRIP) 

PE TURN 

700  CALL  F.NTEP(7, ATRIP) 

RETURN 

1000  CALL  ENTF.R(  10, ATRIR) 

RETURN 

3300  I F (  (XX(  10)  .F.O.O.O)  . AND. ( ATRIR (  1 )  .OF  .  32 .0)  )TMEN 
ATR1  B(  10)=  1  1.6  16/ATRIP(  1  1) 

ENDIF 

IF(  ( ATRIR(  1 ) .EO. 32.0)  . Af:P. (XX(  lO)  .EO.  1.0))RETURN 
IF(  (XX(  10) . F.0 .  1.0)  .  A*’0.  (  ATRIP(  1 )  .C.T.  32.0)  )Ti!FN 
A-rRin(  10)=(2  1 . 554+6 .048) /ATRIR(  1  1) 

ENOIF 


RETITV 


* 


*  interdiction  of  tv: it  *32.  see  r.vc'T  12  for 

k 


R  VPL*  %WTIO\' . 


3400  IF(  (ATRIP/  1)  .F.Q.32.0)  .  AND.  (  ATRT  R''  2  )  .FT.  l.n)  .A'.'O. 

+  (XV(10).E0.  1).ANX/XX(  14)  .EO.O.O)  )T*iEX 
CALI.  SORTIE 
IF(XX/9).EO.O.O)T!1EN 
ATRIP/3)=ATPIP/3)-2 
ATRIP, (  16 )=ATRIB(  16)+DTIMEP/6o.0 
ELSE 

ATRI  P,(  3  )=ATP.IR(  3  )-2 
ATP.IB (  16 )=ATRTR(  16)+OTIVF./60.0 
Ef'DIF 

CALL  MISSILE 
IF/KUMKIL.LE .  3  )THFM 

PTMTS 1=DTMIS 1/60.0  +  0.0 
ATP.IR(  3  )=ATP.TB(  3  )— fJIIMKI  L 
DELCH3=ATRIR( 3) 

IF (  ATRIP (  3 )  .LT . 0)ATRIB( 4  )=ATRIP( 4  )+ATRIP(  3 ) 
ATRTR(fi)=ATRIR/R)-((NUMKIL*37. 3+2*37. 5)/ 1000.0) 
ATRIP/  16)=ATRTR( 16)+DTMIS1 
XX/ 14)=1.0 
VAL3=ATRIB( 16) 

ELSE  I F ( M UMKI L . LF, .  10)TUEN 
DTMIS2=DTMIS 2/60. 0+0.0 
ATRIR( 3 )=ATRIR( 3 )-NUMKI  L 
PELCH3-ATRIB/3) 

IF/ ATRIB/ 3 ) .LT .0)ATRIB( 4 )=ATRIB( 4 )+ATRIP( 3 ) 
ATRIB (  B  )  =ATR I P, (  B  ) - (  ( NUMKI L* 3  7 . 3+2 *3 7 . 5  )  /  1 000 . 0  ) 
ATRIR(  16)=ATRIB(  16)+DTMir>2 
XX/  14)=  1 .0 
VAL3=ATRIB( 16) 

ELSE 

PTMIS 3=DTMIS3/60 . 0+0 .0 
ATRIR( 3 )=ATR IP ( 3 )-KUMKI L 
DELC!l3=ATRIP(  3  ) 

IF/ ATRIP/ 3) .LT.O)ATRIP(4)=ATRIB(4)+ATRIB(3) 
ATRIP/8)=ATRIB(8)-//Nt!.MKIL*37. 3+2*37. 5)/  1000.0) 
ATRIP/  16 )=ATRIB/ 16)+nTMIS3 
XX/  14)=  3 .0 
VAL3=ATRTB/ 16) 

F.NDIF 

ELSE  IF/  /ATP IB/  1)  .F.0. 32.0)  .AMD. /ATRIP./2)  .F.Q. 2.0)  .AND. 
+  (XX/  10).  EO.  1.0). AND. (XX/  14)  ,F.n.  1.0))THF.K 
ATRIB/ 16 )=/ ATRIB/ B) /ATRIB/ 11))+VAL3 
ATRIB/  3  )=r)F.LCH3 
THEAn3=ATRIB(  16) 

FXDIF 

RETURN 
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3  500  PETIT  V 

3500  iF((yy(  }0).rr'.o.r).A':f>.(ATf!in(  i).cE.32.o))T:ir.:- 
ATRI R(  H>)=4  .n/ATRTR(  1  1) 

ENDIF 

IF(  (ATRIP/  1).  Eo.  32.0).  a*.’0.  (yv(  10)  .m.  ].0))TUEN 
ATRIR(  10)=3 . 5  / ATRTP.  (  1  1) 

ENDIF 

IF(  (  ATRIP  (  !).GT.32.0).AND.(XX(  10). EC.  1.0))THEN 
ATP,  IT  (  10)=  1  1.72/ATPIR(  1  1) 

FMDIF 

PETIT': 

3700  ATRIP*/  10)=  1 1 . 72 /ATRIP  (  11) 

RETURN 

3R00  ATRIP- (  10)=  10.3/ATRIR/  11) 

RETURN 

3900  ATRIP.(  10)=2 R . 0/ATRIB/  11) 

RETURN 

* 

*  SEE  EVENT  19  A  ROVE ,  SAME  TUOUCTU  APPLIES. 

* 

4000  IF(ATPIR(2) .F0.2.0)THFN 

ATR I R ( 14)=TNON-ATRIR( 13) 

TF/KNRUN.F.O.  1  )T!1EN 

WRITE/  1 5,  ’  /  31 4 , 3F8. 2  )  '  )3  .N.NRUN,  FOLIC  V,  ATRTP-C  3  )  ,  ATR  I B  (  1 )  , 

+  ATP  IB/  14) 

ENDIF 

XX/2  5)*XX(25)+ATRIE/3) 

IF( ATRIR ( 1).EQ.44.0)THEN 
XX ( 2P)=XX(2S)+1 
XX(2  1)=XX(2 1)+TK0W 
CALL  CAFNCF 
ENDIF 
END  IF 
RETURN- 
END 

* 

* 

*********************************************************************** 

* 

*  ROUTE  FLIT..  THE  FOLLO'.JINC  THREE  SUBROUTINES  PREVENT  UNITS  FROM 

*  PASSING, EACH  OTHER  ONCE  ON  THE  ROUTE  OF  MARCH.  THE  LEAD  UNIT 

*  SETS  THE  COLUMN  PACE.  OTHER  UNIT  RATES  ARE  ADJUSTED  DEPENDING 

*  ON  IF  THEIR  RATE  IS  LESS  THAN  OR  EOUAL.OR  GREATER  THAN  THE 

*  LEAD  UNIT  RATE.  THIS  ASSUMES  UNITS  MAINTAIN  SAME  RATE. 

* 

*********************************************************************** 

* 

* 

SUBROUTINE  RESHUF1 

COMMON/ SOON 1/ATRIR( 100) ,DD( 100), DDL/ 100) ,DTNOV , I I , MFA, MSTOP, NCLNR 
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+  ,  ncp.pr ,  r:  p* nt ,  r: nrun ,  nns r.i ,  nta  rr: ,  s s n  no ) ,  s  si.  ( l on ) t  tnext  ,  r .w ’ , 
+  >: x c  ion) 

co-'MO'.Vroo-'  i/tim^time  ],TF?E.'>,'r-i.ME3,TiMr4,TiuE5,  policy,  r^t? 
+  ,T2,n,j  i,  J2,  J3,nir  ],DiF2,DiF3,nTiMF,rTTurr,r.TM-'iL,PT.Mis], 

+  PT’IS2,RTMIS3,I,  .1 
SAVF./UCOM  1/ 

P HAL  TIME, TIME  1  ,TIME2  .TIMF.3  ,T IME4  ,TIME5 , PI F  1 , DIF2  ,PIF3  , 

+  PTIME  ,  DTIME  P, DTM IS  1 , DTM  I S2  ,  DTMIS3  , TIMEM 

INTEGER  I, I 1,I2,I3,J,J1,J2,J3,F,L,M,N, POLICY, ROUTE, 

+  KUYKIL.Kl'rLAN.MSilOT 


* 

* 

IF( ATRIP ( 1) .FO. 1)THFN 
XX(4 )=ATRIB( 11) 

FLSF.  IF(  (ATRIR(  1).CT.  1)  .  AMD.  (  ATRIB(  1 1 )  .IF .  XX(4  )  )  )T!!F.N 
PIF  1=XX(4 )-ATRIR( 1 1) 

ATP  I P> (  1 1)=ATRIP.(  1  1  )+PIF  1 

ELSE  IF(  (ATRTR(  1).CT.  1 )  .  AND.  (  ATP.IR(  1 1)  ,0T.XX(4)  ))THFN 
DIF  1=ATRTR( 11)-XX(4) 

ATRIR( 1 1)=ATRIF( 1 1)-DIF 1 

fndtf 

RETURN 

END 

* 

* 

♦ROUTE  GREEN 

* 


* 


SUBROUTINE  RESHUF2 

COMMON/SCOM 1/ATRIR( 100) ,DD( 100) , DDL ( 100 ) ,PTNOU , I I ,MFA, MSTOF, NCLNP 
+  ,  NC RDR  ,  N PRUT , NNRUN ,  NNT. ET , NTA FF. ,  S S (  100)  ,£SLC  100 )  , TNEXT ,  TNON  , 

+  XX(  100) 

COMMON/l’COM  1 /T IMF , TIME  1  ,TIMF.2  ,TIME3 ,TIME4 , TIME 5 , POLICY ,  FO'JT E 
+  ,12,13, J1,J2,J3, DIF  1,DIF2,DIF3, DTIME, DTIMEP,NUMKIL, DTMIS 1, 

+  DTMIS2 ,PTMIS3 , I, J 

SAVE/UCON 1/ 

REAL  TIME, TIME  1  ,TIME2 ,TIME3 ,TIME4 , TIME5, DIF  2, DIF2.DIF3, 

+  DT  IMF. ,  DT  IMF  P,  DTM  I S  1 ,  DTM I S 2  ,  DTM IS 3 ,  T IMEM 

INTEGER  I,U,I2,l3,J,Jl,J2,J3,K,L,M,y,  POLICY, ROUTE, 

+  NUMNI L, NU  PI,  AN ,  MS  HOT 

* 

* 


TF(  ATRTRf  1).E0.  16)THF,N 
XX(  5  )=ATPvIB(  1 1) 

ELSE  m(ATRIR(  3).CT.  16)  .AND.  (  ATRI*(  1  1 )  .LE .  XV(  S  )  )  )TI!EN 
DIF2=XN(3)-ATRTR( 11) 

ATRT  P>(  1  1)=ATRIR.(  1 1)+DIF2 

ELSE  IF((ATRIR.(  1).GT.  16)  .  AND.  (ATRIP.  (  1 1)  .OT. XX( 3) ))THEV 
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DTF2=ATRIR(  1  1)-XV(S) 
ATFT?(  1  l)  =  ,\Tnm  1 1)-D1F2 
'•Noir 
urTimr: 
end 

*■ 

* 


* route  red 

* 

* 


* 


* 

* 


* 

* 

* 

* 

* 


+ 

+ 

+ 

+ 


+ 

+ 


SUBROUTINE  RESULTS 

COMMON/SCOM  I/ATIIIBI  100),DD(  100), DDL I  100)  ,DTN0-',II  ,MFA, MSTOP.NCLNR 
,  NC  RPP. ,  K  PRUT ,  NNRUN ,  NNSET ,  NT  APE ,  SS  (  100)  ,SSL(  100)  ,TiNEXT,TNOU  , 
XXfl^O) 


OOMMON/UCOM 1/TIME.TIME 1 ,TIME2 ,TIMC3 ,TIME& .TIME? , POLICY, ROUTE 
, 12 , 13 , J  1 , J2 , J3 ,DIF 1 ,DIF2 , PIF3 , DTIME , DTIMEP, NUMKI L,DTMIS 1 , 
DTMIS2 ,DTMIS3 , I , J 
SAVE/UCOM 1/ 

REAL  TIME, TIME  1 ,TIME2 ,TIME3,TIME4 , TIME 5, DIF  1,DIF2,DIF3, 

DTIME , PTIMFP, DTMIS 1 , DTMIS2 ,DTMIS3 , TIMEM 

INTEGER  I,I1,I2,I3,J,J1,J2,J3,K,L,M,N, POLICY, ROUTE, 

NUMK I L ,  N  U  PLAN ,  MS !  !OT 


IF(  ATRIP>(  1).EQ.32)TIIEN 
XX(6  )=ATRIB(  1 1) 

ELSE  IF(  (ATP  IP.  (  1).GT.32).AND.(ATRIR(  1  1)  .LE  .XX(6)  )  )T!IEN 
DIF3*XX( 6 )-ATRIB( 11) 

ATRIB(  1 1 )=ATRIB( 1  1)+DIF3 

ELSE  I F (  (ATRIB(  1)  .CT.  32)  .AND.(ATRIR(  1 1 )  .GT. XX(6 )  ) )TI!EN 
DT^3=ATRIB(  1  l)-XX(f>) 

ATRIB( 1 1)=ATRIB( 11)-DIF3 
END  IF 
PETURN 
END 


SUBROUTINE  SORTIE 

COMMON/SCOM  1/ATRTP.(  100)  ,DD(  100), DDL!  100)  ,  DTNOW  ,11,  MFA,  MSTOP,  NCLNR 
+  , NCRDR , NPRNT ,NNRUN, NNSET, NTA  PE , SS( 100),SSLC 1 00 ) , T  NE  XT , T  NON , 

+  xx( loo) 

COMMOJI/UCOM 1/TIMF, time ],TIMF.2,TIME3,TIMF4, TIMES,  POLICY, ROUTE 
+  ,  12 , 13  ,.1 1 ,  J2  ,  J3  ,DIF  1  ,DIF2  ,  DIF3 ,  DTIMF. , DTIME P,NUMEIL,DTMTS  1 , 

+  DTMIS2 ,DTMI S3 , I , J 
SAVE/UCOM 1/ 

REAL  TIME, TIME  1  ,TIME2 .TIMF.3 ,TIME4  , TIMES, DIF  1  ,DIF2 , DIF3  , 
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r 

v 


+  PTrMF.  .P^TUF  r,PT'!IS  1,  PTMir2  ,DT''1S3  .TIMEM 

TNTrOFR  T,1  1,1 2,13,.!,  J  l,J2,.n,l’,L,‘V:.  POLIO Y.^OUTE , 

+  ?:r‘,‘’T!./:rn.'‘:,vrnr'T 


xx(Q)=n.n 

i:i’PLAV-vv(7) 

IF(NUrLAN.LE.O.r\)TUF.N 
CO  TO  51 
ENDIF 

* 

*  DELAY  TIME  OMPUTF.P  nEPENDlNO  OX  ROAD  DAMAOE . 

* 

DTIMF.=TRIA0O0.0,34.0,3P.0,5) 

PTIMF.F=TRIAC(  10.0,  12.5,  15.0,  A) 

DO  50  N= 1 , 2 

* 

*  FANDOM  #  DRAWN  AND  TESTED  BELOW.  IF  RN  .LE.  0.35 

*  ROAD  IS  CUT, OTHERWISE  ROAD  NOT  CUT, PUT  PARTIAL  RLOCU. 

* 

PA  SS=" MF  RM( 0 . 0 , 1.0,2) 

WRITE(  11, ’ (I4,FR.2)’ )NNRUN, PASS 
IF( PASS.LE .0,35)THEK 

xx(o)=xx(R)+i 

ELSE 

XX(9)=XX(9)+0.0 

ENDIF 

XX(° )=XX(9 )+0 . 0 

50  CONTINUE 

NU  PLAN=NU  PLAN  -  2 
XY(7)«NTFLAN 

51  RETURN 
END 


SUBROUTINE  MISSILE 

COMMON/SCOM 1/ATRIB(  100) ,DD( ]O0) ,DDL( 100) , DTNON , II , MFA.MSTOP, 

+  , nc rdr , n fr ?rr , n nru n , n ns f.t , nta re , s s (  ioo),ssl(  ioo) ,tnext,tnow, 

+  XX( 100) 

COMMON/UCOMl /TIME, TIME  1  ,TIMF.2  ,TIME3  ,TIMF,4  .TIMF.5  ,  POLICY, ROUTE 
+  , 12 ,13, J  1,J2,J3,DIF 1,DIF2,DIF3,DTIME, DTIMET, NUMEI L , DTMIS 1 , 

+  DTMIS 2 ,DTMIS3 , I , J 
SAVE/UCOM 1 / 

t?E\L  TIME , TIME  1  ,TIMF.2  ,T IMF.3  ,TI’,r.4  ,T IMF 5 , PIF  1 , DIF2  ,PTF 3 , 

+  DTIMF., DTIMFP, DTMIS  1 ,  DTMTS2  ,  DTMIS3  .TIMF.U 

INTECEP  1, 1 1,1 2,13, J,J1,J2,J3,K,L,M,N, POLICY, ROUTE, 


+  N  UMl  I L ,  U  PL  A  r  ,  MS MOT 


VINT=0 .0 
NUMKIL=0 
NUMISS=XX(R) 

IF(NUMTSS.LF.O.O)TUFN 
CO  TO  76 
ENDIF 

DETEKMI NE  DELAY  TIME  DEPENDING  0*:  LEVEL  OF  DAMACE. 

DTMIf,  ]=TRIAC(  15.0,  IE.  5, 25.0,8) 

DTMIS2=TRIAC(22 . 5,30.0,37.5,9) 

DTMIS3=TRIAG(30. 0,4 1.5,50.0,  10) 

MSHOT*XX( 11) 

VEHICLE  INTERVAL  COMPUTED  FOR  EACH  UNIT , ADJUSTMENT 
15  MADE  FOR  INTERVAL  RETUEEN  COMFANIESf  6*3 7 . 5  f.  4*37.5) 

ONCE  AGAIN  AN  AVERAGE  INTERVAL  IS  USED  TO  SUBTRACT  OUT 
INTERVAL  FOR  COMPANIES. 

IF(ATRIP.(  1).E0.  1.0)Ti!EN 

VINT=( ( ATRIR( 8 )* 1 000.0-6*37. 5 ) / ( ATRIB( 3 )+ATRIB( 4 )))/0.304° 
ELSE  IF(ATRIR( 1).E0. 16.0)THEN 

VINT=( ( ATRIB(8)* 1000.0-6*37. 5 ) / ( ATRIB( 3 )+ATRIB( 4 ) ) ) /O . 3048 
F.LSF.  IF(  ATRIR(  1)  ,E0. 32.0)T11EN 

VI NT=( ( ATRIB( 8)* 1000.0-4  *37. 5 )/( ATRID( 3 )+ATRIR( 4 ) ) ) /O . 3048 
ENDTF 

DO  75  N= 1 , MS HOT 

RANDOM  *  DR A1  'N , DE FE NDINC  ON  VEHICLE  INTERVAL 
APPROPRIATE  RRANCH  TAKEN  f.  NUMBER  OF  KILLS 
COMPUTED. 

SHOT=UNFP.M(  0.0,  1.0,7) 

UmITE( 11, ' f I4.FR.2) ' )NNRUN,SHOT 
I F ( V I NT . LE . 120.0)THEN 

IF ( SHOT .LE . 0 . 1074 )THEN 
NU>'KIL*NUMKI  L+0 
ELSE  IF(SH0T.LE .0 .  1 142)THF.N 
NUMKTL-NUMKIL+1 
ELSE  IF( SHOT .LF. . 0 .  1260)THEN 
NUMKI L=NUMK I L+2 
ELSE  IF( SHOT . LE . 0. 1498)THEN 
NUMK I L=NUMK I L+3 
FLSE  I F (  SHOT .LE.O.  1836)TI!EN 
NUUKIL^NUMKI  L+4 
FLSE  I F( SHOT .LF .0 . 2420)THEN 


UUMNI  !.='T‘l'I  1+5 

else  if(seot.le.o.3472)tuen 
uumeii  .=h’mnj  l+a 

ELS E  IF (Siuvr.  Lr  .0.  5754  )THEN 
::umkil=numkil+7 
F.LSF 

kumkil-numkil+s 

endif 

ELSE 

IF( SHOT. LF . 0. 1074 )THEN 
Ml’MKI  L“NUMKI  L+0 
ELSE  IFCSnOT.LE.O.  149S)TEEN 
NUMKI L*NUMX I L+ 1 
ELSE  T F (  S HOT . I.E  . 0 . 2 302  )T ME N 
nitmxil*numkil+2 

ELSE  IF( SHOT. LE. 0.4 5<>2 )Ti!EK 
Kl1MKIL*FL’MKIL+3 
ELSE 

NUMK I  L=N  UMK I L+4 
ENDIF 


rvpjp 


NUMKI L-NUMEIL+O 

75 

CONTINUE 

NU*ISS»N 

,rv1T  SS*Mf»!!OT 

XX(R)-NU 

MISS 

76 

RETURN 

END 

* 

* 

* 

* 

* 

susroitine  cafmgf 

COMMON/SOOM  1/ATRIP- (  100),D0(  100)  ,  DDL  (  1 °0 ) , PT  NOW , 1 1 ,  M  FA , MS  TO  F , NO  LNE 
+  ,  K'CPPR ,  NFELT,  NNPUN,  NNSET ,  NTAFE  ,  SS(  100),SSL(  100)  .TNEXT.TNOU, 

+  XX(  100) 

COMMON/UCOM 1/TIME, TIME 1,TIME2,TIME3,TIME4, TIME  5, POLICY, ROUTE 
+  , 12 , 1 3 , J  1 , J2 , J3 , PIF 1 , DIF2 ,DIF3 ,DTIME , PTIMF.  P,  NIJMKI L,  DTMI S  1 , 

+  PTMIS2 ,DTMIS3 , I , J 

SAVE/UC0M1/ 

REAL  TIME, TIME  1,TI‘?E2, TIMF.3, TIME4, TIME 5, PIF  1,PIF2,PIF3, 

+  DTIME , PTIMEP, DTMIS  1 , DTMIS2  , DTMI  S3 , TIMF.M 

INTEGER  I,I1,I2,I3,J,J1,J2,J3,K,L,M,N, POLICY, ROUTE , 

+  NUMKI L , MU  FLAK , MS  HOT 

* 

* 

* 

IF( ATRI R( 1) .EO.  15.°)THEK 
XX(22)*XX(22)+XX( 10) 

XX(26)=XX(26)+XX(23) 
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rLEF  TF( ATRIR(  1).E0.  ].0)THEN 
>:x(22)-yy(22)+>:v(  19) 
x  x  ( ?  a  ) = x  x  ( 2  g  ) + x  x  ( 2  3 ) 

F.LSF  IF(ATrTP.(  1 )  .FQ . 3  1 . 0 )THEN 

yy(22)=yy(22)+xy(20) 

xy(26)=yx(26)+yy(2A) 

ELSE  IF(ATRir,C  1).E0.  lf>.0)THEM 

yy(22)=yy(22)+xx(20) 

XX(26)=XX(26)+XX(24) 

ELSE  IE(ATRTR(  1 )  .F.0.44 .0)THF.N 
XX(22)*XX(22)+XX(2 1) 

XX(? A)=XX(26)+XX(2  S) 

FNDIF 

IF(XX(2B).F/\3.0)THEN 

XX(27)=XX(27)+(XX(2f>)/(XX(22)/3) ) 

* 

*  APPROPRIATE  GLOBAL  VARIABLES  UPDATED.  REASON  FOR  UNIT 

*  15&1,  10&31.44  IS  THAT  15M.1AS31  FINISH  IN  DIFFERENT 

*  ORDER  DEPENDING  ON  WHETHER  INTERDICTIO!:  OR  NOT, 

*  WHEREAS  44  ALWAYS  FINISHES  LAST  DUE  TO  32  NOT  BEING 

*  DELAYED  LIKE  THE  OTHER  2. 

* 

'•'RITE(23,  ’(214,  SF«.  2)’  )NNRUN,  POLICY,  XXns),XX(  lf>),XX(  17), 
+  XX  ( IS), XX (27) 

* 

*  PART  OF  TAIT.  USED  IN  MODEL  VEP.F  I  CATION. 

* 

WRITF.(  19 , '  (T4  )  '  )NNRUN 

WRITE (  19, * ("FINISH  TIMES  " , 3FR . 2 ) ' )XX( ]0 ) ,XX( 20) , XX( 2 1 ) 
WRITE (ID,' ("AVERAGE  TIME  " ,FR.2) ’ )XX(22 )/3 
WRITE( 19, ' ("ARMORED  VEHICLES  ”,F8.2)')XX(26) 

END  IF 

RETURN 

END 
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The  model  presented  here  was  used  to  get  an  idea  of 
how  to  obtain  a  percentage  figure  of  when  the  road  could 
be  expected  to  be  cut.  After  1000  simulated  bomb  drops, 
the  number  obtained  was  0.35  for  the  parameters  listed  in 
the  program.  After  obtaining  it,  experienced  nersonnel  were 
questioned  about  its  appropriateness.  Based  on  the  general 
consensus  of  those  questioned,  this  number  was  used  as  the 
probability  of  cutting  the  road  in  subroutine  SORTIE  within 
STOPEM . 

The  author  gratefully  acknowledges  the  expert  assis¬ 
tance  of  LTC  Ivy  Cook  in  the  preparation  of  this  model. 


>v.v 


program  cookf.m ( i  n put  ,  out  rur ,  ta  pt,  l ) 

COMMON/UCOM  1/A(  4000) 

peal  a ,  r !  :eta  ,  i  ,  c  r  rr, ,  nr  nd  i ,  vent  i ,  r.  t  cm a::  ,  s  i  c:  •  ay  ,  >:  pr  i  me  i  ,  y  p>  i :  • r.  j , 

+XDIST 1 , YDIST 1 ,XDTST2 , YDIFT2 , XRND2  ,  YP :  ND2  , ND  IT .  I* !  1 ,  YD r>  I’ll,  RRJM2, 
+YD PR  IM?  , X  PR IME2 ,  Y TP. IMP 2 , Y3 , K , X 1 , X2 , Y  1 3 , Y2  3 , 1’AFTCT L , H  AFTCT’ 

integer  n,  k  i  ,  c  ] ,  C2  ,C3 ,  C4  ,c 5 

DATA  A/4000*0.0/ 

TlinTA=30.0 
UR=15.0 
CEPP=  120.0 

N 1=  1 

p=0.0 
01=0 
02=0 
03=0 
C4=0 
05=0 
K.=2  .0 
Y  13=0.0 
Y2 3=0.0 
XDPRIM  l=-30.o 
YD PRIM  1=90.0 
HAFTO,TL=60.0 
HAFTCT"=2  5 . 0 


*  XDrRIMl  &  YDRIPM1  ARP.  the  upper  left  coordinates  of  THE  ROMP 

*  SOX, CORRESPONDING  to  the  2ND  OU ARRANT.  THE  !L\FTGTL  s.  haftgtw 

*  ARE  THE  DIMENSIONS  OF  HALF  OF  THE  ROAD  LENGTH  AND  WIDTH  IN 

*  FEET.  IMSL  SUP, ROUTINE  IS  GALLED  TO  GENERATE  4000  RANDOM 

*  VARIATES  USING  THE  SEED  GIVEN  PF.LOW  AND  PLACING  THEM  IN 

*  ARRAY  A, DIMENSIONED  TO  (1,4000).  A  TRANFORM AT ION  OF 

*  AXIS  IS  MADE  TO  THE  NEW  COORDINATE  SYSTEM,  TN  RELATION 

*  TO  THE  ROAD.  CEP  FORMULA  USED  I"  THIS  FORMULATION  WAS 

*  OSTAINED  FRO’’  (REF  10).  THE  H  VALUE  RF.  FEES  ENTS  THE 

*  RATIO  OF  PEP  TO  DET,  THAT  IS  THE  RET  =  5*PET  DUE 

*  TO  THE  AIRCRAFT  VELOCITY  CONTRIBUTING  THE  GREATEST 

*  ERROR.  THE  ANCLE  OF  ATTACK  USED  IN  THIS  RUN  WAS  30  DEGREES 

*  WHICH  SAYS  THAT  THE  AIRCRAFT  ATTACK  30  DEGREES  OFFSET  FROM 

*  THE  ROAD  AXIS. 


GALL  GGNML(  123457. DO, 4000, A) 

Xl-XD PRIM  I*COSD(THETA)-YD PRIM 1*SIND(THETA) 
Y 1=XD  PRIM 1*SIND(THETA)+YDPRIM 1* COS D( THETA) 
X2*XD  PRIM l*COSD( THETA )+YDPRTM 1*STND(THETA ) 
Y2=X0  PRIM 1*SIND(THETA)-YDPRIM 1*C0SD( THETA) 
SICMAX=GF.rP,/(0.6512  +  K*0.5A4O) 

SIGMAY=K*S ICMAX 

Al=(XnrRIM  1*2 .0)*C0SD(THETA)+HArTGTW 


OP 


QwiwWi! 


IF  THE  ANGLE  IF  LESS  THAN  l.THEN  THE  NEXT  COMPUTATION  IS 
P.YPASSrn  BECAUSE  Till"  ARGUMENT  BLOWS  t'p. 

IF(ARSfTHETA-90.n).LT. l)CO  TO  ] 
P=(HAFTCP.,/XnrRTMl)*TAr:P(T!'ETA)-!:AFT!TL 
]  pr  i  nt* , 

PRINT*, 

PRINT*, 'K  IS 
PRINT*, 'CEP  IS  ' ,  CEPS 
PRINT*, 'THETA  IS  '.THETA 
PRINT* , ' A  1  ’ , A 1 

PRINT*, 'R  ' , R 

PRINT*, 'BOMB  BOX  LENGTH  IS  ’  ,  1RO 

PRINT*, 'BO"R  BOX  WIDTH  IS  ’,60 

PRINT*, 'TARGET  BOX  LENGTH  IS  ' , 2 .0*HAFT0TL 

TUI  NT*,  'TARGET  BOX  WIDTH  IS  '  ,2 .0*HAFTGTW 

PRINT*, 

PRINT*, 

1000  BOMB  DROTS  ARE  SIMULATED. 

no  100  N= 1,1000 

XPvNDl*A(!U)*SIGMAX 
Ml-N  1+1 

YRND 1=A(K !)*SIGMAY 
N 1*11 1+1 

XRKD2*A(N1)*SIGMAX 

Kl-Nl+I 

YRND2=A(N  1)*SICMAY 
N1«N  1+1 

BOMB  PROPS  ADJUSTED  DUE  TO  RANDOM  I M TACTS. 


X2DIST 1  = 
X  1DIST  Is 
Y2DIST 1= 
Y1DIST  Is 
X2PIST25 
X  1PIST2* 
Y2DIST2* 
Y1DIST2= 


X  1+XRND2 
X  1+XRND1 
Y 1+YRN02 
Y 1+YR  ND 1 
X2+XRND2 
X2+XP.ND1 
Y2+YRND2 
Y2+yrt:d  1 


CHECK  MADE  TO  DETERMINE  IMPACT  OF  BOMBS,  TESTED 
AGAINST  THF.  ROAD  DIMENSIONS. 


IF( (X  1DIST 1 . CT . -HAFTGTW) .OR.(X  1DIST2 .LT . A 1) )THEN 
01=0  1+1 
GO  TO  3  SO 

FNDIF 

IF (ABS( THETA) .LT . 1.0)THEN 


300 


IFUYlDIf.Tl.CT.  -OAF  TOT  L)  .  AMP.  (  Y  TTST2  .  LT.  HAFTCTi. '  'T  Hr 
Y]  3=o.O 

Co  to  300 

else 

00  TO  350 

ENDIF 

F.NPIF 

IF(AEB(TUETA-oo.0) .LT.  1)T:>EN 

rp.or>!=YioT5T,]-2.o*>:opaiM  ] 

TF(  (Y  1PICT 1  .LT.HAFTCTL)  .AND.  (  PROF  1  .CT.-I!ArTCTL)  )7HEN 
C2=C2+1 
CO  TO  100 

F.T.SF. 

CO  TO  350 

n>;ni  f 

ENDI*- 

Y  13=Y  1DIBT  l+( (Y1DTST2-Y3DIST1)/(X 1PIST2-X 1DIST1) ) 
*(-HAFTCTV-X  1PTFT1) 

TF(  (Y13.GE.P)  .A\’0.(Y  13.LE.PAFTCTL)  )7HEN 
C2=C2+ 1 
CO  TO  100 
END  IF 


*  THIS  PART  RE  PRESENT?  TEE  SECOND  AIRCRAFT  ATTEMPT  TO 

*  CUT  THE  P.OAO.  PRINT  STATEMENTS  AMP  COLTER  VARIABLES  ARE  PEEP 

*  TO  KF.EP  TRACK  OF  N’lRtprp  OF  CHTS  AMD  MISSES. 

* 

350  CONTI  MU  F. 

TF( (X2PIST1.GT. -HAFTCTV) .OR. ( X2PIST2 .LT . A  1 ) )THEN 
C3=C3+1 
CO  TO  450 

EKDIF 

IF  I ARS (THETA)  .LT.  1.0)T!1EM 

TF(  (Y20IST1  .CT.-HAFTCTL)  .ANP.(  Y20TST2. LT.HAFTCTL’!  )T!LCM 
Y2 3=0.0 
CO  TO  400 

ELSE 

00  TO  4  50 

F.NDIF 

F.NPIF 

IFfAPSfTHETA-OO.O)  .LT.  1  )THrr: 

FP.on2=Y2niSTi-2.o*>;nrr  i-i 

IF(  (  Y2DIST1.  LT.HAFTCTL)  .AMO.  (PR  002.  CT.-HAFTCTL)  )Ti'EM 
C4=C4+ 1 
Co  TO  100 

ELSE 

Co  jo  450 

EHniF 

F.MPTF 

Y23=Y2PISTl+(  fY2niST2-Y2niFT3)/(X2nir.T2-X2Pl?Tl)) 
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4  r \  ■  \ 


*  ( -• ' \v~r'r’ v 21'  J S T  1 ) 


I r ( ( Y2  3 .  r.r r.> .  a-t .  O’?  3 .  it. .  "■ 
04=C4+] 
rn  jo  ]on 

pr'DIF 


450  CONTI NPF 


CSsC'S+l 


ion  CONTI  NIT. 
FRTFT* 
PRINT* 

FR  IMT* 

PT  I  'IT* 

PR  T  N’T* 

pp  j  »;t* 

FTII'II* 
PRINT* 
PRINT* 
PRINT* 
ITT  NT* 
PRINT* 
END 


•  ft 

TIMES 

mot 

Cl’T  DUE 

XP  y-DIPECTION 

*1 

'  ft 

TIMES 

NOT 

CPT  PIT. 

to  x-niRF.CTio:; 

if  2 

»  1: 

TIMES 

cn 

P  v  it  ] 

’  ,C2 

•  ft 

TIMES 

CUT 

F.Y  *2 

\C4 

’TOTAL  CITS  ARE  ’  ,02+0.4 
’ !•  TIMES  NnT  CPT  AT  ALL  ’  ,C5 
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APPENDIX  C 


THREE-WAY  AND  FOUR-WAY  ANOVAS 


The  SPSS  program  for  the  four-way  ANOVA  is  listed  on 
the  next  page.  ANOVAs  are  listed  in  the  following  order: 

(1)  four-way  with  interactions;  (2)  four-way  without  inter¬ 
actions;  (3)  three-way  with  interactions;  and  (4)  three-way 
without  interactions.  The  SPSS  programs  shows  the  ANOVAs 
excluding  the  variable  UINT,  the  all  other  units  interval, 
because  the  reconnaissance  interval  also  varied  by  the  same 
amount.  By  running  SPSS  with  both  of  these  factors,  a 
singular  matrix  resulted  because  both  these  columns  represent 
the  same  thing.  Excluding  this  factor  gave  the  ANOVA  table. 
Tape  1  is  a  complete  listing  of  all  256  simulation  runs.  The 
columns  from  left  to  right  on  tape  1  represent  run  number, 
policy,  length,  rate,  reconnaissance  units  interval,  all 
other  unit  intervals,  and  measure  of  effectiveness. 


mmm 


m 


>'■ 


Kit 

>V.- 


.  *  « 

.A 

i'S\ 

i 


l 

% 


v ' 
V 
r. 
v 

L*  •, 

* 


H 


i 

& 


TM'»*  v,'.\Vff[ 

rAcrsi7r 

cj^r  »j^vr 

VARIABLE 
V*.t’  LAPP], 


.1ST 


recode 
r  eco  nr. 
rf.coop 
peco  nr. 

V  or  CAFES 
INPUT  FORMAT 
LIST  CASES 
ANOVA 

STATISTICS 
REAP  INPUT  DATA 
FINIS!! 


urn  interdiction  analysis 

50 

TAPE  H INTERDICTION  RESULTS) 

POLICY, LENGTH, RATE ,RI  'T ,  ''I NT,  VP  HEAT 
POLICY, INTERDICTION  POLICY/ 

LENCT!', CONVOY  LE”CT!!/ 

TATE, CONVOY  RATE / 

PINT, RECON  UNITS  INTERVAL/ 

PINT, ALL  OTHER  UNITS  INTERVAL/ 

VEHP.AT, ARMORED  VEHICLES  JT.R  HOUR/ 

LENCTH(o=  ] )  (0.  ]=2  ) 

PATE ( 2 5.0=1)  (22.5=2) 

RINT(2 5 .o= ] )  (22.5=2) 

UINT(4=1)  (3.6=2) 

UNKNOWN 

FiyED(4x,F4.0,5F8.2) 

CASF.S=20/VARIARLES=ALL 

VEHRAT  BY  PO LI C Y (  ] ,  4  )  , LE NOT ” (  1,2)  RATZ(  1  *>) 
RINT( 1,2)  *  ’  ’’ 

ALL 


l*--< 
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vVft.VUih.No.V.V 


'  rl 


******* 

*  A  n  a  1 

.  Y  S  T  ?  0  F  V  A 

R  T  ANTE  * 

******* 

VEHRAT 

ARMORED  VEHICLES  PER 

nop 

ry 

POLICY 

INTERDICTION  POLICY 

LENGTH 

CONVOY  LENCTH 

RATE 

CONVOY  RATE 

RINT 

RECON  UNITS  INTERVAL 

******* 

*  *  *  *  * 

*********** 

******* 

******* 

SUM  OF 

MEAN 

SI CN I F 

SPURGE  of  variation 

SOU A RES 

DF  SQUARE 

F  OF  F 

main  fefeots 

5031.203 

6  838.534 

88.481  .OH] 

rOI.TCY 

3  393 .753 

3  1064.584 

113.334  .nn] 

LENGTH 

RATE 


RINT 

2- WAY  INTERACTION? 

rOLICY  LENGTH 
POLICY  RATE 
POLICY  RINT 
LENGTH  RATE 
LENGTH  RUT 
RATE  RINT 

3- WAY  INTERACTIONS 

POLICY  LENGTH 
POLICY  LENGTH 
POLICY  RATE 
LENCTH  RATE 


4-WAY  INTERACTIONS 
POLICY  LENGTH 
RINT 

EXPLAINED 

RESIDUAL 

TOTAL 


16 . 9  IP 
17%. 011 


24.521 

25.555 

1.327 

4.082 

.320 

.006 

8.083 

12.028 

32.083 

.058 

.117 

.350 

31.540 


16.018  1.78! 

1796.0  11  189.51- 


24.521  2.581 


2.  155 
.442 
1.361 
.  1 10 
.006 
8 . 0°3 
12.028 

3.208 
.01° 
.039 
.  120 
31.540 


.227 
.047 
.  144 
.012 
.poi 

.853 

1.26° 


5080.664 
2  122.837 
72  12.50  ] 


31  164.183  17.324  .00! 

224  9.477 

255  28.284 


256  CASES  WERE  PROCESSED. 

0  CASES  (  0  pcT)  WERE  MISSING. 


-WV/'  .N.' 


CV 


*******  ANALYSIS  0  F  Y  A  "  TANG  .  ******** 

vri'PAT  -wom'd  ypiiiglc s  per  popp 

PY  POLICY  INTERDICTION  POLICY 

length  convoy  length 

RATE  CONVOY  P.ATF. 

PINT  RF.OON  UMTS  INTERVAL 

************************************ 


OUR  OF.  OF  VARIATION 
ATM  EFFECTS 


LENGTH 


UiflU 


r mm 

SOUARES 

5031.203 


15.9  IP, 
1796.011 


mea?: 

PF  SQUARE 


SIOMF 
op  F 


6  938.534  95.720  .001 


16.918  1.931 

1796.011  205.019 


PINT 

24.521 

1 

24.521 

2.799 

EXPLAINED 

5031.203 

6 

°3R .  534 

95.720 

RESIDUAL 

2  IP  1.298 

249 

R.760 

TOTAL 

72  12.501 

255 

29.2P4 

256  CASES 
0  CASES 

WERE 

( 

PROCESSED. 

0  PCT)  WERE 

HISSING. 

i?*TT,,J7,T*T*TpTV.  tv  - , .  - 


•  v*  v  -  v^y* _:ns.i.  S'*j7vsm,s'vy'- ’y^.’vyr's 


'».  'j. 1  j  i.i  '  j  .1.' .« .■ 


********  a  n  a  l  v  o  I  s  n  p  v  a  r.  t  a  n  c  e  ******** 

YE  HP.  AT  VEHICLE1)  rT.p  POUR 

r.Y  mi. ICY  INTERDICTION  POLICY 

RATE  pn-Tnv  RATE 

RINT  pp.coy  n»’iT2  INTERVAL 

************************************* 


SOURCE  OF  VARIATION 


SUN  OF 
SOUARES 


DF 


MTT 

SCARE 


SI  C*'  IF 
OF  F 


■ 

MAIN  EFFECTS 

5014.285 

5 

1002.857 

1  10.335 

.oo] 

POLICY 

3  103.7  53 

3 

1064.584 

117.  126 

.  on  ] 

7k 

RATE 

3705.01] 

1 

1796.011 

197. 5°7 

.o°] 

& 

T>  TV*p 

24.52] 

1 

24.521 

2.698 

.  102 

2-MAY  INTERACTIONS 

16.43° 

7 

2.348 

.2.58 

.960 

POLICY  RATE 

4.082 

3 

1.361 

.  150 

.030 

POLICY  PINT 

.32° 

3 

.  1 10 

.012 

,°0P 

Kv- 

RATE  PINT 

12.028 

1 

12.028 

1.323 

.251 

3-0 AV  INTERACTIONS 

.35° 

3 

.  1 20 

.013 

.oqp 

POLICY  RATE 

RTNT  .359 

3 

.  120 

.0  13 

.008 

P 

EXPLAINED 

5031.083 

15 

335.406 

36.90  1 

.001 

fc> 

RESIDUAL 

218  1.4  18 

240 

o.08° 

| 

TOTAL 

72  12.501 

255 

28.284 

P 

2  56  CASES  ME  RE 

PROCESSED. 

ly> 

0  CASES  ( 

0  PCT)  MF.RF  MISSINC. 

********  A  v  1 

!  Y  S  I  S  n  F  V  A 

n  t  / 

\  c,  p  * 

★  *  *  * 

*  *  * 

VEHRAT 

AR  10RF.!-  VESICLES  PER 

PY  POLICY 

INTEPDICTIOU  pnLICY’ 

PATE 

CONVOY  RATE 

PINT 

RZCON  UNITS  TNTERVAL 

************ 

********** 

*  *  * 

★  *  ★  *  •k 

*  *  *  * 

★  ★  ★ 

SUM  op 

Mrv: 

SICN1F 

SOURCE  OF  VARY ‘"ION 

SOUAEES 

DF 

squahh 

r 

OF  F 

MAIN  EFFECTS 

5014.285 

5 

1002.857 

1 14.064 

.001 

ror.ICY 

3193.753 

3 

1064.584 

121.074 

.001 

RATE 

1706.01] 

1 

1796.011 

204.258 

.on] 

RINT 

24.521 

1 

24.52 1 

2.780 

.096 

EXPLAINED 

5014.285 

5 

1002.857 

1  14.054 

.00  3 

RESIDUAL 

2  198.2  16 

250 

8.793 

TOTAL 

7212.501 

255 

28.284 

256  CASES  '.'ERF.  PROCESSED. 

0  CASES  (  o  PCT)  '..’ERE  MISSING. 
IMHO  INTERDICTION  ANALYSIS 


03/09/83  17.22.00. 


PACE 


Tape  1 


25. on  25 
25.no  25 
2  5.no  2  5 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  22 
25.00  22 
25.0 o  22 
25.00  22 
25.00  22 
25.00  22 
25.no  22 
25.00  22 

22.50  25 

22.50  25 

22.50  25 

22.50  25 

22.50  25 

22.50  25 

22.50  25 

22.50  25 

22.50  22 
22.5n  22 

22.50  22 

22.50  22 

22.50  22 

22.50  22 

22.50  22 

22.50  22 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  25 
25.00  22 
25.00  22 
25.00  22 
25.00  22 
25.00  22 
25.00  22 
25.00  22 
25.00  22 

22.50  25 


,00 

4.00 

,00 

4. no 

,00 

4.oo 

no 

4.00 

no 

4.00 

no 

4.00 

,00 

4  .no 

00 

4.00 

50 

3.50 

50 

3.50 

50 

3.50 

50 

3.50 

50 

3.60 

,50 

3.60 

,50 

3.60 

,50 

3.60 

,00 

4. 00 

,00 

4.00 

,00 

4.00 

,00 

4.00 

00 

4.00 

,00 

4.00 

00 

4.00 

00 

4.00 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

00 

4.00 

00 

4  .00 

00 

4.00 

00 

4.00 

00 

4.00 

00 

4. no 

00 

4. oo 

00 

4.00 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

50 

3.60 

00 

4.00 

\"av'.  ‘'.v.1.  v 


■^  •>  _-w 


23 

24 

25 

26 
2.1 
28 

29 

30 
3] 

32 

33 

34 

35 
34 

37 

38 

39 

40 

41 

42 

43 

44 

45 
44 
47 
4P 

49 

50 
f\  1 

52 

53 

54 

55 
54 

57 

58 

59 

40 

41 

42 

43 

44 
1 
2 

3 

4 

5 
4 

7 

8 


3 

o.oo 

22.50 

25.no 

4.00 

70.35 

3 

n.OO 

22.50 

25.no 

4.00 

68.20 

3 

0.00 

22.50 

22.  so 

3.40 

69.61 

3 

0.00 

2  2.50 

22.50 

3.40 

66.22 

3 

0.00 

22.50 

22.50 

3.40 

74.57 

3 

0.00 

22.50 

22.50 

3.40 

65.30 

3 

0.00 

22 . 50 

22.50 

3.40 

66.08 

3 

0.00 

22.50 

22.50 

3.40 

69 . 14 

3 

0.00 

22.50 

22.50 

3.40 

67.36 

3 

0.00 

22.50 

22.50 

3.40 

73.99 

3 

.  10 

25.00 

25.00 

4.00 

72. 4  P 

3 

.  10 

25.00 

2  5.00 

4  .00 

71.75 

3 

.  10 

25.00 

25.00 

4.00 

68.60 

3 

.  10 

25.00 

25.00 

4.00 

77.43 

3 

.  10 

25.00 

25.00 

4.00 

77.61 

3 

.  10 

25.00 

25.00 

4  .00 

78.78 

3 

.  10 

25.00 

25.00 

4.00 

74.45 

3 

.  10 

2  5.00 

25.00 

4.00 

78.  12 

3 

.  10 

25.00 

22.50 

3.40 

73.  10 

3 

.  10 

25.00 

22.50 

3.40 

73.32 

3 

.  10 

25.oo 

22.50 

3.60 

73.06 

3 

.  10 

25.00 

22.50 

3.60 

72.  16 

3 

.  ]n 

25.00 

22.50 

3.60 

73.89 

3 

.  10 

25.00 

22.50 

3.60 

77.11 

3 

.  10 

25.00 

22.50 

3.60 

78.29 

3 

.  10 

25.00 

22.50 

3.60 

76.65 

3 

.  10 

22.50 

25.00 

4.00 

67.20 

3 

.  10 

22.50 

25.no 

4. 00 

67.38 

3 

.  10 

22.50 

25.00 

4.00 

66.33 

3 

.  10 

22.50 

25.00 

4.00 

72.81 

3 

.  10 

22.50 

25.00 

4.00 

69.57 

3 

.  10 

22.50 

25.00 

4.00 

69.  12 

3 

.  10 

22.50 

2  5.00 

4.00 

67.02 

3 

.  10 

22.50 

2  5.00 

4.00 

67.88, 

3 

.  10 

22.50 

22.50 

3.60 

68.75 

3 

.  10 

22.50 

22.50 

3.60 

75.29 

3 

.  in 

22.50 

22.50 

3.60 

7  1.84 

3 

.  10 

22 . 50 

22.50 

3.60 

70.  10 

3 

.  10 

22.50 

32.50 

3.60 

67.29 

3 

.  10 

22.50 

22.50 

3.60 

71.32 

3 

.  10 

22.50 

22.50 

3.60 

70 .  np 

3 

.  10 

22.50 

22.50 

3.60 

70.00 

4 

0.00 

25.00 

25.00 

',.00 

66.46 

4 

0.00 

25.00 

25.00 

4.00 

70.02 

4 

o.on 

25.00 

25.00 

4.00 

73.92 

4 

o.nn 

25.00 

25.00 

4.00 

72.37 

4 

0.00 

25.on 

25.00 

4.00 

72.33 

4 

O.o  o 

2  5.00 

25.00 

4.no 

76.30 

4 

o.on 

25.00 

25.00 

4.00 

72.34 

4 

o.oo 

25.00 

25.00 

4.00 

78.30 
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O.nn 

n.rn 

n. no 
O.no 
O.no 
0.00 
0.00 

o. oo 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 

.  10 
.  10 
.  10 
.  10 
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APPENDIX  D 


TESTS  OF  STATISTICAL  DISTRIBUTIONS 


SPSS  tests  for  uniform  distribution  are  presented 
in  this  Appendix.  The  uniform  distribution  used  in  SORTIE 
and  MISSILE  was  tested  using  a  Kolnogorov-Smirnov  (K-S) 
test.  The  SPSS  program  for  this  analysis  precedes  the 
results.  An  alpha  of  0.05  was  used  in  testing  a  sample  of 
576  random  variates.  A  sample  of  100  of  these  576  numbers, 
tape  12,  is  shown  following  the  SPSS  results 

Distribution  of  SORTIE  and  MISSILE 
The  K-S  test  done  tested  the  following: 


H  :  The  576  random  variates  came  from  a  Uniform 
(0,1)  distribution. 

H1 :  The  576  random  variates  did  not  come  from  a 
Uniform  (0,1)  distribution. 


Dcrit  "  D. 05, 576 


0.057 


The  underlined  portion  of  the  SPSS  printout  shows  that  the 
computed  value  was  0.0253. 

0.0253  <  0.057 

I D 1  max  <  Dcrit  for  alPha  equals  to  0.05 


Therefore,  the  null  hypothesis  cannot  be  rejected. 
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VERIFICATION  DATA 
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The  following  outputs  from  the  model  demonstrate  the 
model  is  working  as  intended.  The  first  output  shows  the 
information  used  in  developing  the  measure  of  effectiveness. 
Tape  18  shows  no  interdiction  result  on  total  number  of 
armored  vehicles,  the  finish  time  for  the  last  entity  on 
each  route,  and  the  MRD  average  finish  time.  This  tape 
shows  a  constant  number  of  armored  vehicles,  but  varying 
finish  times.  Tape  21,  representing  Policy  four,  shows  the 
variation  of  both  time  and  vehicles.  The  last  two  tapes, 
tapes  14  and  17,  show  the  effect  on  order  of  finish  and  time 
of  finish  for  no  interdiction  and  interdiction.  Tape  14 
represents  no  interdiction  and  17  is  interdiction.  The 
columns  on  tapes  14  and  17  from  left  to  right  represent 
route  number,  run  number,  policy,  number  of  armored  vehicles, 
unit  number,  and  finish  time.  Thus,  the  model  is  maintaining 
convoy  order  as  originally  intended.  The  last  three  pages 
in  this  appendix  represent  selected  print  statements  included 
at  strategic  portions  of  the  model  to  verify  the  interdiction 
subroutines.  These  show  the  sortie  and  missile  attacks 
results.  All  global  variables  are  as  explained  at  start  of 
SLAM  coded  portion  of  network.  Appendix  A.  The  'AT'  repre¬ 
sents  shorthand  for  attribute.  The  printed  values  show  that 
the  model  is  working  as  intended. 
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